(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 9.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[     74454,       1283]
NotebookOptionsPosition[     72418,       1220]
NotebookOutlinePosition[     73015,       1242]
CellTagsIndexPosition[     72972,       1239]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell["\<\
(*NH2COH Formamide, with \"orthogonal\" \[Sigma],\[Pi] LPs, 13.01.2013, not \
final optimization: slight overlaps *)\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{
  3.5668390630515456`*^9, {3.5668394752510695`*^9, 3.5668394888386936`*^9}, {
   3.566840005402401*^9, 3.5668400293172426`*^9}, {3.566840906678384*^9, 
   3.566840923027212*^9}, {3.5668417775723133`*^9, 3.5668417895999346`*^9}, {
   3.566841995894697*^9, 3.5668420196535387`*^9}, {3.566842057031204*^9, 
   3.5668420654552193`*^9}, {3.5668422519691467`*^9, 3.5668422563371544`*^9}, 
   3.5668422887072115`*^9, {3.566842337504097*^9, 3.566842360326937*^9}, {
   3.5668424117758274`*^9, 3.5668424218690453`*^9}, {3.566842627885008*^9, 
   3.566842684060706*^9}, {3.566842771030859*^9, 3.5668427769432697`*^9}, {
   3.5668428297025623`*^9, 3.5668428906206694`*^9}, {3.5668429241763277`*^9, 
   3.566842930135538*^9}, {3.5668871695283947`*^9, 3.5668871699963956`*^9}, {
   3.5668872245808916`*^9, 3.566887225407693*^9}, {3.5668875348030367`*^9, 
   3.566887535037037*^9}, {3.566887597327946*^9, 3.566887637981618*^9}, 
   3.5668876998201265`*^9, {3.5668877336253853`*^9, 3.5668877352165885`*^9}, {
   3.566887814495928*^9, 3.5668878152135286`*^9}, {3.566887857926404*^9, 
   3.566887961448186*^9}, 3.5668880086070685`*^9, {3.5668880951092205`*^9, 
   3.5668881570569296`*^9}, {3.5668882853983545`*^9, 
   3.5668884070317683`*^9}, {3.566888509258748*^9, 3.5668885097423487`*^9}, {
   3.566888582719277*^9, 3.566888583795679*^9}, {3.566888668176227*^9, 
   3.566888671935834*^9}, {3.5668887042746906`*^9, 3.566888705070292*^9}, {
   3.566888771822809*^9, 3.566888808857274*^9}, {3.5668889463559155`*^9, 
   3.566888956839134*^9}, {3.5669685950612297`*^9, 3.5669685991952367`*^9}, {
   3.566968722700654*^9, 3.5669687698127365`*^9}, {3.5669688330864477`*^9, 
   3.566968848015674*^9}, {3.5669689606010714`*^9, 3.566968973471094*^9}, {
   3.5669881655485907`*^9, 3.5669881659853916`*^9}, {3.566996986024498*^9, 
   3.5669969910321074`*^9}, {3.5669973015818524`*^9, 
   3.5669973626247597`*^9}, {3.566997996329073*^9, 3.566998005049488*^9}, {
   3.56705755136038*^9, 3.5670575543243856`*^9}, {3.5670590414593973`*^9, 
   3.5670590456714053`*^9}, {3.567061276787324*^9, 3.567061372633892*^9}, {
   3.567061756144966*^9, 3.5670617569249673`*^9}, {3.567061904329626*^9, 
   3.567061915546046*^9}, {3.567062309945938*^9, 3.5670623507556105`*^9}, {
   3.567062381019663*^9, 3.5670625124810944`*^9}, {3.5670638527106485`*^9, 
   3.5670638629598665`*^9}, 3.5670687979573345`*^9, {3.567078692613513*^9, 
   3.56707870788594*^9}}],

Cell["\<\
Clear[ZC,ZN,ZO,sig1,sig2,sig3,sig4,sig5,sig6,sig7,sig8,sig9,k1,k2,k3,k4k5,k6,\
k7,k8,k9,k10,k11,nc,w1,w2,w3,pnh,pch,vee,vne,vnn,
xn,yn,zn,xc,yc,zc,c1,c2,c3,cs1,ss1,cs2,ss2,cs3,ss3,pi,i,j,d,t];
Clear[Cr,Nr,Oa,NH,CH,CO,CN,OL,NM,OC,OM,w1,w2,w3,pch,pnh];
ZC=6.0; ZN=7.0; ZO=8.0; sig1=0.3; sig2=0.3; sig3=0.3; sig4=0.3; sig5=0.3; \
sig6=0.3;
sig7=0.3; sig8=0.3; sig9=0.3;
nc=15; nn=6; pi=0.0;
k1=1.02; k2=1.02; k3=1.02; k4=1.2; k5=1.25; k6=1.56; k7=1.26; k8=1.4; k9=1.1; \
k10=1.1; k11=1.0;  
bohr=0.529177; rad=57.29578;

(* 16 Variables: Cr,Nr,Or,NH,CH,CN,CO,OL,NM,OC,OM,pnh,pch,w1,w2,w3 *)\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{
  3.5668390630515456`*^9, {3.5668394752510695`*^9, 3.5668394888386936`*^9}, {
   3.566840005402401*^9, 3.5668400293172426`*^9}, {3.566840906678384*^9, 
   3.566840923027212*^9}, {3.5668417775723133`*^9, 3.5668417895999346`*^9}, {
   3.566841995894697*^9, 3.5668420196535387`*^9}, {3.566842057031204*^9, 
   3.5668420654552193`*^9}, {3.5668422519691467`*^9, 3.5668422563371544`*^9}, 
   3.5668422887072115`*^9, {3.566842337504097*^9, 3.566842360326937*^9}, {
   3.5668424117758274`*^9, 3.5668424218690453`*^9}, {3.566842627885008*^9, 
   3.566842684060706*^9}, {3.566842771030859*^9, 3.5668427769432697`*^9}, {
   3.5668428297025623`*^9, 3.5668428906206694`*^9}, {3.5668429241763277`*^9, 
   3.566842930135538*^9}, {3.5668871695283947`*^9, 3.5668871699963956`*^9}, {
   3.5668872245808916`*^9, 3.566887225407693*^9}, {3.5668875348030367`*^9, 
   3.566887535037037*^9}, {3.566887597327946*^9, 3.566887637981618*^9}, 
   3.5668876998201265`*^9, {3.5668877336253853`*^9, 3.5668877352165885`*^9}, {
   3.566887814495928*^9, 3.5668878152135286`*^9}, {3.566887857926404*^9, 
   3.566887961448186*^9}, 3.5668880086070685`*^9, {3.5668880951092205`*^9, 
   3.5668881570569296`*^9}, {3.5668882853983545`*^9, 
   3.5668884070317683`*^9}, {3.566888509258748*^9, 3.5668885097423487`*^9}, {
   3.566888582719277*^9, 3.566888583795679*^9}, {3.566888668176227*^9, 
   3.566888671935834*^9}, {3.5668887042746906`*^9, 3.566888705070292*^9}, {
   3.566888771822809*^9, 3.566888808857274*^9}, {3.5668889463559155`*^9, 
   3.566888956839134*^9}, {3.5669685950612297`*^9, 3.5669685991952367`*^9}, {
   3.566968722700654*^9, 3.5669687698127365`*^9}, {3.5669688330864477`*^9, 
   3.566968848015674*^9}, {3.5669689606010714`*^9, 3.566968973471094*^9}, {
   3.5669881655485907`*^9, 3.5669881659853916`*^9}, {3.566996986024498*^9, 
   3.5669969910321074`*^9}, {3.5669973015818524`*^9, 
   3.5669973626247597`*^9}, {3.566997996329073*^9, 3.566998005049488*^9}, {
   3.56705755136038*^9, 3.5670575543243856`*^9}, {3.5670590414593973`*^9, 
   3.5670590456714053`*^9}, {3.567061276787324*^9, 3.567061372633892*^9}, {
   3.567061756144966*^9, 3.5670617569249673`*^9}, {3.567061904329626*^9, 
   3.567061915546046*^9}, {3.567062309945938*^9, 3.5670623507556105`*^9}, {
   3.567062381019663*^9, 3.5670625124810944`*^9}, {3.5670638527106485`*^9, 
   3.5670638629598665`*^9}, 3.5670687979573345`*^9, {3.567069322009055*^9, 
   3.5670693806495576`*^9}, {3.567069419290826*^9, 3.567069451957283*^9}, {
   3.5670694881649466`*^9, 3.5670695381162343`*^9}, {3.56708248268017*^9, 
   3.567082487516179*^9}, {3.567085712191843*^9, 3.56708586995492*^9}, {
   3.5670863390009437`*^9, 3.5670863409197474`*^9}, {3.5670863724006023`*^9, 
   3.5670863736018047`*^9}, {3.567086415363078*^9, 3.56708641629908*^9}, {
   3.5670864537547455`*^9, 3.567086586245778*^9}, {3.56708706580622*^9, 
   3.567087066742222*^9}, {3.5670871956452484`*^9, 3.567087197376851*^9}, {
   3.567087313846656*^9, 3.5670873147670574`*^9}, {3.567087346825114*^9, 
   3.5670873478703156`*^9}, {3.5670873954503994`*^9, 3.567087396292801*^9}, {
   3.5670874292868586`*^9, 3.5670874532017007`*^9}, {3.567087493090971*^9, 
   3.5670874935745716`*^9}, {3.5670979031336823`*^9, 3.5670979250517206`*^9}}],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{
   RowBox[{"Cr", "=", "0.26"}], ";", " ", 
   RowBox[{"Nr", "=", "0.22"}], ";", " ", 
   RowBox[{"Oa", "=", "0.186"}], ";", " ", 
   RowBox[{"NH", "=", "1.2"}], ";", " ", 
   RowBox[{"CH", "=", "1.25"}], ";", " ", 
   RowBox[{"CO", "=", "0.928"}], ";", " ", 
   RowBox[{"CN", "=", "1.04"}], ";", " ", 
   RowBox[{"OL", "=", "1.0"}], ";", "\[IndentingNewLine]", 
   RowBox[{"NM", "=", "1.1"}], ";", " ", 
   RowBox[{"OC", "=", "1.4"}], ";", " ", 
   RowBox[{"OM", "=", "1.05"}], ";", " ", 
   RowBox[{"w1", "=", "1.9"}], ";", " ", 
   RowBox[{"w2", "=", "1.0"}], ";", " ", 
   RowBox[{"w3", "=", "1.0"}], ";", " ", 
   RowBox[{"pch", "=", "0.58"}], ";", 
   RowBox[{"pnh", "=", "0.5"}], ";"}], " ", "*)"}]], "Input",
 CellChangeTimes->{{3.5670689879656677`*^9, 3.5670691959452333`*^9}, 
   3.5670698299303465`*^9, {3.567077287748646*^9, 3.5670773058290777`*^9}, {
   3.567077396527637*^9, 3.567077402346447*^9}, {3.5670803249623804`*^9, 
   3.5670804001545124`*^9}, 3.5670804667198296`*^9, 3.5670805071239004`*^9, {
   3.5670810452624454`*^9, 3.567081083794513*^9}, {3.5670811174593725`*^9, 
   3.5670811463974233`*^9}, {3.5670811875190954`*^9, 
   3.5670812332895756`*^9}, {3.5670812684832373`*^9, 3.567081269356839*^9}, {
   3.5670852872782965`*^9, 3.567085293970708*^9}}],

Cell["\<\
(* C He shell *)
Ekin = 2.25*k1/Cr^2;
vee=3.0*sig1/Cr;
vne=-3.0*ZC/Cr;

(* N He shell *)
Ekin = Ekin+2.25*k2/Nr^2;
vee=vee+3.0*sig2/Nr;
vne=vne-3.0*ZN/Nr;

(* O He shell *)
Ekin = Ekin+2.25*k3/Oa^2;
vee=vee+3.0*sig3/Oa;
vne=vne-3.0*ZO/Oa;\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5668390630515456`*^9, 3.5668390682151546`*^9}, {
  3.566971938334302*^9, 3.566971944714713*^9}, {3.566972152507078*^9, 
  3.5669721651899004`*^9}, {3.5669970287685733`*^9, 3.5669972828462195`*^9}, {
  3.566997377709986*^9, 3.5669974183012576`*^9}, {3.5669979768290386`*^9, 
  3.566997987171857*^9}, {3.566998792414071*^9, 3.5669988056272945`*^9}, {
  3.567063874659887*^9, 3.567063885189905*^9}, {3.5670698395087633`*^9, 
  3.567069848260379*^9}}],

Cell["\<\
(* bonding and lone pairs *)
Ekin = Ekin + \
2.25*(2*k4/NH^2+k5/CN^2+k6/CO^2+k7/CH^2+k8/NM^2+k9/OL^2+k10/OM^2+k11/OC^2);

vee=vee+3.0*(2*sig4/NH+sig5/CN+sig6/CO+sig7/CH+sig9/OL);\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5668390630515456`*^9, 3.566839078760773*^9}, {
   3.566839318626795*^9, 3.5668393238528037`*^9}, {3.566839506201524*^9, 
   3.5668395108503323`*^9}, {3.566839541207986*^9, 3.566839568664034*^9}, {
   3.566839684650237*^9, 3.56683970887708*^9}, {3.5668397799352045`*^9, 
   3.5668398201988754`*^9}, {3.566839877872177*^9, 3.5668398909606*^9}, {
   3.5668404106911125`*^9, 3.566840417539525*^9}, {3.5668404815464373`*^9, 
   3.566840691460406*^9}, {3.566840732519678*^9, 3.5668408608923035`*^9}, {
   3.566840930000425*^9, 3.5668409412168446`*^9}, {3.566841374686006*^9, 
   3.566841375388007*^9}, {3.5668421835006266`*^9, 3.5668421845770283`*^9}, {
   3.566887398037596*^9, 3.566887402078003*^9}, {3.5668874322328568`*^9, 
   3.5668874455864797`*^9}, {3.566887493166563*^9, 3.5668874939309645`*^9}, {
   3.5668882488786907`*^9, 3.5668882511094947`*^9}, {3.56688847075788*^9, 
   3.566888479322295*^9}, {3.5668885903476906`*^9, 3.5668885954956994`*^9}, {
   3.566971858103361*^9, 3.566971918927868*^9}, {3.566987354128765*^9, 
   3.5669873624903803`*^9}, {3.566987488164201*^9, 3.566987535619484*^9}, {
   3.566987795843541*^9, 3.5669878015843515`*^9}, {3.566987856168847*^9, 
   3.56698786325126*^9}, {3.566988264733565*^9, 3.566988269819174*^9}, {
   3.566988366133743*^9, 3.5669883719525533`*^9}, {3.5669885431160536`*^9, 
   3.566988549309265*^9}, {3.566991946182942*^9, 3.5669919483981457`*^9}, 
   3.5669919866182127`*^9, {3.5669920226542764`*^9, 3.5669920505783253`*^9}, {
   3.566992103743219*^9, 3.5669921079396257`*^9}, 3.5669921962201815`*^9, {
   3.5669922303998413`*^9, 3.5669922340658474`*^9}, {3.566992291177548*^9, 
   3.5669922926751504`*^9}, {3.566992333890423*^9, 3.5669923348264246`*^9}, 
   3.566992409675356*^9, {3.566997472043352*^9, 3.566997587171554*^9}, {
   3.566997674297707*^9, 3.5669977903463106`*^9}, {3.566997868580448*^9, 
   3.5669978780652647`*^9}, {3.566998052270771*^9, 3.566998071786405*^9}, {
   3.5669981665409718`*^9, 3.5669981730305834`*^9}, {3.56699824522751*^9, 
   3.566998326394452*^9}, {3.5669983623837156`*^9, 3.5669984017425847`*^9}, {
   3.566998470788306*^9, 3.5669984830811276`*^9}, {3.566998534405218*^9, 
   3.5669985995821323`*^9}, {3.5669988153461113`*^9, 3.566998936995125*^9}, 
   3.566998967134378*^9, 3.5669990224832754`*^9, {3.5669990973790064`*^9, 
   3.56699912179305*^9}, {3.5669993047345705`*^9, 3.5669993241410046`*^9}, {
   3.566999416274767*^9, 3.566999426742385*^9}, {3.566999497301309*^9, 
   3.5669995032137194`*^9}, {3.566999549733001*^9, 3.566999551433404*^9}, {
   3.5670567651502*^9, 3.567056829110312*^9}, {3.5670569077344503`*^9, 
   3.5670570287126627`*^9}, {3.56705707303234*^9, 3.5670571842137356`*^9}, {
   3.5670572582982655`*^9, 3.567057292103525*^9}, {3.56705732891959*^9, 
   3.5670574906294737`*^9}, {3.5670578301328697`*^9, 3.56705795512029*^9}, {
   3.5670580197668033`*^9, 3.567058030437222*^9}, {3.5670580810281105`*^9, 
   3.5670581611654515`*^9}, {3.5670581932547083`*^9, 3.567058240616391*^9}, {
   3.5670582801312604`*^9, 3.5670582990540934`*^9}, {3.567058343560972*^9, 
   3.5670584018894744`*^9}, {3.56705845049916*^9, 3.5670585083128614`*^9}, {
   3.567058578731385*^9, 3.567058593645011*^9}, {3.5670586348134837`*^9, 
   3.5670586577143235`*^9}, {3.567058689007979*^9, 3.5670587629209085`*^9}, {
   3.567058832777831*^9, 3.5670588394702425`*^9}, {3.567058955628047*^9, 
   3.5670589616496572`*^9}, {3.5670592284569263`*^9, 3.567059354442747*^9}, {
   3.5670593898548098`*^9, 3.567059560456709*^9}, {3.5670596222796183`*^9, 
   3.5670596362884426`*^9}, {3.5670596726677065`*^9, 3.567059775347087*^9}, {
   3.567059838449198*^9, 3.5670598453756094`*^9}, {3.5670598764976645`*^9, 
   3.5670598805380716`*^9}, 3.567059911254525*^9, {3.5670599744502363`*^9, 
   3.567060007818695*^9}, {3.567060059127185*^9, 3.5670600818564253`*^9}, {
   3.5670601143824825`*^9, 3.567060314593234*^9}, {3.5670607181347427`*^9, 
   3.5670607311607656`*^9}, {3.567060785573661*^9, 3.567061232998047*^9}, {
   3.5670626282800975`*^9, 3.56706267537658*^9}, 3.567062795371991*^9, {
   3.567069851801585*^9, 3.567069879460434*^9}, {3.567070767429593*^9, 
   3.5670707823588195`*^9}, {3.5670805354535503`*^9, 
   3.5670806312689185`*^9}, {3.567083393113369*^9, 3.5670834586022844`*^9}, {
   3.567083501299559*^9, 3.5670835804696984`*^9}, {3.5670836156477604`*^9, 
   3.5670836215757704`*^9}, {3.567083656504232*^9, 3.5670836570034328`*^9}, 
   3.5670837376399746`*^9}],

Cell["\<\
(*  1  2  3  4  5  6  7  8  9  10 11 12 13 14 15  *)
(*  C  N  O  Hn Hn Hc B  B  L2 L1 L1 Bp Bp L1 L1  *)  
rr={Cr,Nr,Oa,NH,NH,CH,CN,CO,OL,NM,NM,OC,OC,OM,OM};
oc={-2,-2,-2,-2,-2,-2,-2,-2,-2,-1,-1,-1,-1,-1,-1};
ch={ZC,ZN,ZO, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0};
csw1=Cos[w1/2];  ssw1=Sin[w1/2];
csw2=Cos[w2];    ssw2=Sin[w2];
csw3=Cos[w3];    ssw3=Sin[w3];\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5668390630515456`*^9, 3.566839078760773*^9}, {
   3.566839318626795*^9, 3.5668393238528037`*^9}, {3.566839506201524*^9, 
   3.5668395108503323`*^9}, {3.566839541207986*^9, 3.566839568664034*^9}, {
   3.566839684650237*^9, 3.56683970887708*^9}, {3.5668397799352045`*^9, 
   3.5668398201988754`*^9}, {3.566839877872177*^9, 3.5668398909606*^9}, {
   3.5668404106911125`*^9, 3.566840417539525*^9}, {3.5668404815464373`*^9, 
   3.566840691460406*^9}, {3.566840732519678*^9, 3.5668408608923035`*^9}, {
   3.566840930000425*^9, 3.5668409412168446`*^9}, {3.566841374686006*^9, 
   3.566841375388007*^9}, {3.5668421835006266`*^9, 3.5668421845770283`*^9}, {
   3.566887398037596*^9, 3.566887402078003*^9}, {3.5668874322328568`*^9, 
   3.5668874455864797`*^9}, {3.566887493166563*^9, 3.5668874939309645`*^9}, {
   3.5668882488786907`*^9, 3.5668882511094947`*^9}, {3.56688847075788*^9, 
   3.566888479322295*^9}, {3.5668885903476906`*^9, 3.5668885954956994`*^9}, {
   3.566971858103361*^9, 3.566971918927868*^9}, {3.566987354128765*^9, 
   3.5669873624903803`*^9}, {3.566987488164201*^9, 3.566987535619484*^9}, {
   3.566987795843541*^9, 3.5669878015843515`*^9}, {3.566987856168847*^9, 
   3.56698786325126*^9}, {3.566988264733565*^9, 3.566988269819174*^9}, {
   3.566988366133743*^9, 3.5669883719525533`*^9}, {3.5669885431160536`*^9, 
   3.566988549309265*^9}, {3.566991946182942*^9, 3.5669919483981457`*^9}, 
   3.5669919866182127`*^9, {3.5669920226542764`*^9, 3.5669920505783253`*^9}, {
   3.566992103743219*^9, 3.5669921079396257`*^9}, 3.5669921962201815`*^9, {
   3.5669922303998413`*^9, 3.5669922340658474`*^9}, {3.566992291177548*^9, 
   3.5669922926751504`*^9}, {3.566992333890423*^9, 3.5669923348264246`*^9}, 
   3.566992409675356*^9, {3.566997472043352*^9, 3.566997587171554*^9}, {
   3.566997674297707*^9, 3.5669977903463106`*^9}, {3.566997868580448*^9, 
   3.5669978780652647`*^9}, {3.566998052270771*^9, 3.566998071786405*^9}, {
   3.5669981665409718`*^9, 3.5669981730305834`*^9}, {3.56699824522751*^9, 
   3.566998326394452*^9}, {3.5669983623837156`*^9, 3.5669984017425847`*^9}, {
   3.566998470788306*^9, 3.5669984830811276`*^9}, {3.566998534405218*^9, 
   3.5669985995821323`*^9}, {3.5669988153461113`*^9, 3.566998936995125*^9}, 
   3.566998967134378*^9, 3.5669990224832754`*^9, {3.5669990973790064`*^9, 
   3.56699912179305*^9}, {3.5669993047345705`*^9, 3.5669993241410046`*^9}, {
   3.566999416274767*^9, 3.566999426742385*^9}, {3.566999497301309*^9, 
   3.5669995032137194`*^9}, {3.566999549733001*^9, 3.566999551433404*^9}, {
   3.5670567651502*^9, 3.567056829110312*^9}, {3.5670569077344503`*^9, 
   3.5670570287126627`*^9}, {3.56705707303234*^9, 3.5670571842137356`*^9}, {
   3.5670572582982655`*^9, 3.567057292103525*^9}, {3.56705732891959*^9, 
   3.5670574906294737`*^9}, {3.5670578301328697`*^9, 3.56705795512029*^9}, {
   3.5670580197668033`*^9, 3.567058030437222*^9}, {3.5670580810281105`*^9, 
   3.5670581611654515`*^9}, {3.5670581932547083`*^9, 3.567058240616391*^9}, {
   3.5670582801312604`*^9, 3.5670582990540934`*^9}, {3.567058343560972*^9, 
   3.5670584018894744`*^9}, {3.56705845049916*^9, 3.5670585083128614`*^9}, {
   3.567058578731385*^9, 3.567058593645011*^9}, {3.5670586348134837`*^9, 
   3.5670586577143235`*^9}, {3.567058689007979*^9, 3.5670587629209085`*^9}, {
   3.567058832777831*^9, 3.5670588394702425`*^9}, {3.567058955628047*^9, 
   3.5670589616496572`*^9}, {3.5670592284569263`*^9, 3.567059354442747*^9}, {
   3.5670593898548098`*^9, 3.567059560456709*^9}, {3.5670596222796183`*^9, 
   3.5670596362884426`*^9}, {3.5670596726677065`*^9, 3.567059775347087*^9}, {
   3.567059838449198*^9, 3.5670598453756094`*^9}, {3.5670598764976645`*^9, 
   3.5670598805380716`*^9}, 3.567059911254525*^9, {3.5670599744502363`*^9, 
   3.567060007818695*^9}, {3.567060059127185*^9, 3.5670600818564253`*^9}, {
   3.5670601143824825`*^9, 3.567060314593234*^9}, {3.5670607181347427`*^9, 
   3.5670607311607656`*^9}, {3.567060785573661*^9, 3.567061232998047*^9}, {
   3.5670626282800975`*^9, 3.56706267537658*^9}, 3.567062795371991*^9, {
   3.567069851801585*^9, 3.567069879460434*^9}, {3.567070767429593*^9, 
   3.5670707823588195`*^9}, {3.5670805354535503`*^9, 
   3.5670806312689185`*^9}, {3.567083393113369*^9, 3.5670834586022844`*^9}, {
   3.567083501299559*^9, 3.5670835804696984`*^9}, {3.5670836156477604`*^9, 
   3.5670836215757704`*^9}, {3.567083656504232*^9, 3.5670836570034328`*^9}, 
   3.5670837376399746`*^9}],

Cell["\<\
dcnx=Cr+Nr+2*CN;             dcox=(Cr+Oa+2*CO)*csw2;  dnhx=dcnx+(Nr+NH)*csw1; \
dchx=(Cr+CH)*csw3;           dccnx=Cr+CN; 
pnhx=dcnx+(Nr+NH+pnh)*csw1;  pchx=(Cr+CH+pch)*csw3;
dcoy=(Cr+Oa+2*CO)*ssw2;      dnhy=(Nr+NH)*ssw1;       dchy=(Cr+CH)*ssw3;  
pnhy=(Nr+NH+pnh)*ssw1;       pchy=(Cr+CH+pch)*ssw3;
dccox=(Cr+CO)*csw2;          dccoy=(Cr+CO)*ssw2;      
dolx=(Cr+2*CO+2*Oa+OL)*csw2; doly=(Cr+2*CO+2*Oa+OL)*ssw2;
dom1x=dcox-(Oa+OM)*ssw2;     dom2x=dcox+(Oa+OM)*ssw2;
dom1y=dcoy+(Oa+OM)*csw2;     dom2y=dcoy-(Oa+OM)*csw2;\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5668390630515456`*^9, 3.566839078760773*^9}, {
   3.566839318626795*^9, 3.5668393238528037`*^9}, {3.566839506201524*^9, 
   3.5668395108503323`*^9}, {3.566839541207986*^9, 3.566839568664034*^9}, {
   3.566839684650237*^9, 3.56683970887708*^9}, {3.5668397799352045`*^9, 
   3.5668398201988754`*^9}, {3.566839877872177*^9, 3.5668398909606*^9}, {
   3.5668404106911125`*^9, 3.566840417539525*^9}, {3.5668404815464373`*^9, 
   3.566840691460406*^9}, {3.566840732519678*^9, 3.5668408608923035`*^9}, {
   3.566840930000425*^9, 3.5668409412168446`*^9}, {3.566841374686006*^9, 
   3.566841375388007*^9}, {3.5668421835006266`*^9, 3.5668421845770283`*^9}, {
   3.566887398037596*^9, 3.566887402078003*^9}, {3.5668874322328568`*^9, 
   3.5668874455864797`*^9}, {3.566887493166563*^9, 3.5668874939309645`*^9}, {
   3.5668882488786907`*^9, 3.5668882511094947`*^9}, {3.56688847075788*^9, 
   3.566888479322295*^9}, {3.5668885903476906`*^9, 3.5668885954956994`*^9}, {
   3.566971858103361*^9, 3.566971918927868*^9}, {3.566987354128765*^9, 
   3.5669873624903803`*^9}, {3.566987488164201*^9, 3.566987535619484*^9}, {
   3.566987795843541*^9, 3.5669878015843515`*^9}, {3.566987856168847*^9, 
   3.56698786325126*^9}, {3.566988264733565*^9, 3.566988269819174*^9}, {
   3.566988366133743*^9, 3.5669883719525533`*^9}, {3.5669885431160536`*^9, 
   3.566988549309265*^9}, {3.566991946182942*^9, 3.5669919483981457`*^9}, 
   3.5669919866182127`*^9, {3.5669920226542764`*^9, 3.5669920505783253`*^9}, {
   3.566992103743219*^9, 3.5669921079396257`*^9}, 3.5669921962201815`*^9, {
   3.5669922303998413`*^9, 3.5669922340658474`*^9}, {3.566992291177548*^9, 
   3.5669922926751504`*^9}, {3.566992333890423*^9, 3.5669923348264246`*^9}, 
   3.566992409675356*^9, {3.566997472043352*^9, 3.566997587171554*^9}, {
   3.566997674297707*^9, 3.5669977903463106`*^9}, {3.566997868580448*^9, 
   3.5669978780652647`*^9}, {3.566998052270771*^9, 3.566998071786405*^9}, {
   3.5669981665409718`*^9, 3.5669981730305834`*^9}, {3.56699824522751*^9, 
   3.566998326394452*^9}, {3.5669983623837156`*^9, 3.5669984017425847`*^9}, {
   3.566998470788306*^9, 3.5669984830811276`*^9}, {3.566998534405218*^9, 
   3.5669985995821323`*^9}, {3.5669988153461113`*^9, 3.566998936995125*^9}, 
   3.566998967134378*^9, 3.5669990224832754`*^9, {3.5669990973790064`*^9, 
   3.56699912179305*^9}, {3.5669993047345705`*^9, 3.5669993241410046`*^9}, {
   3.566999416274767*^9, 3.566999426742385*^9}, {3.566999497301309*^9, 
   3.5669995032137194`*^9}, {3.566999549733001*^9, 3.566999551433404*^9}, {
   3.5670567651502*^9, 3.567056829110312*^9}, {3.5670569077344503`*^9, 
   3.5670570287126627`*^9}, {3.56705707303234*^9, 3.5670571842137356`*^9}, {
   3.5670572582982655`*^9, 3.567057292103525*^9}, {3.56705732891959*^9, 
   3.5670574906294737`*^9}, {3.5670578301328697`*^9, 3.56705795512029*^9}, {
   3.5670580197668033`*^9, 3.567058030437222*^9}, {3.5670580810281105`*^9, 
   3.5670581611654515`*^9}, {3.5670581932547083`*^9, 3.567058240616391*^9}, {
   3.5670582801312604`*^9, 3.5670582990540934`*^9}, {3.567058343560972*^9, 
   3.5670584018894744`*^9}, {3.56705845049916*^9, 3.5670585083128614`*^9}, {
   3.567058578731385*^9, 3.567058593645011*^9}, {3.5670586348134837`*^9, 
   3.5670586577143235`*^9}, {3.567058689007979*^9, 3.5670587629209085`*^9}, {
   3.567058832777831*^9, 3.5670588394702425`*^9}, {3.567058955628047*^9, 
   3.5670589616496572`*^9}, {3.5670592284569263`*^9, 3.567059354442747*^9}, {
   3.5670593898548098`*^9, 3.567059560456709*^9}, {3.5670596222796183`*^9, 
   3.5670596362884426`*^9}, {3.5670596726677065`*^9, 3.567059775347087*^9}, {
   3.567059838449198*^9, 3.5670598453756094`*^9}, {3.5670598764976645`*^9, 
   3.5670598805380716`*^9}, 3.567059911254525*^9, {3.5670599744502363`*^9, 
   3.567060007818695*^9}, {3.567060059127185*^9, 3.5670600818564253`*^9}, {
   3.5670601143824825`*^9, 3.567060314593234*^9}, {3.5670607181347427`*^9, 
   3.5670607311607656`*^9}, {3.567060785573661*^9, 3.567061232998047*^9}, {
   3.5670626282800975`*^9, 3.56706267537658*^9}, 3.567062795371991*^9, {
   3.567069851801585*^9, 3.567069879460434*^9}, {3.567070767429593*^9, 
   3.5670707823588195`*^9}, {3.5670805354535503`*^9, 
   3.5670806312689185`*^9}, {3.567083393113369*^9, 3.5670834586022844`*^9}, {
   3.567083501299559*^9, 3.5670835804696984`*^9}, {3.5670836156477604`*^9, 
   3.5670836215757704`*^9}, {3.567083656504232*^9, 3.5670836570034328`*^9}, {
   3.5670837376399746`*^9, 3.56708374668799*^9}}],

Cell["\<\
(* 6 nuclei: coord *)
xn={0,-dcnx,dcox,-pnhx,-pnhx,pchx};
yn={0,0,dcoy,pnhy,-pnhy,-pchy};
zn={0,0,0,0,0,0};\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5668390630515456`*^9, 3.566839078760773*^9}, {
   3.566839318626795*^9, 3.5668393238528037`*^9}, {3.566839506201524*^9, 
   3.5668395108503323`*^9}, {3.566839541207986*^9, 3.566839568664034*^9}, {
   3.566839684650237*^9, 3.56683970887708*^9}, {3.5668397799352045`*^9, 
   3.5668398201988754`*^9}, {3.566839877872177*^9, 3.5668398909606*^9}, {
   3.5668404106911125`*^9, 3.566840417539525*^9}, {3.5668404815464373`*^9, 
   3.566840691460406*^9}, {3.566840732519678*^9, 3.5668408608923035`*^9}, {
   3.566840930000425*^9, 3.5668409412168446`*^9}, {3.566841374686006*^9, 
   3.566841375388007*^9}, {3.5668421835006266`*^9, 3.5668421845770283`*^9}, {
   3.566887398037596*^9, 3.566887402078003*^9}, {3.5668874322328568`*^9, 
   3.5668874455864797`*^9}, {3.566887493166563*^9, 3.5668874939309645`*^9}, {
   3.5668882488786907`*^9, 3.5668882511094947`*^9}, {3.56688847075788*^9, 
   3.566888479322295*^9}, {3.5668885903476906`*^9, 3.5668885954956994`*^9}, {
   3.566971858103361*^9, 3.566971918927868*^9}, {3.566987354128765*^9, 
   3.5669873624903803`*^9}, {3.566987488164201*^9, 3.566987535619484*^9}, {
   3.566987795843541*^9, 3.5669878015843515`*^9}, {3.566987856168847*^9, 
   3.56698786325126*^9}, {3.566988264733565*^9, 3.566988269819174*^9}, {
   3.566988366133743*^9, 3.5669883719525533`*^9}, {3.5669885431160536`*^9, 
   3.566988549309265*^9}, {3.566991946182942*^9, 3.5669919483981457`*^9}, 
   3.5669919866182127`*^9, {3.5669920226542764`*^9, 3.5669920505783253`*^9}, {
   3.566992103743219*^9, 3.5669921079396257`*^9}, 3.5669921962201815`*^9, {
   3.5669922303998413`*^9, 3.5669922340658474`*^9}, {3.566992291177548*^9, 
   3.5669922926751504`*^9}, {3.566992333890423*^9, 3.5669923348264246`*^9}, 
   3.566992409675356*^9, {3.566997472043352*^9, 3.566997587171554*^9}, {
   3.566997674297707*^9, 3.5669977903463106`*^9}, {3.566997868580448*^9, 
   3.5669978780652647`*^9}, {3.566998052270771*^9, 3.566998071786405*^9}, {
   3.5669981665409718`*^9, 3.5669981730305834`*^9}, {3.56699824522751*^9, 
   3.566998326394452*^9}, {3.5669983623837156`*^9, 3.5669984017425847`*^9}, {
   3.566998470788306*^9, 3.5669984830811276`*^9}, {3.566998534405218*^9, 
   3.5669985995821323`*^9}, {3.5669988153461113`*^9, 3.566998936995125*^9}, 
   3.566998967134378*^9, 3.5669990224832754`*^9, {3.5669990973790064`*^9, 
   3.56699912179305*^9}, {3.5669993047345705`*^9, 3.5669993241410046`*^9}, {
   3.566999416274767*^9, 3.566999426742385*^9}, {3.566999497301309*^9, 
   3.5669995032137194`*^9}, {3.566999549733001*^9, 3.566999551433404*^9}, {
   3.5670567651502*^9, 3.567056829110312*^9}, {3.5670569077344503`*^9, 
   3.5670570287126627`*^9}, {3.56705707303234*^9, 3.5670571842137356`*^9}, {
   3.5670572582982655`*^9, 3.567057292103525*^9}, {3.56705732891959*^9, 
   3.5670574906294737`*^9}, {3.5670578301328697`*^9, 3.56705795512029*^9}, {
   3.5670580197668033`*^9, 3.567058030437222*^9}, {3.5670580810281105`*^9, 
   3.5670581611654515`*^9}, {3.5670581932547083`*^9, 3.567058240616391*^9}, {
   3.5670582801312604`*^9, 3.5670582990540934`*^9}, {3.567058343560972*^9, 
   3.5670584018894744`*^9}, {3.56705845049916*^9, 3.5670585083128614`*^9}, {
   3.567058578731385*^9, 3.567058593645011*^9}, {3.5670586348134837`*^9, 
   3.5670586577143235`*^9}, {3.567058689007979*^9, 3.5670587629209085`*^9}, {
   3.567058832777831*^9, 3.5670588394702425`*^9}, {3.567058955628047*^9, 
   3.5670589616496572`*^9}, {3.5670592284569263`*^9, 3.567059354442747*^9}, {
   3.5670593898548098`*^9, 3.567059560456709*^9}, {3.5670596222796183`*^9, 
   3.5670596362884426`*^9}, {3.5670596726677065`*^9, 3.567059775347087*^9}, {
   3.567059838449198*^9, 3.5670598453756094`*^9}, {3.5670598764976645`*^9, 
   3.5670598805380716`*^9}, 3.567059911254525*^9, {3.5670599744502363`*^9, 
   3.567060007818695*^9}, {3.567060059127185*^9, 3.5670600818564253`*^9}, {
   3.5670601143824825`*^9, 3.567060314593234*^9}, {3.5670607181347427`*^9, 
   3.5670607311607656`*^9}, {3.567060785573661*^9, 3.567061232998047*^9}, {
   3.5670626282800975`*^9, 3.56706267537658*^9}, 3.567062795371991*^9, {
   3.567069851801585*^9, 3.567069879460434*^9}, {3.567070767429593*^9, 
   3.5670707823588195`*^9}, {3.5670805354535503`*^9, 
   3.5670806312689185`*^9}, {3.567083393113369*^9, 3.5670834586022844`*^9}, {
   3.567083501299559*^9, 3.5670835804696984`*^9}, {3.5670836156477604`*^9, 
   3.5670836215757704`*^9}, {3.567083656504232*^9, 3.5670836570034328`*^9}, {
   3.5670837376399746`*^9, 3.567083757514409*^9}}],

Cell["\<\
(* 15 clouds: coord *)        
xc={0,-dcnx,dcox,-dnhx,-dnhx,dchx,-dccnx,dccox,dolx,-dcnx,-dcnx,dccox,dccox,\
dom1x,dom2x};
yc={0,0,dcoy,dnhy,-dnhy,-dchy,0,dccoy,doly,0,0,dccoy,dccoy,dom1y,dom2y};
zc={0,0,0,0,0,0,0,0,0,NM,-NM,OC,-OC,0,0};\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5668390630515456`*^9, 3.566839078760773*^9}, {
   3.566839318626795*^9, 3.5668393238528037`*^9}, {3.566839506201524*^9, 
   3.5668395108503323`*^9}, {3.566839541207986*^9, 3.566839568664034*^9}, {
   3.566839684650237*^9, 3.56683970887708*^9}, {3.5668397799352045`*^9, 
   3.5668398201988754`*^9}, {3.566839877872177*^9, 3.5668398909606*^9}, {
   3.5668404106911125`*^9, 3.566840417539525*^9}, {3.5668404815464373`*^9, 
   3.566840691460406*^9}, {3.566840732519678*^9, 3.5668408608923035`*^9}, {
   3.566840930000425*^9, 3.5668409412168446`*^9}, {3.566841374686006*^9, 
   3.566841375388007*^9}, {3.5668421835006266`*^9, 3.5668421845770283`*^9}, {
   3.566887398037596*^9, 3.566887402078003*^9}, {3.5668874322328568`*^9, 
   3.5668874455864797`*^9}, {3.566887493166563*^9, 3.5668874939309645`*^9}, {
   3.5668882488786907`*^9, 3.5668882511094947`*^9}, {3.56688847075788*^9, 
   3.566888479322295*^9}, {3.5668885903476906`*^9, 3.5668885954956994`*^9}, {
   3.566971858103361*^9, 3.566971918927868*^9}, {3.566987354128765*^9, 
   3.5669873624903803`*^9}, {3.566987488164201*^9, 3.566987535619484*^9}, {
   3.566987795843541*^9, 3.5669878015843515`*^9}, {3.566987856168847*^9, 
   3.56698786325126*^9}, {3.566988264733565*^9, 3.566988269819174*^9}, {
   3.566988366133743*^9, 3.5669883719525533`*^9}, {3.5669885431160536`*^9, 
   3.566988549309265*^9}, {3.566991946182942*^9, 3.5669919483981457`*^9}, 
   3.5669919866182127`*^9, {3.5669920226542764`*^9, 3.5669920505783253`*^9}, {
   3.566992103743219*^9, 3.5669921079396257`*^9}, 3.5669921962201815`*^9, {
   3.5669922303998413`*^9, 3.5669922340658474`*^9}, {3.566992291177548*^9, 
   3.5669922926751504`*^9}, {3.566992333890423*^9, 3.5669923348264246`*^9}, 
   3.566992409675356*^9, {3.566997472043352*^9, 3.566997587171554*^9}, {
   3.566997674297707*^9, 3.5669977903463106`*^9}, {3.566997868580448*^9, 
   3.5669978780652647`*^9}, {3.566998052270771*^9, 3.566998071786405*^9}, {
   3.5669981665409718`*^9, 3.5669981730305834`*^9}, {3.56699824522751*^9, 
   3.566998326394452*^9}, {3.5669983623837156`*^9, 3.5669984017425847`*^9}, {
   3.566998470788306*^9, 3.5669984830811276`*^9}, {3.566998534405218*^9, 
   3.5669985995821323`*^9}, {3.5669988153461113`*^9, 3.566998936995125*^9}, 
   3.566998967134378*^9, 3.5669990224832754`*^9, {3.5669990973790064`*^9, 
   3.56699912179305*^9}, {3.5669993047345705`*^9, 3.5669993241410046`*^9}, {
   3.566999416274767*^9, 3.566999426742385*^9}, {3.566999497301309*^9, 
   3.5669995032137194`*^9}, {3.566999549733001*^9, 3.566999551433404*^9}, {
   3.5670567651502*^9, 3.567056829110312*^9}, {3.5670569077344503`*^9, 
   3.5670570287126627`*^9}, {3.56705707303234*^9, 3.5670571842137356`*^9}, {
   3.5670572582982655`*^9, 3.567057292103525*^9}, {3.56705732891959*^9, 
   3.5670574906294737`*^9}, {3.5670578301328697`*^9, 3.56705795512029*^9}, {
   3.5670580197668033`*^9, 3.567058030437222*^9}, {3.5670580810281105`*^9, 
   3.5670581611654515`*^9}, {3.5670581932547083`*^9, 3.567058240616391*^9}, {
   3.5670582801312604`*^9, 3.5670582990540934`*^9}, {3.567058343560972*^9, 
   3.5670584018894744`*^9}, {3.56705845049916*^9, 3.5670585083128614`*^9}, {
   3.567058578731385*^9, 3.567058593645011*^9}, {3.5670586348134837`*^9, 
   3.5670586577143235`*^9}, {3.567058689007979*^9, 3.5670587629209085`*^9}, {
   3.567058832777831*^9, 3.5670588394702425`*^9}, {3.567058955628047*^9, 
   3.5670589616496572`*^9}, {3.5670592284569263`*^9, 3.567059354442747*^9}, {
   3.5670593898548098`*^9, 3.567059560456709*^9}, {3.5670596222796183`*^9, 
   3.5670596362884426`*^9}, {3.5670596726677065`*^9, 3.567059775347087*^9}, {
   3.567059838449198*^9, 3.5670598453756094`*^9}, {3.5670598764976645`*^9, 
   3.5670598805380716`*^9}, 3.567059911254525*^9, {3.5670599744502363`*^9, 
   3.567060007818695*^9}, {3.567060059127185*^9, 3.5670600818564253`*^9}, {
   3.5670601143824825`*^9, 3.567060314593234*^9}, {3.5670607181347427`*^9, 
   3.5670607311607656`*^9}, {3.567060785573661*^9, 3.567061232998047*^9}, {
   3.5670626282800975`*^9, 3.56706267537658*^9}, 3.567062795371991*^9, {
   3.567069851801585*^9, 3.567069879460434*^9}, {3.567070767429593*^9, 
   3.5670707823588195`*^9}, {3.5670805354535503`*^9, 
   3.5670806312689185`*^9}, {3.567083393113369*^9, 3.5670834586022844`*^9}, {
   3.567083501299559*^9, 3.5670835804696984`*^9}, {3.5670836156477604`*^9, 
   3.5670836215757704`*^9}, {3.567083656504232*^9, 3.5670836570034328`*^9}, {
   3.5670837376399746`*^9, 3.5670837651584225`*^9}}],

Cell["\<\
(* epot terms *)
vne=vne-2*(3-(pnh/NH)^2)/NH;
vne=vne-(3-(pch/CH)^2)/CH;
vnn = 0.0;\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5668390630515456`*^9, 3.566839078760773*^9}, {
   3.566839318626795*^9, 3.5668393238528037`*^9}, {3.566839506201524*^9, 
   3.5668395108503323`*^9}, {3.566839541207986*^9, 3.566839568664034*^9}, {
   3.566839684650237*^9, 3.56683970887708*^9}, {3.5668397799352045`*^9, 
   3.5668398201988754`*^9}, {3.566839877872177*^9, 3.5668398909606*^9}, {
   3.5668404106911125`*^9, 3.566840417539525*^9}, {3.5668404815464373`*^9, 
   3.566840691460406*^9}, {3.566840732519678*^9, 3.5668408608923035`*^9}, {
   3.566840930000425*^9, 3.5668409412168446`*^9}, {3.566841374686006*^9, 
   3.566841375388007*^9}, {3.5668421835006266`*^9, 3.5668421845770283`*^9}, {
   3.566887398037596*^9, 3.566887402078003*^9}, {3.5668874322328568`*^9, 
   3.5668874455864797`*^9}, {3.566887493166563*^9, 3.5668874939309645`*^9}, {
   3.5668882488786907`*^9, 3.5668882511094947`*^9}, {3.56688847075788*^9, 
   3.566888479322295*^9}, {3.5668885903476906`*^9, 3.5668885954956994`*^9}, {
   3.566971858103361*^9, 3.566971918927868*^9}, {3.566987354128765*^9, 
   3.5669873624903803`*^9}, {3.566987488164201*^9, 3.566987535619484*^9}, {
   3.566987795843541*^9, 3.5669878015843515`*^9}, {3.566987856168847*^9, 
   3.56698786325126*^9}, {3.566988264733565*^9, 3.566988269819174*^9}, {
   3.566988366133743*^9, 3.5669883719525533`*^9}, {3.5669885431160536`*^9, 
   3.566988549309265*^9}, {3.566991946182942*^9, 3.5669919483981457`*^9}, 
   3.5669919866182127`*^9, {3.5669920226542764`*^9, 3.5669920505783253`*^9}, {
   3.566992103743219*^9, 3.5669921079396257`*^9}, 3.5669921962201815`*^9, {
   3.5669922303998413`*^9, 3.5669922340658474`*^9}, {3.566992291177548*^9, 
   3.5669922926751504`*^9}, {3.566992333890423*^9, 3.5669923348264246`*^9}, 
   3.566992409675356*^9, {3.566997472043352*^9, 3.566997587171554*^9}, {
   3.566997674297707*^9, 3.5669977903463106`*^9}, {3.566997868580448*^9, 
   3.5669978780652647`*^9}, {3.566998052270771*^9, 3.566998071786405*^9}, {
   3.5669981665409718`*^9, 3.5669981730305834`*^9}, {3.56699824522751*^9, 
   3.566998326394452*^9}, {3.5669983623837156`*^9, 3.5669984017425847`*^9}, {
   3.566998470788306*^9, 3.5669984830811276`*^9}, {3.566998534405218*^9, 
   3.5669985995821323`*^9}, {3.5669988153461113`*^9, 3.566998936995125*^9}, 
   3.566998967134378*^9, 3.5669990224832754`*^9, {3.5669990973790064`*^9, 
   3.56699912179305*^9}, {3.5669993047345705`*^9, 3.5669993241410046`*^9}, {
   3.566999416274767*^9, 3.566999426742385*^9}, {3.566999497301309*^9, 
   3.5669995032137194`*^9}, {3.566999549733001*^9, 3.566999551433404*^9}, {
   3.5670567651502*^9, 3.567056829110312*^9}, {3.5670569077344503`*^9, 
   3.5670570287126627`*^9}, {3.56705707303234*^9, 3.5670571842137356`*^9}, {
   3.5670572582982655`*^9, 3.567057292103525*^9}, {3.56705732891959*^9, 
   3.5670574906294737`*^9}, {3.5670578301328697`*^9, 3.56705795512029*^9}, {
   3.5670580197668033`*^9, 3.567058030437222*^9}, {3.5670580810281105`*^9, 
   3.5670581611654515`*^9}, {3.5670581932547083`*^9, 3.567058240616391*^9}, {
   3.5670582801312604`*^9, 3.5670582990540934`*^9}, {3.567058343560972*^9, 
   3.5670584018894744`*^9}, {3.56705845049916*^9, 3.5670585083128614`*^9}, {
   3.567058578731385*^9, 3.567058593645011*^9}, {3.5670586348134837`*^9, 
   3.5670586577143235`*^9}, {3.567058689007979*^9, 3.5670587629209085`*^9}, {
   3.567058832777831*^9, 3.5670588394702425`*^9}, {3.567058955628047*^9, 
   3.5670589616496572`*^9}, {3.5670592284569263`*^9, 3.567059354442747*^9}, {
   3.5670593898548098`*^9, 3.567059560456709*^9}, {3.5670596222796183`*^9, 
   3.5670596362884426`*^9}, {3.5670596726677065`*^9, 3.567059775347087*^9}, {
   3.567059838449198*^9, 3.5670598453756094`*^9}, {3.5670598764976645`*^9, 
   3.5670598805380716`*^9}, 3.567059911254525*^9, {3.5670599744502363`*^9, 
   3.567060007818695*^9}, {3.567060059127185*^9, 3.5670600818564253`*^9}, {
   3.5670601143824825`*^9, 3.567060314593234*^9}, {3.5670607181347427`*^9, 
   3.5670607311607656`*^9}, {3.567060785573661*^9, 3.567061232998047*^9}, {
   3.5670626282800975`*^9, 3.56706267537658*^9}, 3.567062795371991*^9, {
   3.567069851801585*^9, 3.567069879460434*^9}, {3.567070767429593*^9, 
   3.5670707823588195`*^9}, {3.5670805354535503`*^9, 
   3.5670806312689185`*^9}, {3.567083393113369*^9, 3.5670834586022844`*^9}, {
   3.567083501299559*^9, 3.5670835804696984`*^9}, {3.5670836156477604`*^9, 
   3.5670836215757704`*^9}, {3.567083656504232*^9, 3.5670836570034328`*^9}, {
   3.5670837376399746`*^9, 3.5670837651584225`*^9}}],

Cell["\<\
(*ww*)
For[i = 1, i < nc, i++,
  For[j = i+1, j < nc+1, j++,
    vee = vee + \
oc[[i]]*oc[[j]]/Sqrt[(xc[[i]]-xc[[j]])^2+(yc[[i]]-yc[[j]])^2+(zc[[i]]-zc[[j]])\
^2]]]\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5668390630515456`*^9, 3.566839083924382*^9}, {
  3.5668410500738354`*^9, 3.5668411059531336`*^9}, {3.5668412393245683`*^9, 
  3.5668412393713684`*^9}, {3.5668412896190567`*^9, 3.5668413118646955`*^9}, {
  3.566841703815384*^9, 3.5668417040025845`*^9}}],

Cell["\<\
(*nn*)
For[i = 1, i < nn, i++,
  For[j = i+1, j < nn+1, j++,
    vnn = vnn + \
ch[[i]]*ch[[j]]/Sqrt[(xn[[i]]-xn[[j]])^2+(yn[[i]]-yn[[j]])^2+(zn[[i]]-zn[[j]])\
^2]]]\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5668390630515456`*^9, 3.5668390890567913`*^9}, {
  3.566841524165468*^9, 3.566841530218279*^9}, {3.5670590745470557`*^9, 
  3.567059089304682*^9}}],

Cell["\<\
(*nw*)
For[i = 1, i < nc+1, i++,
  For[j = 1, j < nn+1, j++,
    If[i \[NotEqual] j,
      vne = vne + \
oc[[i]]*ch[[j]]/Sqrt[(xc[[i]]-xn[[j]])^2+(yc[[i]]-yn[[j]])^2+(zc[[i]]-zn[[j]])\
^2]]]]\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5668390630515456`*^9, 3.5668390979332066`*^9}, 
   3.5668412394337683`*^9, {3.566841289743857*^9, 3.566841311927096*^9}, 
   3.5668415580955276`*^9, {3.5668416231788425`*^9, 3.5668416346760626`*^9}, {
   3.5669003862430363`*^9, 3.5669003862430363`*^9}, {3.5670591033447065`*^9, 
   3.5670591078063145`*^9}}],

Cell[CellGroupData[{

Cell["\<\
Epot=vne+vee+vnn;
func=Ekin+Epot;
t = FindMinimum[func,{Cr,0.26},{Nr,0.22},{Oa,0.186},{NH,1.11},{CH,1.12},{CO,0.\
92},{CN,1.03},{OL,0.905},{NM,0.765},{OC,1.36},{OM,0.86},{w1,2.15},{w2,1.0},{\
w3,1.12},{pch,0.51},{pnh,0.5},{Method -> \"Newton\"}, {MaxIterations -> 800}]\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5668390630515456`*^9, 3.566839104968819*^9}, {
   3.5668402974973135`*^9, 3.5668403782742558`*^9}, {3.5668414049812593`*^9, 
   3.5668414375229163`*^9}, 3.566841756605877*^9, {3.5668429499163733`*^9, 
   3.566842978277223*^9}, {3.5668880399631233`*^9, 3.5668880687763743`*^9}, {
   3.56697197120356*^9, 3.5669719847755833`*^9}, {3.5669875623735313`*^9, 
   3.566987566460738*^9}, 3.566987810538767*^9, {3.566987882142893*^9, 
   3.5669878871349015`*^9}, {3.566987967677843*^9, 3.5669880000323*^9}, {
   3.5669882132846746`*^9, 3.5669882156558785`*^9}, 3.5669882806455927`*^9, {
   3.5669883912029867`*^9, 3.5669884039482093`*^9}, 3.5669885600420837`*^9, {
   3.5670614496356277`*^9, 3.5670616095515084`*^9}, {3.56706164472957*^9, 
   3.567061671780018*^9}, {3.5670617748181987`*^9, 3.567061891771604*^9}, {
   3.5670620082414083`*^9, 3.5670620167278233`*^9}, {3.5670640067297187`*^9, 
   3.567064010270925*^9}, {3.5670640525313993`*^9, 3.567064055854205*^9}, 
   3.567064094339473*^9, {3.5670689484819984`*^9, 3.5670689607748203`*^9}, 
   3.5670698952320614`*^9, {3.56708497059774*^9, 3.567084977087352*^9}, {
   3.567085213599367*^9, 3.5670852679342623`*^9}, {3.567085494899061*^9, 
   3.567085499485469*^9}, {3.5670859700758953`*^9, 3.5670859708870974`*^9}, {
   3.567086006080759*^9, 3.5670861172933545`*^9}, {3.5670866938391666`*^9, 
   3.567086753244071*^9}, {3.5670871411699524`*^9, 3.5670871501399684`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "168.06916898729622`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"Cr", "\[Rule]", "0.2597231856403463`"}], ",", 
     RowBox[{"Nr", "\[Rule]", "0.22405135815353058`"}], ",", 
     RowBox[{"Oa", "\[Rule]", "0.1916417225081855`"}], ",", 
     RowBox[{"NH", "\[Rule]", "1.1619297579240278`"}], ",", 
     RowBox[{"CH", "\[Rule]", "1.2532410401053755`"}], ",", 
     RowBox[{"CO", "\[Rule]", "0.9166962134031293`"}], ",", 
     RowBox[{"CN", "\[Rule]", "1.0305980699531003`"}], ",", 
     RowBox[{"OL", "\[Rule]", "0.8932960411034288`"}], ",", 
     RowBox[{"NM", "\[Rule]", "0.7472908218886891`"}], ",", 
     RowBox[{"OC", "\[Rule]", "1.5679221435272637`"}], ",", 
     RowBox[{"OM", "\[Rule]", "0.8346241451827114`"}], ",", 
     RowBox[{"w1", "\[Rule]", "2.1564408390806467`"}], ",", 
     RowBox[{"w2", "\[Rule]", "0.9700759023513378`"}], ",", 
     RowBox[{"w3", "\[Rule]", "1.1571282820932742`"}], ",", 
     RowBox[{"pch", "\[Rule]", "0.572080058406644`"}], ",", 
     RowBox[{"pnh", "\[Rule]", "0.5183973604881592`"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.567085316559548*^9, 3.5670855076442833`*^9, 3.5670858936357613`*^9, 
   3.5670859776419086`*^9, 3.5670861256081686`*^9, 3.5670863498429627`*^9, 
   3.5670863799354153`*^9, 3.5670864253782954`*^9, {3.5670864623503604`*^9, 
   3.567086520866063*^9}, {3.5670865519881177`*^9, 3.5670865924389887`*^9}, 
   3.5670867631032887`*^9, 3.567087074027435*^9, 3.567087206019266*^9, 
   3.567087321334669*^9, 3.5670873552335286`*^9, 3.5670874041084146`*^9, {
   3.5670874360884705`*^9, 3.5670874603309135`*^9}, 3.5670874999393826`*^9, {
   3.5670979102472944`*^9, 3.5670979311981316`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\

Ekin /. t[[2]]
vne /. t[[2]]
vee /. t[[2]]
vnn /. t[[2]]
-Epot/Ekin /. t[[2]]
w1*rad /. t[[2]]
w2*rad /. t[[2]]
w3*rad /. t[[2]]
(Cr+Oa+2*CO)*bohr /. t[[2]]
(Cr+Nr+2*CN)*bohr /. t[[2]]
(Cr+CH+pch)*bohr /. t[[2]]
(Nr+NH+pnh)*bohr /. t[[2]]\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5668390630515456`*^9, 3.566839112846833*^9}, {
  3.566842810779729*^9, 3.566842811091729*^9}, {3.5669879418285975`*^9, 
  3.5669879480686083`*^9}, {3.5670621290636206`*^9, 3.567062152463662*^9}, {
  3.567063905688341*^9, 3.5670639108675504`*^9}, {3.5670689307603674`*^9, 
  3.5670689312127686`*^9}, {3.5670696234639845`*^9, 3.5670696694216647`*^9}, {
  3.5670795267235785`*^9, 3.5670795309823856`*^9}, {3.567079591604092*^9, 
  3.567079592524494*^9}, {3.567079628810158*^9, 3.567079717433913*^9}, {
  3.567082843118803*^9, 3.567082848547613*^9}, {3.567082895987296*^9, 
  3.567082897079298*^9}, {3.567085039347061*^9, 3.5670851279708166`*^9}}],

Cell[BoxData["168.06916898729625`"], "Output",
 CellChangeTimes->{
  3.5670784509066887`*^9, 3.5670785564876738`*^9, 3.5670793607324867`*^9, 
   3.567079418989789*^9, 3.567079722519522*^9, {3.5670803370212016`*^9, 
   3.5670804052089214`*^9}, 3.567080471727438*^9, 3.567081090393325*^9, 
   3.567081122451381*^9, 3.5670811538698363`*^9, 3.5670812010599194`*^9, 
   3.5670812378915834`*^9, 3.567081274380048*^9, 3.5670824912601852`*^9, 
   3.5670828555988255`*^9, 3.567082903334909*^9, 3.5670836938662977`*^9, {
   3.5670844599392433`*^9, 3.567084518517346*^9}, 3.567084593257077*^9, 
   3.5670846788388276`*^9, 3.567085316575148*^9, 3.5670855076598835`*^9, 
   3.5670858936669617`*^9, 3.567085977657509*^9, 3.567086125623769*^9, 
   3.5670863498585625`*^9, 3.5670863799510155`*^9, 3.567086425393895*^9, {
   3.5670864623815603`*^9, 3.5670865208816633`*^9}, {3.567086552003718*^9, 
   3.567086592470189*^9}, 3.567086763118889*^9, 3.5670870740430346`*^9, 
   3.5670872060348663`*^9, 3.5670873213502693`*^9, 3.5670873552491283`*^9, 
   3.567087404139614*^9, {3.5670874361040707`*^9, 3.5670874603465133`*^9}, 
   3.5670874999549828`*^9, {3.5670979102628946`*^9, 3.567097931213732*^9}}],

Cell[BoxData[
 RowBox[{"-", "538.3611060012901`"}]], "Output",
 CellChangeTimes->{
  3.5670784509066887`*^9, 3.5670785564876738`*^9, 3.5670793607324867`*^9, 
   3.567079418989789*^9, 3.567079722519522*^9, {3.5670803370212016`*^9, 
   3.5670804052089214`*^9}, 3.567080471727438*^9, 3.567081090393325*^9, 
   3.567081122451381*^9, 3.5670811538698363`*^9, 3.5670812010599194`*^9, 
   3.5670812378915834`*^9, 3.567081274380048*^9, 3.5670824912601852`*^9, 
   3.5670828555988255`*^9, 3.567082903334909*^9, 3.5670836938662977`*^9, {
   3.5670844599392433`*^9, 3.567084518517346*^9}, 3.567084593257077*^9, 
   3.5670846788388276`*^9, 3.567085316575148*^9, 3.5670855076598835`*^9, 
   3.5670858936669617`*^9, 3.567085977657509*^9, 3.567086125623769*^9, 
   3.5670863498585625`*^9, 3.5670863799510155`*^9, 3.567086425393895*^9, {
   3.5670864623815603`*^9, 3.5670865208816633`*^9}, {3.567086552003718*^9, 
   3.567086592470189*^9}, 3.567086763118889*^9, 3.5670870740430346`*^9, 
   3.5670872060348663`*^9, 3.5670873213502693`*^9, 3.5670873552491283`*^9, 
   3.567087404139614*^9, {3.5670874361040707`*^9, 3.5670874603465133`*^9}, 
   3.5670874999549828`*^9, {3.5670979102628946`*^9, 3.567097931213732*^9}}],

Cell[BoxData["130.60649695673055`"], "Output",
 CellChangeTimes->{
  3.5670784509066887`*^9, 3.5670785564876738`*^9, 3.5670793607324867`*^9, 
   3.567079418989789*^9, 3.567079722519522*^9, {3.5670803370212016`*^9, 
   3.5670804052089214`*^9}, 3.567080471727438*^9, 3.567081090393325*^9, 
   3.567081122451381*^9, 3.5670811538698363`*^9, 3.5670812010599194`*^9, 
   3.5670812378915834`*^9, 3.567081274380048*^9, 3.5670824912601852`*^9, 
   3.5670828555988255`*^9, 3.567082903334909*^9, 3.5670836938662977`*^9, {
   3.5670844599392433`*^9, 3.567084518517346*^9}, 3.567084593257077*^9, 
   3.5670846788388276`*^9, 3.567085316575148*^9, 3.5670855076598835`*^9, 
   3.5670858936669617`*^9, 3.567085977657509*^9, 3.567086125623769*^9, 
   3.5670863498585625`*^9, 3.5670863799510155`*^9, 3.567086425393895*^9, {
   3.5670864623815603`*^9, 3.5670865208816633`*^9}, {3.567086552003718*^9, 
   3.567086592470189*^9}, 3.567086763118889*^9, 3.5670870740430346`*^9, 
   3.5670872060348663`*^9, 3.5670873213502693`*^9, 3.5670873552491283`*^9, 
   3.567087404139614*^9, {3.5670874361040707`*^9, 3.5670874603465133`*^9}, 
   3.5670874999549828`*^9, {3.5670979102628946`*^9, 3.567097931213732*^9}}],

Cell[BoxData["71.6162710699671`"], "Output",
 CellChangeTimes->{
  3.5670784509066887`*^9, 3.5670785564876738`*^9, 3.5670793607324867`*^9, 
   3.567079418989789*^9, 3.567079722519522*^9, {3.5670803370212016`*^9, 
   3.5670804052089214`*^9}, 3.567080471727438*^9, 3.567081090393325*^9, 
   3.567081122451381*^9, 3.5670811538698363`*^9, 3.5670812010599194`*^9, 
   3.5670812378915834`*^9, 3.567081274380048*^9, 3.5670824912601852`*^9, 
   3.5670828555988255`*^9, 3.567082903334909*^9, 3.5670836938662977`*^9, {
   3.5670844599392433`*^9, 3.567084518517346*^9}, 3.567084593257077*^9, 
   3.5670846788388276`*^9, 3.567085316575148*^9, 3.5670855076598835`*^9, 
   3.5670858936669617`*^9, 3.567085977657509*^9, 3.567086125623769*^9, 
   3.5670863498585625`*^9, 3.5670863799510155`*^9, 3.567086425393895*^9, {
   3.5670864623815603`*^9, 3.5670865208816633`*^9}, {3.567086552003718*^9, 
   3.567086592470189*^9}, 3.567086763118889*^9, 3.5670870740430346`*^9, 
   3.5670872060348663`*^9, 3.5670873213502693`*^9, 3.5670873552491283`*^9, 
   3.567087404139614*^9, {3.5670874361040707`*^9, 3.5670874603465133`*^9}, 
   3.5670874999549828`*^9, {3.5670979102628946`*^9, 3.567097931213732*^9}}],

Cell[BoxData["1.9999999999999996`"], "Output",
 CellChangeTimes->{
  3.5670784509066887`*^9, 3.5670785564876738`*^9, 3.5670793607324867`*^9, 
   3.567079418989789*^9, 3.567079722519522*^9, {3.5670803370212016`*^9, 
   3.5670804052089214`*^9}, 3.567080471727438*^9, 3.567081090393325*^9, 
   3.567081122451381*^9, 3.5670811538698363`*^9, 3.5670812010599194`*^9, 
   3.5670812378915834`*^9, 3.567081274380048*^9, 3.5670824912601852`*^9, 
   3.5670828555988255`*^9, 3.567082903334909*^9, 3.5670836938662977`*^9, {
   3.5670844599392433`*^9, 3.567084518517346*^9}, 3.567084593257077*^9, 
   3.5670846788388276`*^9, 3.567085316575148*^9, 3.5670855076598835`*^9, 
   3.5670858936669617`*^9, 3.567085977657509*^9, 3.567086125623769*^9, 
   3.5670863498585625`*^9, 3.5670863799510155`*^9, 3.567086425393895*^9, {
   3.5670864623815603`*^9, 3.5670865208816633`*^9}, {3.567086552003718*^9, 
   3.567086592470189*^9}, 3.567086763118889*^9, 3.5670870740430346`*^9, 
   3.5670872060348663`*^9, 3.5670873213502693`*^9, 3.5670873552491283`*^9, 
   3.567087404139614*^9, {3.5670874361040707`*^9, 3.5670874603465133`*^9}, 
   3.5670874999549828`*^9, {3.5670979102628946`*^9, 3.5670979312293315`*^9}}],

Cell[BoxData["123.55495989898013`"], "Output",
 CellChangeTimes->{
  3.5670784509066887`*^9, 3.5670785564876738`*^9, 3.5670793607324867`*^9, 
   3.567079418989789*^9, 3.567079722519522*^9, {3.5670803370212016`*^9, 
   3.5670804052089214`*^9}, 3.567080471727438*^9, 3.567081090393325*^9, 
   3.567081122451381*^9, 3.5670811538698363`*^9, 3.5670812010599194`*^9, 
   3.5670812378915834`*^9, 3.567081274380048*^9, 3.5670824912601852`*^9, 
   3.5670828555988255`*^9, 3.567082903334909*^9, 3.5670836938662977`*^9, {
   3.5670844599392433`*^9, 3.567084518517346*^9}, 3.567084593257077*^9, 
   3.5670846788388276`*^9, 3.567085316575148*^9, 3.5670855076598835`*^9, 
   3.5670858936669617`*^9, 3.567085977657509*^9, 3.567086125623769*^9, 
   3.5670863498585625`*^9, 3.5670863799510155`*^9, 3.567086425393895*^9, {
   3.5670864623815603`*^9, 3.5670865208816633`*^9}, {3.567086552003718*^9, 
   3.567086592470189*^9}, 3.567086763118889*^9, 3.5670870740430346`*^9, 
   3.5670872060348663`*^9, 3.5670873213502693`*^9, 3.5670873552491283`*^9, 
   3.567087404139614*^9, {3.5670874361040707`*^9, 3.5670874603465133`*^9}, 
   3.5670874999549828`*^9, {3.5670979102628946`*^9, 3.5670979312293315`*^9}}],

Cell[BoxData["55.581255484423735`"], "Output",
 CellChangeTimes->{
  3.5670784509066887`*^9, 3.5670785564876738`*^9, 3.5670793607324867`*^9, 
   3.567079418989789*^9, 3.567079722519522*^9, {3.5670803370212016`*^9, 
   3.5670804052089214`*^9}, 3.567080471727438*^9, 3.567081090393325*^9, 
   3.567081122451381*^9, 3.5670811538698363`*^9, 3.5670812010599194`*^9, 
   3.5670812378915834`*^9, 3.567081274380048*^9, 3.5670824912601852`*^9, 
   3.5670828555988255`*^9, 3.567082903334909*^9, 3.5670836938662977`*^9, {
   3.5670844599392433`*^9, 3.567084518517346*^9}, 3.567084593257077*^9, 
   3.5670846788388276`*^9, 3.567085316575148*^9, 3.5670855076598835`*^9, 
   3.5670858936669617`*^9, 3.567085977657509*^9, 3.567086125623769*^9, 
   3.5670863498585625`*^9, 3.5670863799510155`*^9, 3.567086425393895*^9, {
   3.5670864623815603`*^9, 3.5670865208816633`*^9}, {3.567086552003718*^9, 
   3.567086592470189*^9}, 3.567086763118889*^9, 3.5670870740430346`*^9, 
   3.5670872060348663`*^9, 3.5670873213502693`*^9, 3.5670873552491283`*^9, 
   3.567087404139614*^9, {3.5670874361040707`*^9, 3.5670874603465133`*^9}, 
   3.5670874999549828`*^9, {3.5670979102628946`*^9, 3.5670979312293315`*^9}}],

Cell[BoxData["66.29856748259418`"], "Output",
 CellChangeTimes->{
  3.5670784509066887`*^9, 3.5670785564876738`*^9, 3.5670793607324867`*^9, 
   3.567079418989789*^9, 3.567079722519522*^9, {3.5670803370212016`*^9, 
   3.5670804052089214`*^9}, 3.567080471727438*^9, 3.567081090393325*^9, 
   3.567081122451381*^9, 3.5670811538698363`*^9, 3.5670812010599194`*^9, 
   3.5670812378915834`*^9, 3.567081274380048*^9, 3.5670824912601852`*^9, 
   3.5670828555988255`*^9, 3.567082903334909*^9, 3.5670836938662977`*^9, {
   3.5670844599392433`*^9, 3.567084518517346*^9}, 3.567084593257077*^9, 
   3.5670846788388276`*^9, 3.567085316575148*^9, 3.5670855076598835`*^9, 
   3.5670858936669617`*^9, 3.567085977657509*^9, 3.567086125623769*^9, 
   3.5670863498585625`*^9, 3.5670863799510155`*^9, 3.567086425393895*^9, {
   3.5670864623815603`*^9, 3.5670865208816633`*^9}, {3.567086552003718*^9, 
   3.567086592470189*^9}, 3.567086763118889*^9, 3.5670870740430346`*^9, 
   3.5670872060348663`*^9, 3.5670873213502693`*^9, 3.5670873552491283`*^9, 
   3.567087404139614*^9, {3.5670874361040707`*^9, 3.5670874603465133`*^9}, 
   3.5670874999549828`*^9, {3.5670979102628946`*^9, 3.5670979312293315`*^9}}],

Cell[BoxData["1.2090410322393712`"], "Output",
 CellChangeTimes->{
  3.5670784509066887`*^9, 3.5670785564876738`*^9, 3.5670793607324867`*^9, 
   3.567079418989789*^9, 3.567079722519522*^9, {3.5670803370212016`*^9, 
   3.5670804052089214`*^9}, 3.567080471727438*^9, 3.567081090393325*^9, 
   3.567081122451381*^9, 3.5670811538698363`*^9, 3.5670812010599194`*^9, 
   3.5670812378915834`*^9, 3.567081274380048*^9, 3.5670824912601852`*^9, 
   3.5670828555988255`*^9, 3.567082903334909*^9, 3.5670836938662977`*^9, {
   3.5670844599392433`*^9, 3.567084518517346*^9}, 3.567084593257077*^9, 
   3.5670846788388276`*^9, 3.567085316575148*^9, 3.5670855076598835`*^9, 
   3.5670858936669617`*^9, 3.567085977657509*^9, 3.567086125623769*^9, 
   3.5670863498585625`*^9, 3.5670863799510155`*^9, 3.567086425393895*^9, {
   3.5670864623815603`*^9, 3.5670865208816633`*^9}, {3.567086552003718*^9, 
   3.567086592470189*^9}, 3.567086763118889*^9, 3.5670870740430346`*^9, 
   3.5670872060348663`*^9, 3.5670873213502693`*^9, 3.5670873552491283`*^9, 
   3.567087404139614*^9, {3.5670874361040707`*^9, 3.5670874603465133`*^9}, 
   3.5670874999549828`*^9, {3.5670979102628946`*^9, 3.5670979312293315`*^9}}],

Cell[BoxData["1.346739951488356`"], "Output",
 CellChangeTimes->{
  3.5670784509066887`*^9, 3.5670785564876738`*^9, 3.5670793607324867`*^9, 
   3.567079418989789*^9, 3.567079722519522*^9, {3.5670803370212016`*^9, 
   3.5670804052089214`*^9}, 3.567080471727438*^9, 3.567081090393325*^9, 
   3.567081122451381*^9, 3.5670811538698363`*^9, 3.5670812010599194`*^9, 
   3.5670812378915834`*^9, 3.567081274380048*^9, 3.5670824912601852`*^9, 
   3.5670828555988255`*^9, 3.567082903334909*^9, 3.5670836938662977`*^9, {
   3.5670844599392433`*^9, 3.567084518517346*^9}, 3.567084593257077*^9, 
   3.5670846788388276`*^9, 3.567085316575148*^9, 3.5670855076598835`*^9, 
   3.5670858936669617`*^9, 3.567085977657509*^9, 3.567086125623769*^9, 
   3.5670863498585625`*^9, 3.5670863799510155`*^9, 3.567086425393895*^9, {
   3.5670864623815603`*^9, 3.5670865208816633`*^9}, {3.567086552003718*^9, 
   3.567086592470189*^9}, 3.567086763118889*^9, 3.5670870740430346`*^9, 
   3.5670872060348663`*^9, 3.5670873213502693`*^9, 3.5670873552491283`*^9, 
   3.567087404139614*^9, {3.5670874361040707`*^9, 3.5670874603465133`*^9}, 
   3.5670874999549828`*^9, {3.5670979102628946`*^9, 3.5670979312293315`*^9}}],

Cell[BoxData["1.1033574791548966`"], "Output",
 CellChangeTimes->{
  3.5670784509066887`*^9, 3.5670785564876738`*^9, 3.5670793607324867`*^9, 
   3.567079418989789*^9, 3.567079722519522*^9, {3.5670803370212016`*^9, 
   3.5670804052089214`*^9}, 3.567080471727438*^9, 3.567081090393325*^9, 
   3.567081122451381*^9, 3.5670811538698363`*^9, 3.5670812010599194`*^9, 
   3.5670812378915834`*^9, 3.567081274380048*^9, 3.5670824912601852`*^9, 
   3.5670828555988255`*^9, 3.567082903334909*^9, 3.5670836938662977`*^9, {
   3.5670844599392433`*^9, 3.567084518517346*^9}, 3.567084593257077*^9, 
   3.5670846788388276`*^9, 3.567085316575148*^9, 3.5670855076598835`*^9, 
   3.5670858936669617`*^9, 3.567085977657509*^9, 3.567086125623769*^9, 
   3.5670863498585625`*^9, 3.5670863799510155`*^9, 3.567086425393895*^9, {
   3.5670864623815603`*^9, 3.5670865208816633`*^9}, {3.567086552003718*^9, 
   3.567086592470189*^9}, 3.567086763118889*^9, 3.5670870740430346`*^9, 
   3.5670872060348663`*^9, 3.5670873213502693`*^9, 3.5670873552491283`*^9, 
   3.567087404139614*^9, {3.5670874361040707`*^9, 3.5670874603465133`*^9}, 
   3.5670874999549828`*^9, {3.5670979102628946`*^9, 3.5670979312449317`*^9}}],

Cell[BoxData["1.0077532890936167`"], "Output",
 CellChangeTimes->{
  3.5670784509066887`*^9, 3.5670785564876738`*^9, 3.5670793607324867`*^9, 
   3.567079418989789*^9, 3.567079722519522*^9, {3.5670803370212016`*^9, 
   3.5670804052089214`*^9}, 3.567080471727438*^9, 3.567081090393325*^9, 
   3.567081122451381*^9, 3.5670811538698363`*^9, 3.5670812010599194`*^9, 
   3.5670812378915834`*^9, 3.567081274380048*^9, 3.5670824912601852`*^9, 
   3.5670828555988255`*^9, 3.567082903334909*^9, 3.5670836938662977`*^9, {
   3.5670844599392433`*^9, 3.567084518517346*^9}, 3.567084593257077*^9, 
   3.5670846788388276`*^9, 3.567085316575148*^9, 3.5670855076598835`*^9, 
   3.5670858936669617`*^9, 3.567085977657509*^9, 3.567086125623769*^9, 
   3.5670863498585625`*^9, 3.5670863799510155`*^9, 3.567086425393895*^9, {
   3.5670864623815603`*^9, 3.5670865208816633`*^9}, {3.567086552003718*^9, 
   3.567086592470189*^9}, 3.567086763118889*^9, 3.5670870740430346`*^9, 
   3.5670872060348663`*^9, 3.5670873213502693`*^9, 3.5670873552491283`*^9, 
   3.567087404139614*^9, {3.5670874361040707`*^9, 3.5670874603465133`*^9}, 
   3.5670874999549828`*^9, {3.5670979102628946`*^9, 3.5670979312449317`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
plot1=Graphics[{Circle[{xc[[1]],yc[[1]]},Cr],      \
Circle[{xc[[2]],yc[[2]]},Nr],Circle[{xc[[3]],yc[[3]]},Oa],Circle[{xc[[4]],yc[[\
4]]},NH],
Circle[{xc[[5]],yc[[5]]},NH],Circle[{xc[[6]],yc[[6]]},CH],Circle[{xc[[7]],yc[[\
7]]},CN],Circle[{xc[[8]],yc[[8]]},CO],Circle[{xc[[9]],yc[[9]]},OL],Circle[{xc[\
[10]],yc[[10]]},NM],Circle[{xc[[11]],yc[[11]]},NM],Circle[{xc[[12]],yc[[12]]},\
OC],Circle[{xc[[13]],yc[[13]]},OC],Circle[{xc[[14]],yc[[14]]},OM],Circle[{xc[[\
15]],yc[[15]]},OM],Disk[{xn[[1]],yn[[1]]},0.08], \
Disk[{xn[[2]],yn[[2]]},0.08], \
Disk[{xn[[3]],yn[[3]]},0.08],Disk[{xn[[4]],yn[[4]]},0.08], \
Disk[{xn[[5]],yn[[5]]},0.08], Disk[{xn[[6]],yn[[6]]},0.08]} ] /. t[[2]];

Show[plot1,{AspectRatio \[Rule] Automatic,Axes -> True,GridLines -> \
Automatic, PlotRange \[Rule] Automatic, Frame -> True}]\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5668391703953342`*^9, 3.5668391835773573`*^9}, {
   3.5668412395117683`*^9, 3.566841239776969*^9}, {3.5668412898218565`*^9, 
   3.566841290087057*^9}, {3.56684186652367*^9, 3.566841909423745*^9}, 
   3.5670689312751684`*^9, {3.5670699737625995`*^9, 3.567070022403485*^9}, {
   3.5670703462132535`*^9, 3.5670703473676558`*^9}, {3.56707040660096*^9, 
   3.567070410220166*^9}, 3.567071022770842*^9, {3.567071058073704*^9, 
   3.567071077339738*^9}, {3.567071109569394*^9, 3.56707112992743*^9}, {
   3.567071166259894*^9, 3.5670711898627357`*^9}, {3.56707606019129*^9, 
   3.567076240699207*^9}, {3.567076605131447*^9, 3.5670766430707135`*^9}, {
   3.5670849912053766`*^9, 3.567085001423394*^9}}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{0, 0}, 0.2597231856403463], 
   CircleBox[{-2.5449706837000776`, 0}, 0.22405135815353058`], 
   CircleBox[{1.291429196139679, 1.884761711991131}, 0.1916417225081855], 
   CircleBox[{-3.200397178358574, 1.2212124156850317`}, 1.1619297579240278`], 
   CircleBox[{-3.200397178358574, -1.2212124156850317`}, 1.1619297579240278`],
    CircleBox[{0.6081672572221739, -1.3853495355430014`}, 
    1.2532410401053755`], 
   CircleBox[{-1.2903212555934467`, 0}, 1.0305980699531003`], 
   CircleBox[{0.6649556762928184, 0.9704620296599826}, 0.9166962134031293], 
   CircleBox[{1.9046760746972449`, 2.77975792251382}, 0.8932960411034288], 
   CircleBox[{-2.5449706837000776`, 0}, 0.7472908218886891], 
   CircleBox[{-2.5449706837000776`, 0}, 0.7472908218886891], 
   CircleBox[{0.6649556762928184, 0.9704620296599826}, 1.5679221435272637`], 
   CircleBox[{0.6649556762928184, 0.9704620296599826}, 1.5679221435272637`], 
   CircleBox[{0.44483311166516193`, 2.4648450688802646`}, 0.8346241451827114],
    CircleBox[{2.138025280614196, 1.3046783551019976`}, 0.8346241451827114], 
   DiskBox[{0, 0}, 0.08], DiskBox[{-2.5449706837000776`, 0}, 0.08], 
   DiskBox[{1.291429196139679, 1.884761711991131}, 0.08], 
   DiskBox[{-3.4455458033171533`, 1.6779814766359766`}, 0.08], 
   DiskBox[{-3.4455458033171533`, -1.6779814766359766`}, 0.08], 
   DiskBox[{0.8381266667787244, -1.9091760938454152`}, 0.08]},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->Automatic]], "Output",
 CellChangeTimes->{
  3.5670784509222884`*^9, 3.567078556518874*^9, 3.5670793607636867`*^9, 
   3.567079419005389*^9, 3.567079722535122*^9, {3.5670803370368013`*^9, 
   3.567080405224521*^9}, 3.567080471743038*^9, 3.567081090408925*^9, 
   3.567081122482581*^9, 3.5670811539010363`*^9, 3.5670812010911193`*^9, 
   3.567081237922784*^9, 3.567081274395648*^9, 3.567082491291385*^9, 
   3.567082855630025*^9, 3.567082903366109*^9, 3.5670836938974977`*^9, {
   3.5670844599704432`*^9, 3.5670845185485463`*^9}, 3.5670845932882776`*^9, 
   3.5670846788700275`*^9, 3.567085316606348*^9, 3.567085507691084*^9, 
   3.5670858936981616`*^9, 3.5670859776887093`*^9, 3.567086125654969*^9, 
   3.567086349889763*^9, 3.567086379982216*^9, 3.5670864254250956`*^9, {
   3.5670864624127607`*^9, 3.5670865209128633`*^9}, {3.567086552050518*^9, 
   3.5670865925013885`*^9}, 3.567086763150089*^9, 3.567087074074235*^9, 
   3.5670872060660667`*^9, 3.5670873213814692`*^9, 3.567087355295929*^9, 
   3.5670874041708145`*^9, {3.567087436150871*^9, 3.567087460377713*^9}, 
   3.567087500001783*^9, {3.567097910294095*^9, 3.5670979312449317`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
plot2=Graphics[{Circle[{xc[[1]],zc[[1]]},Cr],      \
Circle[{xc[[2]],zc[[2]]},Nr],Circle[{xc[[3]],zc[[3]]},Oa],Circle[{xc[[4]],zc[[\
4]]},NH],
Circle[{xc[[5]],zc[[5]]},NH],Circle[{xc[[6]],zc[[6]]},CH],Circle[{xc[[7]],zc[[\
7]]},CN],Circle[{xc[[8]],zc[[8]]},CO],Circle[{xc[[9]],zc[[9]]},OL],Circle[{xc[\
[10]],zc[[10]]},NM],Circle[{xc[[11]],zc[[11]]},NM],Circle[{xc[[12]],zc[[12]]},\
OC],Circle[{xc[[13]],zc[[13]]},OC],Circle[{xc[[14]],zc[[14]]},OM],Circle[{xc[[\
15]],zc[[15]]},OM],Disk[{xn[[1]],zn[[1]]},0.08], \
Disk[{xn[[2]],zn[[2]]},0.08], \
Disk[{xn[[3]],zn[[3]]},0.08],Disk[{xn[[4]],zn[[4]]},0.08], \
Disk[{xn[[5]],zn[[5]]},0.08], Disk[{xn[[6]],zn[[6]]},0.08]} ] /. t[[2]];

Show[plot2,{AspectRatio \[Rule] Automatic,Axes -> True,GridLines -> \
Automatic, PlotRange \[Rule] Automatic, Frame -> True}]\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5668391703953342`*^9, 3.5668391881949654`*^9}, {
   3.5668412398393693`*^9, 3.5668412401201696`*^9}, {3.566841311459095*^9, 
   3.5668413117242956`*^9}, {3.5668419284557786`*^9, 
   3.5668419511070185`*^9}, {3.566842124953724*^9, 3.566842139087349*^9}, 
   3.5670689313531685`*^9, {3.5670705170023537`*^9, 3.567070548654809*^9}, {
   3.567076321179748*^9, 3.567076324596154*^9}, {3.5670763612250185`*^9, 
   3.567076444107964*^9}, {3.5670764854792366`*^9, 3.5670764972884574`*^9}, {
   3.567076691212398*^9, 3.5670767228336535`*^9}, {3.5670850083498063`*^9, 
   3.5670850144962173`*^9}}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{0, 0}, 0.2597231856403463], 
   CircleBox[{-2.5449706837000776`, 0}, 0.22405135815353058`], 
   CircleBox[{1.291429196139679, 0}, 0.1916417225081855], 
   CircleBox[{-3.200397178358574, 0}, 1.1619297579240278`], 
   CircleBox[{-3.200397178358574, 0}, 1.1619297579240278`], 
   CircleBox[{0.6081672572221739, 0}, 1.2532410401053755`], 
   CircleBox[{-1.2903212555934467`, 0}, 1.0305980699531003`], 
   CircleBox[{0.6649556762928184, 0}, 0.9166962134031293], 
   CircleBox[{1.9046760746972449`, 0}, 0.8932960411034288], 
   CircleBox[{-2.5449706837000776`, 0.7472908218886891}, 0.7472908218886891], 
   CircleBox[{-2.5449706837000776`, -0.7472908218886891}, 0.7472908218886891],
    CircleBox[{0.6649556762928184, 1.5679221435272637`}, 1.5679221435272637`],
    CircleBox[{0.6649556762928184, -1.5679221435272637`}, 
    1.5679221435272637`], 
   CircleBox[{0.44483311166516193`, 0}, 0.8346241451827114], 
   CircleBox[{2.138025280614196, 0}, 0.8346241451827114], 
   DiskBox[{0, 0}, 0.08], DiskBox[{-2.5449706837000776`, 0}, 0.08], 
   DiskBox[{1.291429196139679, 0}, 0.08], 
   DiskBox[{-3.4455458033171533`, 0}, 0.08], 
   DiskBox[{-3.4455458033171533`, 0}, 0.08], 
   DiskBox[{0.8381266667787244, 0}, 0.08]},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->Automatic]], "Output",
 CellChangeTimes->{
  3.5670784509378886`*^9, 3.567078556550074*^9, 3.567079360779287*^9, 
   3.567079419036589*^9, 3.5670797225663223`*^9, {3.5670803370524015`*^9, 
   3.5670804052401214`*^9}, 3.5670804717586384`*^9, 3.567081090424525*^9, 
   3.567081122482581*^9, 3.5670811539166365`*^9, 3.5670812010911193`*^9, 
   3.567081237938384*^9, 3.567081274411248*^9, 3.567082491291385*^9, 
   3.5670828556612253`*^9, 3.5670829033973093`*^9, 3.5670836939130974`*^9, {
   3.567084460001643*^9, 3.5670845185797462`*^9}, 3.567084593303878*^9, 
   3.5670846788856277`*^9, 3.5670853166219482`*^9, 3.5670855077222834`*^9, 
   3.5670858937137613`*^9, 3.567085977719909*^9, 3.567086125670569*^9, 
   3.567086349905363*^9, 3.5670863799978156`*^9, 3.567086425440696*^9, {
   3.5670864624127607`*^9, 3.567086520928463*^9}, {3.567086552050518*^9, 
   3.5670865925169888`*^9}, 3.5670867631656885`*^9, 3.5670870740898347`*^9, 
   3.567087206081667*^9, 3.567087321397069*^9, 3.567087355295929*^9, 
   3.5670874041708145`*^9, {3.567087436150871*^9, 3.5670874603933134`*^9}, 
   3.567087500001783*^9, {3.567097910309695*^9, 3.5670979312605314`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
plot3=Graphics3D[{
Sphere[{xc[[1]],yc[[1]],zc[[1]]},Cr],      
Sphere[{xc[[2]],yc[[2]],zc[[2]]},Nr],
Sphere[{xc[[3]],yc[[3]],zc[[3]]},Oa],
{Opacity[0.6],
Sphere[{xc[[4]],yc[[4]],zc[[4]]},NH],
Sphere[{xc[[5]],yc[[5]],zc[[5]]},NH],
Sphere[{xc[[6]],yc[[6]],zc[[6]]},CH],
Sphere[{xc[[7]],yc[[7]],zc[[7]]},CN],
Sphere[{xc[[8]],yc[[8]],zc[[8]]},CO],
{Lighter[Blue],
Sphere[{xc[[9]],yc[[9]],zc[[9]]},OL],
Sphere[{xc[[10]],yc[[10]],zc[[10]]},NM],
Sphere[{xc[[11]],yc[[11]],zc[[11]]},NM],
Sphere[{xc[[14]],yc[[14]],zc[[14]]},OM],
Sphere[{xc[[15]],yc[[15]],zc[[15]]},OM]}},
{Opacity[0.4],{Lighter[Yellow],
Sphere[{xc[[12]],yc[[12]],zc[[12]]},OC],
Sphere[{xc[[13]],yc[[13]],zc[[13]]},OC]}},
{Darker[Red,1],
Sphere[{xn[[1]],yc[[1]],zn[[1]]},0.08],
Sphere[{xn[[2]],yc[[2]],zn[[2]]},0.08],
Sphere[{xn[[3]],yc[[3]],zn[[3]]},0.08],
Sphere[{xn[[4]],yc[[4]],zn[[4]]},0.08],
Sphere[{xn[[5]],yc[[5]],zn[[5]]},0.08],
Sphere[{xn[[6]],yc[[6]],zn[[6]]},0.08]}}] /. t[[2]];

Show[plot3,{AspectRatio \[Rule] Automatic,Boxed->False}]\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5668391703953342`*^9, 3.566839177774147*^9}, {
   3.5669295568515253`*^9, 3.5669295979107976`*^9}, {3.5669296458028817`*^9, 
   3.5669298487904377`*^9}, {3.5669298934845166`*^9, 
   3.5669299137177525`*^9}, {3.566929958552231*^9, 3.566929984853877*^9}, {
   3.5669300180975356`*^9, 3.5669300244935465`*^9}, {3.5669681367168245`*^9, 
   3.566968174578091*^9}, {3.5669682084145503`*^9, 3.5669682238429775`*^9}, {
   3.56696839045127*^9, 3.566968403352493*^9}, {3.5669886931259174`*^9, 
   3.5669887406436005`*^9}, {3.56698877437086*^9, 3.5669887786296673`*^9}, {
   3.5669888178481364`*^9, 3.5669888516221952`*^9}, {3.5669888927126675`*^9, 
   3.5669888937422695`*^9}, {3.56698893940355*^9, 3.566988940932352*^9}, {
   3.566989011850077*^9, 3.5669890235968976`*^9}, {3.5669890790549955`*^9, 
   3.566989179831172*^9}, {3.5669892163196363`*^9, 3.5669892198764424`*^9}, {
   3.566989533171793*^9, 3.5669895449498134`*^9}, 3.567068931415569*^9, {
   3.5670704486430335`*^9, 3.567070500840725*^9}, {3.5670768428602643`*^9, 
   3.567076863427901*^9}, {3.567076934439225*^9, 3.5670771373487816`*^9}, {
   3.567077191402877*^9, 3.567077229716544*^9}, {3.5670776494196815`*^9, 
   3.5670776603397007`*^9}, 3.567077706489382*^9, 3.5670779441365986`*^9, {
   3.5670781500413604`*^9, 3.567078154253368*^9}, {3.567078187263026*^9, 
   3.5670782620127573`*^9}, {3.5670783597557287`*^9, 3.567078441453072*^9}, {
   3.5670784741819296`*^9, 3.5670785488436604`*^9}, {3.567079302263584*^9, 
   3.567079352339672*^9}, 3.5670794136233797`*^9, {3.567084448847624*^9, 
   3.567084513977738*^9}, {3.567084586159065*^9, 3.5670845877970676`*^9}, {
   3.567084647170772*^9, 3.5670846660780053`*^9}, {3.5670850221714306`*^9, 
   3.567085027023039*^9}}],

Cell[BoxData[
 Graphics3DBox[{SphereBox[{0, 0, 0}, 0.2597231856403463], 
   SphereBox[{-2.5449706837000776`, 0, 0}, 0.22405135815353058`], 
   SphereBox[{1.291429196139679, 1.884761711991131, 0}, 0.1916417225081855], 
   {Opacity[0.6], 
    SphereBox[{-3.200397178358574, 1.2212124156850317`, 0}, 
     1.1619297579240278`], 
    SphereBox[{-3.200397178358574, -1.2212124156850317`, 0}, 
     1.1619297579240278`], 
    SphereBox[{0.6081672572221739, -1.3853495355430014`, 0}, 
     1.2532410401053755`], 
    SphereBox[{-1.2903212555934467`, 0, 0}, 1.0305980699531003`], 
    SphereBox[{0.6649556762928184, 0.9704620296599826, 0}, 
     0.9166962134031293], 
    {RGBColor[
      NCache[
       Rational[1, 3], 0.3333333333333333], 
      NCache[
       Rational[1, 3], 0.3333333333333333], 1], 
     SphereBox[{1.9046760746972449`, 2.77975792251382, 0}, 
      0.8932960411034288], 
     SphereBox[{-2.5449706837000776`, 0, 0.7472908218886891}, 
      0.7472908218886891], 
     SphereBox[{-2.5449706837000776`, 0, -0.7472908218886891}, 
      0.7472908218886891], 
     SphereBox[{0.44483311166516193`, 2.4648450688802646`, 0}, 
      0.8346241451827114], 
     SphereBox[{2.138025280614196, 1.3046783551019976`, 0}, 
      0.8346241451827114]}}, 
   {Opacity[0.4], 
    {RGBColor[1, 1, 
      NCache[
       Rational[1, 3], 0.3333333333333333]], 
     SphereBox[{0.6649556762928184, 0.9704620296599826, 1.5679221435272637`}, 
      1.5679221435272637`], 
     SphereBox[{0.6649556762928184, 0.9704620296599826, -1.5679221435272637`},
       1.5679221435272637`]}}, 
   {RGBColor[0, 0, 0], SphereBox[{0, 0, 0}, 0.08], 
    SphereBox[{-2.5449706837000776`, 0, 0}, 0.08], 
    SphereBox[{1.291429196139679, 1.884761711991131, 0}, 0.08], 
    SphereBox[{-3.4455458033171533`, 1.2212124156850317`, 0}, 0.08], 
    SphereBox[{-3.4455458033171533`, -1.2212124156850317`, 0}, 0.08], 
    SphereBox[{0.8381266667787244, -1.3853495355430014`, 0}, 0.08]}},
  AspectRatio->Automatic,
  Boxed->False,
  ImageSize->{463.76936334873824`, 466.},
  ViewPoint->{-0.8360838006883977, -1.8692398322265755`, 2.6938645711772566`},
  ViewVertical->{0.07609704531206406, 0.34258655569640434`, 
   1.3537210602769851`}]], "Output",
 CellChangeTimes->{
  3.5670784509534883`*^9, 3.567078556565674*^9, 3.5670793608572865`*^9, 
   3.56707941927059*^9, 3.5670797225819225`*^9, {3.5670803370680017`*^9, 
   3.567080405255721*^9}, 3.567080471774238*^9, 3.5670810904401245`*^9, 
   3.5670811224981813`*^9, 3.5670811539166365`*^9, 3.567081201106719*^9, 
   3.567081237938384*^9, 3.567081274426848*^9, 3.5670824913069854`*^9, 
   3.5670828556768255`*^9, 3.5670829034129095`*^9, 3.5670836939286976`*^9, {
   3.5670844600172434`*^9, 3.5670845185953465`*^9}, 3.5670845933194776`*^9, 
   3.567084678901228*^9, 3.567085316637548*^9, 3.5670855077378836`*^9, 
   3.5670858937293615`*^9, 3.567085977735509*^9, 3.567086125701769*^9, 
   3.567086349920963*^9, 3.5670863800134153`*^9, 3.5670864254562955`*^9, {
   3.56708646244396*^9, 3.567086520944063*^9}, {3.5670865520661182`*^9, 
   3.567086592532589*^9}, 3.567086763196889*^9, 3.5670870740898347`*^9, 
   3.567087206097267*^9, 3.567087321412669*^9, 3.5670873553115287`*^9, 
   3.567087404202015*^9, {3.5670874361664705`*^9, 3.567087460408913*^9}, 
   3.5670875000173826`*^9, {3.5670979103252945`*^9, 3.5670979312761316`*^9}}]
}, Open  ]]
},
WindowToolbars->"EditBar",
WindowSize->{830, 964},
WindowMargins->{{Automatic, 525}, {Automatic, 0}},
PrintingCopies->1,
PrintingPageRange->{5, 6},
PrintingOptions->{"Magnification"->1.,
"PaperOrientation"->"Portrait",
"PaperSize"->{595.2, 841.8}},
PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}},
FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (November 20, 2012)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[557, 20, 2599, 38, 48, "Input"],
Cell[3159, 60, 3877, 59, 218, "Input"],
Cell[7039, 121, 1319, 26, 72, "Input"],
Cell[8361, 149, 758, 23, 252, "Input"],
Cell[9122, 174, 4686, 67, 82, "Input"],
Cell[13811, 243, 4863, 70, 150, "Input"],
Cell[18677, 315, 5053, 71, 167, "Input"],
Cell[23733, 388, 4639, 66, 82, "Input"],
Cell[28375, 456, 4771, 67, 82, "Input"],
Cell[33149, 525, 4617, 66, 82, "Input"],
Cell[37769, 593, 489, 12, 99, "Input"],
Cell[38261, 607, 384, 11, 99, "Input"],
Cell[38648, 620, 571, 14, 116, "Input"],
Cell[CellGroupData[{
Cell[39244, 638, 1756, 26, 116, "Input"],
Cell[41003, 666, 1721, 31, 72, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[42761, 702, 955, 24, 235, "Input"],
Cell[43719, 728, 1181, 16, 31, "Output"],
Cell[44903, 746, 1197, 17, 31, "Output"],
Cell[46103, 765, 1181, 16, 31, "Output"],
Cell[47287, 783, 1179, 16, 31, "Output"],
Cell[48469, 801, 1183, 16, 31, "Output"],
Cell[49655, 819, 1183, 16, 31, "Output"],
Cell[50841, 837, 1183, 16, 31, "Output"],
Cell[52027, 855, 1182, 16, 31, "Output"],
Cell[53212, 873, 1183, 16, 31, "Output"],
Cell[54398, 891, 1182, 16, 31, "Output"],
Cell[55583, 909, 1183, 16, 31, "Output"],
Cell[56769, 927, 1183, 16, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[57989, 948, 1572, 26, 218, "Input"],
Cell[59564, 976, 2659, 42, 328, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[62260, 1023, 1466, 25, 218, "Input"],
Cell[63729, 1050, 2491, 42, 327, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[66257, 1097, 2803, 54, 490, "Input"],
Cell[69063, 1153, 3339, 64, 481, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)