(* [ HCCH computation;  
        coordinates and radii from H2C20001 . exe,  
        06.02  .2012 *) 
 (* 
      pi clouds off axis and towards N .  
          pi optimization produces nonsense! *) 

t = {{0.00000, 0.00000, 
            0.00000, 0.25080}, {2.35323, 0.00000, 0.00000, 
            0.25080}, {-1.19937, 0.00000, 0.00000, 0.94858}, {1.17662, 
            0.00000, 0.00000, 0.92582}, {1.17662, 2.24758, 0.00000, 
            2.24758}, {1.17662, 0.00000, 2.24758, 
            2.24758}, {1.17662, -2.24758, 0.00000, 2.24758}, {1.17662, 
            0.00000, -2.24758, 2.24758}, {3.55261, 0.00000, 0.00000, 
            0.94858}};

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

plot2 = 
      Graphics[{Circle[{x[[1]], y[[1]]}, 
            r[[1]]], 
Circle[{x[[2]], y[[2]]}, 
            r[[2]]], Circle[{x[[9]], y[[9]]}, r[[9]]], 
          Circle[{x[[3]], y[[3]]}, r[[3]]], 
          Circle[{x[[4]], y[[4]]}, r[[4]]], 
          Circle[{x[[5]], y[[5]]}, r[[5]]], 
          Circle[{x[[6]], y[[6]]}, r[[6]]], 
          Circle[{x[[7]], y[[7]]}, r[[7]]], 
          Circle[{x[[8]], y[[8]]}, r[[8]]], 
          Disk[{x[[1]], y[[1]]}, 0.08],  
          Disk[{x[[2]], y[[2]]}, 0.08],  
          Disk[{-1.6829, 0}, 0.08],  
          Disk[{x[[2]] + 1.6829, 0}, 
            0.08]} ]

Show[
      plot2, {AspectRatio  ->  Automatic, Axes  ->  True, 
        GridLines  ->  Automatic,  
        PlotRange  ->  {{-3, 5}, {-4, 4}},  
        Frame  ->  True}]

  
plot3 = 
      Graphics[{Circle[{y[[1]], z[[1]]}, 
            r[[1]]], 
Circle[{sq*y[[2]], 
              z[[2]]}, r[[2]]], 
          Circle[{sq*y[[3]], z[[3]]}, r[[3]]], 
          Circle[{sq*y[[4]], z[[4]]}, r[[4]]], 
          Circle[{sq*y[[5]], z[[5]]}, r[[5]]], 
          Circle[{sq*y[[6]], z[[6]]}, r[[6]]], 
          Circle[{sq*y[[7]], z[[7]]}, r[[7]]], 
          Circle[{sq*y[[8]], z[[8]]}, r[[8]]], 
          Circle[{sq*y[[9]], z[[9]]}, r[[9]]], 
          Disk[{sq*x[[1]], z[[1]]}, 
            0.08]}]

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