(* SiF4 computation; coordinates and radii from SiF40005.exe, 23.01.2012 *)

t={{0.00000,0.00000,0.00000,0.59534},{0.97878,0.97878,0.97878,1.09996},{1.70145,1.70145,1.70145,0.15174},{2.55097,1.06684,2.55097,0.90865},{2.55097,2.55097,1.06684,0.90865},{1.06684,2.55097,2.55097,0.90865},{-0.97878,-0.97878,0.97878,1.09996},{-1.70145,-1.70145,1.70145,0.15174},{-2.55097,-2.55097,1.06684,0.90865},{-1.06684,-2.55097,2.55097,0.90865},{-2.55097,-1.06684,2.55097,0.90865},{0.97878,-0.97878,-0.97878,1.09996},{1.70145,-1.70145,-1.70145,0.15174},{1.06684,-2.55097,-2.55097,0.90865},{2.55097,-1.06684,-2.55097,0.90865},{2.55097,-2.55097,-1.06684,0.90865},{-0.97878,0.97878,-0.97878,1.09996},{-1.70145,1.70145,-1.70145,0.15174},{-2.55097,1.06684,-2.55097,0.90865},{-2.55097,2.55097,-1.06684,0.90865},{-1.06684,2.55097,-2.55097,0.90865}};

x=Table[t[[i,1]],{i,1,21}];
y=Table[t[[i,2]],{i,1,21}];
z=Table[t[[i,3]],{i,1,21}];
r=Table[t[[i,4]],{i,1,21}];
sq=Sqrt[2]/2;

(* xy plane *)

plot2=Table[Graphics[{Circle[{x[[i]],y[[i]]},r[[i]]]}],{i,1,21}];

plot2a=Graphics[{Disk[{x[[3]],y[[3]]},0.08], Disk[{x[[8]],y[[8]]},0.08], Disk[{x[[13]],y[[13]]},0.08], Disk[{x[[18]],y[[18]]},0.08], Disk[{0,0},0.08],{Thickness[0.006],Line[{{x[[8]],y[[8]]},{x[[3]],y[[3]]}}],Line[{{x[[18]],y[[18]]},{x[[13]],y[[13]]}}] }} ];

Show[plot2,plot2a,{AspectRatio -> Automatic,Axes -> True,GridLines -> Automatic, PlotRange -> {{-4.75,4.75},{-4.75,4.75}}, Frame -> True}]

plot2b=Table[Graphics3D[{Sphere[{x[[i]],y[[i]],z[[i]]},r[[i]]]}],{i,1,21}];
Show[plot2b,(*plot2a,*){AspectRatio -> Automatic}]

t1=Table[{x[[i]],y[[i]],z[[i]]},{i,1,21}];
t=t1.RotationMatrix[π/4,{0,1,0}];

plot3=Table[Graphics[{Circle[{t[[i,1]],t[[i,2]]},r[[i]]]}],{i,1,21}];

plot3a=Graphics[{Disk[{t[[3,1]],t[[3,2]]},0.08], Disk[{t[[8,1]],t[[8,2]]},0.08], Disk[{t[[13,1]],t[[13,2]]},0.08], Disk[{t[[18,1]],t[[18,2]]},0.08], Disk[{0,0},0.1],{Thickness[0.006],Line[{{t[[8,1]],t[[8,2]]},{t[[1,1]],t[[1,2]]}}],Line[{{t[[1,1]],t[[1,2]]},{t[[3,1]],t[[3,2]]}}],Line[{{t[[1,1]],t[[1,2]]},{t[[18,1]],t[[18,2]]}}],Line[{{t[[1,1]],t[[1,2]]},{t[[13,1]],t[[13,2]]}}] }} ];

Show[plot3,plot3a,{AspectRatio -> Automatic,Axes -> True,GridLines -> Automatic, PlotRange -> {{-4.75,4.75},{-4.75,4.75}}, Frame -> True}]