(* [Fe(CN)6]---- computation; coordinates and radii from FeCN6II8.exe, 02.10.2012/15.02.2013 *)

t={{0.00000,0.00000,0.00000,2.49901},{2.33848,0.00000,0.00000,1.69653},{4.29533,0.00000,0.00000,0.26032},{5.40207,0.00000,0.00000,0.84643},{5.40207,1.89008,0.00000,1.89008},{5.40207,0.00000,1.89008,1.89008},{5.40207,-1.89008,0.00000,1.89008},{5.40207,0.00000,-1.89008,1.89008},{6.47226,0.00000,0.00000,0.22376},{7.79852,0.00000,0.00000,1.10250},{0.00000,2.33848,0.00000,1.69653},{0.00000,4.29533,0.00000,0.26032},
{0.00000,5.40207,0.00000,0.84643},{1.89008,5.40207,0.00000,1.89008},{0.00000,5.40207,1.89008,1.89008},
{-1.89008,5.40207,0.00000,1.89008},{0.00000,5.40207,-1.89008,1.89008},{0.00000,6.47226,0.00000,0.22376},
{0.00000,7.79852,0.00000,1.10250},{-2.33848,0.00000,0.00000,1.69653},{-4.29533,0.00000,0.00000,0.26032},{-5.40207,0.00000,0.00000,0.84643},{-5.40207,-1.89008,0.00000,1.89008},{-5.40207,0.00000,-1.89008,1.89008},{-5.40207,1.89008,0.00000,1.89008},{-5.40207,0.00000,1.89008,1.89008},{-6.47226,0.00000,0.00000,0.22376},{-7.79852,0.00000,0.00000,1.10250}};

x=Table[t[[i,1]],{i,1,28}];
y=Table[t[[i,2]],{i,1,28}];
z=Table[t[[i,3]],{i,1,28}];
r=Table[t[[i,4]],{i,1,28}];

(* one of the three sigmah planes *)
plot2=Table[Graphics[{Circle[{x[[i]],y[[i]]},r[[i]]]}],{i,1,28}];

plot21=Table[Graphics[{Circle[{x[[i]],-y[[i]]},r[[i]]]}],{i,11,19}];
plot22=Graphics[{Circle[{x[[1]],y[[1]]},0.05],Circle[{x[[1]],y[[1]]},0.25],Circle[{x[[1]],y[[1]]},0.636],
Disk[{x[[3]],y[[3]]},0.1], Disk[{x[[9]],y[[9]]},0.1], 
Disk[{x[[12]],y[[12]]},0.1], Disk[{x[[12]],-y[[12]]},0.1],Disk[{x[[18]],y[[18]]},0.1],Disk[{x[[18]],-y[[18]]},0.1], Disk[{x[[21]],y[[21]]},0.1], Disk[{x[[27]],y[[27]]},0.1], 
Disk[{0,0},0.16] } ];

Show[plot2,plot21,plot22,{AspectRatio -> Automatic,Axes -> True,GridLines -> Automatic, PlotRange -> {{-10,10},{-10,10}}, Frame -> True}]

plot3=Table[Graphics3D[{Opacity[0.6],Sphere[{x[[i]],y[[i]],z[[i]]},r[[i]]]}],{i,1,28}];
plot31=Table[Graphics3D[{Opacity[0.3],Lighter[Blue,0.3],Sphere[{x[[i]],y[[i]],z[[i]]},r[[i]]]}],{i,5,8}];
plot32=Table[Graphics3D[{Opacity[0.3],Lighter[Blue,0.3],Sphere[{x[[i]],y[[i]],z[[i]]},r[[i]]]}],{i,14,17}];
plot33=Table[Graphics3D[{Opacity[0.3],Lighter[Blue,0.3],Sphere[{x[[i]],y[[i]],z[[i]]},r[[i]]]}],{i,23,26}];

plot4=Table[Graphics3D[{Opacity[0.6],Sphere[{x[[i+10]],-y[[i+10]],z[[i+10]]},r[[i+10]]]}],{i,1,9}];
plot41=Table[Graphics3D[{Opacity[0.3],Lighter[Blue,0.3],Sphere[{x[[i+10]],-y[[i+10]],z[[i+10]]},r[[i+10]]]}],{i,4,7}];

plot5 = Graphics3D[{Darker[Red,1.0],Sphere[{x[[1]],y[[1]],z[[1]]},0.05],Sphere[{x[[1]],y[[1]],z[[1]]},0.25],Sphere[{x[[1]],y[[1]],z[[1]]},0.636]}];

plot6=Table[Graphics3D[{Opacity[0.6],{Lighter[Yellow],Sphere[{x[[i+10]],-z[[i+10]],y[[i+10]]},r[[i+10]]]}}],{i,1,9}];
plot61=Table[Graphics3D[{Opacity[0.3],{Lighter[Blue,0.3],Sphere[{x[[i+10]],-z[[i+10]],y[[i+10]]},r[[i+10]]]}}],{i,4,7}];

plot7=Table[Graphics3D[{Opacity[0.6],{Lighter[Yellow],Sphere[{x[[i+10]],-z[[i+10]],-y[[i+10]]},r[[i+10]]]}}],{i,1,9}];
plot71=Table[Graphics3D[{Opacity[0.3],{Lighter[Blue,0.3],Sphere[{x[[i+10]],-z[[i+10]],-y[[i+10]]},r[[i+10]]]}}],{i,4,7}];

plot8=Graphics3D[{Opacity[1.0],{Darker[Brown,0.5],Sphere[{x[[3]],y[[3]],z[[3]]},r[[3]]], Sphere[{x[[12]],y[[12]],z[[12]]},r[[12]]], Sphere[{x[[12]],-y[[12]],z[[12]]},r[[12]]], Sphere[{x[[12]],-z[[12]],y[[12]]},r[[12]]], Sphere[{x[[12]],-z[[12]],-y[[12]]},r[[12]]], Sphere[{x[[21]],y[[21]],z[[21]]},r[[21]]]},{Lighter[Blue],Sphere[{x[[9]],y[[9]],z[[9]]},r[[9]]], Sphere[{x[[18]],y[[18]],z[[18]]},r[[18]]],Sphere[{x[[18]],-y[[18]],z[[18]]},r[[18]]], Sphere[{x[[18]],-z[[18]],y[[18]]},r[[18]]], Sphere[{x[[18]],-z[[18]],-y[[18]]},r[[18]]],  Sphere[{x[[27]],y[[27]],z[[27]]},r[[27]]] }} ];

plot9=Graphics3D[{Opacity[1.0],{Darker[Brown,0.5],Tube[{{0,0,0},{x[[3]],y[[3]],z[[3]]}},0.15], Tube[{{0,0,0},{x[[12]],y[[12]],z[[12]]}},0.15], Tube[{{0,0,0},{x[[12]],-y[[12]],z[[12]]}},0.15], Tube[{{0,0,0},{x[[12]],-z[[12]],y[[12]]}},0.15],Tube[{{0,0,0},{x[[12]],-z[[12]],-y[[12]]}},0.15], Tube[{{0,0,0},{x[[21]],y[[21]],z[[21]]}},0.15]},{Darker[Blue],Tube[{{x[[3]],y[[3]],z[[3]]},{x[[9]],y[[9]],z[[9]]}},0.15], Tube[{{x[[12]],y[[12]],z[[12]]},{x[[18]],y[[18]],z[[18]]}},0.15],Tube[{{x[[12]],-y[[12]],z[[12]]},{x[[18]],-y[[18]],z[[18]]}},0.15], Tube[{{x[[12]],-z[[12]],y[[12]]},{x[[18]],-z[[18]],y[[18]]}},0.15], Tube[{{x[[12]],-z[[12]],-y[[12]]},{x[[18]],-z[[18]],-y[[18]]}},0.15], Tube[{{x[[21]],y[[21]],z[[21]]},{x[[27]],y[[27]],z[[27]]}},0.15] }} ];
plot10=Graphics3D[{Opacity[0.8],Lighter[Magenta,0.8],Sphere[{x[[1]],y[[1]],z[[1]]},r[[1]]]}];

Show[plot10,plot5,plot3,plot31,plot32,plot33,plot4,plot41,plot6,plot61,plot7,plot71,plot8,plot9,{AspectRatio -> Automatic,Boxed->False, SphericalRegion->True}]