> restart:
> with(plots):

>ShowForce:=proc(FileName::string,AtomSize,AtomColor,ArrowLength,ArrowShape)
	local data,Position,i,Force,D0:
	data:=readdata(FileName,6):
	Position:=[]:
	for i from 1 to nops(data) do
		Position:=[op(Position),[data[i][1],data[i][2],data[i][3]]]:
		Force[i]:=arrow(<data[i][1],data[i][2],data[i][3]>;,
		<data[i][4]*ArrowLength,data[i][5]*ArrowLength,data[i][6]*ArrowLength>,
		shape=ArrowShape,color=AtomColor,width=0.25):
	od:
	D0:=polyhedraplot(Position,polyscale=AtomSize,polytype=icosahedron,scaling=constrained,color=AtomColor,axes=BOXED):
	RETURN([D0,seq(Force[i],i=1..nops(data))]):
 end:

># Usage: display(ShowForce("Filename",Atomsize(int),AtomColor("red","blue","green",...),ArrowLength(float),arrow));