Wall clock time and date at job start Wed Apr 1 2020 05:18:53
AMSOL-version 7.1
by G. D. Hawkins, D. J. Giesen, G. C. Lynch, C. C. Chambers,
I. Rossi, J. W. Storer, J. Li, J. D. Thompson, P. Winget,
B. J. Lynch, D. Rinaldi, D. A. Liotard, C. J. Cramer, and D. G. Truhlar
Copyright 2004 by Regents of the University of Minnesota.
All rights reserved.
Notice: recipients of this code are asked to comply
with the user agreement in Section 1 of the documentation file.
*******************************************************************************
* 1SCF - SCF CALCULATION WITHOUT GEOMETRY OPTIMIZATION
* - USE EF ROUTINE FOR MINIMUM SEARCH (DEFAULT)
* GEO-OK - OVERRIDE INTERATOMIC DISTANCE CHECK
* TLIMIT= - A TIME OF 15. SECONDS REQUESTED
* CHARGE - CHARGE ON SYSTEM= -1
* AM1 - THE AM1 HAMILTONIAN TO BE USED
* DEV - DEVELOPER OPTIONS ARE ALLOWED
* SM5.42R - SM5.42R CALCULATIONS WILL BE PERFORMED
* - CHARGE MODEL 2 WILL BE USED (CM2A)
* - THE SOLVENT IS USER-SPECIFIED
* DIELEC - THE SOLVENT DIELECTRIC CONSTANT IS 2.06
* IOFR - THE SOLVENT INDEX OF REFRACTION IS 1.4345
* ALPHA - THE SOLVENT ALPHA IS 0.00
* BETA - THE SOLVENT BETA IS 0.00
* GAMMA - THE SOLVENT MACROSCOPIC SURFACE TENSION IS
* 38.93 CAL MOL^-1 ANGSTROM^-2
* FACARB - THE FRACTION OF AROMATIC CARBONS IS 0.00
* FEHALO - THE FRACTION OF ELECTRONEGATIVE HALOGENS IS 0.00
*******************************************************************************
Atom NO. Chemical Bond length Bond angle Dihedral angle
number (k) symbol (angstroms) (degrees) (degrees)
(I) NA:I NB:NA:I NC:NB:NA:I NA NB NC
1 1 C
2 2 C 1.53005 * 1
3 3 H 1.08999 * 109.47159 * 2 1
4 4 C 1.50703 * 109.46659 * 119.99425 * 2 1 3
5 5 O 1.21288 * 120.00006 * 94.20799 * 4 2 1
6 6 N 1.34773 * 120.00287 * 274.20124 * 4 2 1
7 7 C 1.47418 * 125.65162 * 174.69382 * 6 4 2
8 8 C 1.54908 * 104.83348 * 204.16688 * 7 6 4
9 9 H 1.08998 * 111.01010 * 81.06404 * 8 7 6
10 10 C 1.50706 * 111.00316 * 204.91255 * 8 7 6
11 11 C 1.34674 * 125.98871 * 26.91437 * 10 8 7
12 12 N 1.36884 * 106.81351 * 179.84236 * 11 10 8
13 13 C 1.46502 * 126.39146 * 180.11556 * 12 11 10
14 14 C 1.35037 * 107.20623 * 0.41341 * 12 11 10
15 15 N 1.30497 * 108.70125 * 359.73916 * 14 12 11
16 16 C 1.55151 * 101.58224 * 322.98770 * 8 7 6
17 17 H 1.09005 * 110.87309 * 277.07198 * 16 8 7
18 18 N 1.46503 * 110.72170 * 153.86563 * 16 8 7
19 19 C 1.36940 * 119.99607 * 160.02521 * 18 16 8
20 20 H 1.07997 * 120.00481 * 0.02562 * 19 18 16
21 21 C 1.38819 * 119.99655 * 180.02562 * 19 18 16
22 22 N 1.31619 * 119.99462 * 359.97438 * 21 19 18
23 23 Si 1.86798 * 120.00055 * 179.72444 * 21 19 18
24 24 H 1.48499 * 109.46992 * 359.97438 * 23 21 19
25 25 H 1.48496 * 109.47086 * 119.99395 * 23 21 19
26 26 H 1.48497 * 109.46847 * 239.99595 * 23 21 19
27 27 C 1.47025 * 125.64644 * 354.73622 * 6 4 2
28 28 C 1.52991 * 109.47311 * 239.99348 * 2 1 3
29 29 C 1.53781 * 113.61501 * 159.13690 * 28 2 1
30 30 C 1.53953 * 87.02146 * 220.04951 * 29 28 2
31 31 C 1.53776 * 87.11312 * 25.36408 * 30 29 28
32 32 H 1.08996 * 109.46823 * 64.82310 * 1 2 3
33 33 H 1.08997 * 109.47308 * 184.81767 * 1 2 3
34 34 H 1.08994 * 109.47215 * 304.81963 * 1 2 3
35 35 H 1.09000 * 110.36428 * 323.00333 * 7 6 4
36 36 H 1.08996 * 110.36926 * 85.21451 * 7 6 4
37 37 H 1.07994 * 126.59075 * 359.96047 * 11 10 8
38 38 H 1.09003 * 109.47156 * 270.32126 * 13 12 11
39 39 H 1.08998 * 109.47256 * 30.32722 * 13 12 11
40 40 H 1.09005 * 109.46501 * 150.32716 * 13 12 11
41 41 H 1.08008 * 125.65016 * 179.71539 * 14 12 11
42 42 H 0.96994 * 119.99793 * 340.02877 * 18 16 8
43 43 H 0.97008 * 119.99276 * 179.97438 * 22 21 19
44 44 H 1.09001 * 109.88186 * 59.48671 * 27 6 4
45 45 H 1.09003 * 109.88180 * 298.32147 * 27 6 4
46 46 H 1.09004 * 112.85459 * 290.33477 * 28 2 1
47 47 H 1.08997 * 113.61621 * 334.57111 * 29 28 2
48 48 H 1.08998 * 113.61468 * 105.47262 * 29 28 2
49 49 H 1.08996 * 113.54043 * 270.80943 * 30 29 28
50 50 H 1.09002 * 113.74511 * 139.94275 * 30 29 28
51 51 H 1.09001 * 113.61828 * 220.09190 * 31 30 29
52 52 H 1.08995 * 113.69443 * 89.19363 * 31 30 29
Note: An asterisk (*) indicates this parameter was or will be optimized.
Cartesian coordinates (angstroms)
NO. Atom X Y Z
1 6 0.0000 0.0000 0.0000
2 6 1.5300 0.0000 0.0000
3 1 1.8934 1.0276 0.0000
4 6 2.0323 -0.7103 1.2306
5 8 2.3070 -1.8906 1.1799
6 7 2.1762 -0.0341 2.3875
7 6 2.5493 -0.6231 3.6863
8 6 1.9890 0.3607 4.7436
9 1 0.9257 0.1898 4.9116
10 6 2.7601 0.2630 6.0347
11 6 4.0433 -0.1254 6.1626
12 7 4.3436 -0.0667 7.4968
13 6 5.6306 -0.3978 8.1135
14 6 3.2300 0.3467 8.1391
15 7 2.2871 0.5420 7.2583
16 6 2.2301 1.7295 4.0540
17 1 3.2532 2.0699 4.2137
18 7 1.2724 2.7282 4.5353
19 6 1.5425 4.0626 4.3877
20 1 2.4612 4.3787 3.9162
21 6 0.6353 5.0089 4.8442
22 7 -0.4847 4.6236 5.4184
23 14 1.0071 6.8292 4.6499
24 1 2.3234 6.9971 3.9833
25 1 -0.0509 7.4652 3.8246
26 1 1.0447 7.4722 5.9879
27 6 1.9864 1.4132 2.5633
28 6 2.0401 -0.7213 -1.2491
29 6 3.4868 -0.3577 -1.6229
30 6 3.0235 -0.3622 -3.0910
31 6 1.6122 -0.0511 -2.5654
32 1 -0.3633 0.4372 -0.9300
33 1 -0.3634 -1.0240 0.0863
34 1 -0.3633 0.5868 0.8436
35 1 2.0920 -1.6057 3.8026
36 1 3.6334 -0.6976 3.7718
37 1 4.7082 -0.4236 5.3656
38 1 5.6479 -1.4564 8.3725
39 1 6.4360 -0.1823 7.4113
40 1 5.7655 0.1992 9.0155
41 1 3.1353 0.4915 9.2053
42 1 0.4471 2.4443 4.9584
43 1 -1.1184 5.2849 5.7377
44 1 0.9689 1.6898 2.2872
45 1 2.7012 1.9584 1.9469
46 1 1.8521 -1.7948 -1.2245
47 1 3.7964 0.6244 -1.2657
48 1 4.2080 -1.1409 -1.3894
49 1 3.1057 -1.3375 -3.5705
50 1 3.4534 0.4411 -3.6894
51 1 0.8223 -0.5870 -3.0916
52 1 1.4069 1.0147 -2.4660
RHF calculation, no. of doubly occupied orbitals= 66
REFERENCES FOR PARAMETERS IN GAS-PHASE HAMILTONIAN:
H: (AM1): M.J.S. DEWAR ET AL, J. AM. CHEM. SOC. 107 3902-3909 (1985)
C: (AM1): M.J.S. DEWAR ET AL, J. AM. CHEM. SOC. 107 3902-3909 (1985)
N: (AM1): M.J.S. DEWAR ET AL, J. AM. CHEM. SOC. 107 3902-3909 (1985)
O: (AM1): M.J.S. DEWAR ET AL, J. AM. CHEM. SOC. 107 3902-3909 (1985)
Si: (AM1): M.J.S.DEWAR, C. JIE, ORGANOMETALLICS, 6, 1486-1490 (1987).
REFERENCE FOR CHARGE MODEL 2:
J. Li, J. Xing, C. J. Cramer, and D. G. Truhlar, J. Chem. Phys. 111 (1999) 885.
REFERENCE FOR THE PARAMETERIZATION OF THE SM5.42R SOLVATION MODEL:
J. Li, T. Zhu, G. D. Hawkins, P. Winget, D. A. Liotard, C. J. Cramer, and
D. G. Truhlar, Theor. Chem. Acc. 103 (1999) 9-63
P. Winget, J. D. Thompson C. J. Cramer,and D. G. Truhlar, J. Phys. Chem. B. submitted.
CHARGE=-1 AM1 1SCF TLIMIT=15 GEO-OK SM5.42R
& SOLVNT=GENORG IOFR=1.4345 ALPHA=0.00 BETA=0.00 GAMMA=38.93
& DIELEC=2.06 FACARB=0.00 FEHALO=0.00 DEV
ZINC001070505771.mol2 52
1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED
SCF FIELD WAS ACHIEVED
Wall clock time and date at job start Wed Apr 1 2020 05:18:53
Heat of formation + Delta-G solvation = 61.845233 kcal
Electronic energy + Delta-G solvation = -31299.143978 eV
Core-core repulsion = 27300.697670 eV
Total energy + Delta-G solvation = -3998.446307 eV
No. of doubly occupied orbitals = 66
Molecular weight
(most abundant/longest-lived isotopes) = 346.211 amu
Computer time = 1.21 seconds
Orbital eigenvalues (eV)
-41.84422 -40.20164 -38.65087 -37.17268 -34.06578 -33.90277 -32.94002 -31.51793
-30.57585 -28.26517 -27.37199 -26.74156 -26.16890 -25.26650 -24.31619 -23.34359
-22.02540 -20.94561 -20.67887 -19.49514 -19.33895 -18.67186 -17.38688 -16.89314
-16.64478 -15.78928 -15.62092 -15.47739 -15.13758 -14.66592 -14.51874 -14.17157
-14.01259 -13.70786 -13.39761 -13.10821 -13.04649 -12.91812 -12.82618 -12.71008
-12.58754 -12.53702 -12.41921 -11.93259 -11.86313 -11.64252 -11.48573 -11.33905
-11.08318 -10.84059 -10.75997 -10.52188 -10.39982 -10.35037 -10.14163 -10.03843
-9.73030 -9.54558 -9.43764 -9.26376 -8.94796 -8.27481 -8.06126 -6.97759
-5.75590 -3.05590 2.25904 2.89339 2.95280 3.26422 3.42644 3.46550
3.48396 3.50640 4.46413 4.56181 4.58766 4.62804 4.65768 4.73935
4.78790 4.95091 4.97903 5.00240 5.16809 5.20760 5.21464 5.26164
5.40155 5.43810 5.45267 5.62819 5.64716 5.66450 5.69190 5.74142
5.77680 5.80962 5.84327 5.87715 5.90517 5.96594 5.99779 6.05546
6.20557 6.24995 6.31628 6.39375 6.43517 6.48487 6.50878 6.63954
6.90210 7.01675 7.39714 7.69685 7.86826 8.40669 8.43166 8.46030
9.53455 10.08737 10.44312 11.46113
Molecular weight = 346.21amu
Principal moments of inertia in cm(-1)
A = 0.007098 B = 0.003884 C = 0.002744
Principal moments of inertia in units of 10**(-40)*gram-cm**2
A = 3943.632266 B = 7207.365836 C =10201.101907
Net atomic charges, atomic populations, and dipole contributions
calculated with CM2
Atom NO. Type Charge No. of electrons
1 C -0.139 4.139
2 C -0.097 4.097
3 H 0.091 0.909
4 C 0.535 3.465
5 O -0.555 6.555
6 N -0.612 5.612
7 C 0.120 3.880
8 C -0.027 4.027
9 H 0.096 0.904
10 C 0.040 3.960
11 C -0.066 4.066
12 N -0.468 5.468
13 C 0.076 3.924
14 C 0.177 3.823
15 N -0.468 5.468
16 C 0.187 3.813
17 H 0.063 0.937
18 N -0.705 5.705
19 C -0.244 4.244
20 H 0.070 0.930
21 C -0.006 4.006
22 N -0.924 5.924
23 Si 0.903 3.097
24 H -0.281 1.281
25 H -0.291 1.291
26 H -0.292 1.292
27 C 0.072 3.928
28 C -0.091 4.091
29 C -0.138 4.138
30 C -0.142 4.142
31 C -0.135 4.135
32 H 0.067 0.933
33 H 0.058 0.942
34 H 0.065 0.935
35 H 0.074 0.926
36 H 0.064 0.936
37 H 0.170 0.830
38 H 0.072 0.928
39 H 0.082 0.918
40 H 0.085 0.915
41 H 0.201 0.799
42 H 0.390 0.610
43 H 0.255 0.745
44 H 0.077 0.923
45 H 0.059 0.941
46 H 0.099 0.901
47 H 0.079 0.921
48 H 0.066 0.934
49 H 0.077 0.923
50 H 0.065 0.935
51 H 0.066 0.934
52 H 0.082 0.918
For the charges calculated by CM2:
Dipole moment (debyes) X Y Z Total
from point charges 10.496 -11.976 -5.548 16.863
Note: The Mulliken population analysis charges presented below were not
used in the solvation calculation but are provided for completeness.
The chosen solvation model uses CM2 partial charges for calculating
solvation energies.
Net atomic charges, atomic populations, and dipole contributions
using Mulliken population analysis
Atom NO. Type Charge No. of electrons
1 C -0.197 4.197
2 C -0.118 4.118
3 H 0.109 0.891
4 C 0.324 3.676
5 O -0.433 6.433
6 N -0.345 5.345
7 C -0.003 4.003
8 C -0.048 4.048
9 H 0.114 0.886
10 C -0.090 4.090
11 C -0.198 4.198
12 N -0.152 5.152
13 C -0.067 4.067
14 C -0.105 4.105
15 N -0.185 5.185
16 C 0.074 3.926
17 H 0.080 0.920
18 N -0.345 5.345
19 C -0.375 4.375
20 H 0.088 0.912
21 C -0.203 4.203
22 N -0.571 5.571
23 Si 0.734 3.266
24 H -0.206 1.206
25 H -0.219 1.219
26 H -0.219 1.219
27 C -0.054 4.054
28 C -0.110 4.110
29 C -0.175 4.175
30 C -0.179 4.179
31 C -0.173 4.173
32 H 0.086 0.914
33 H 0.077 0.923
34 H 0.084 0.916
35 H 0.092 0.908
36 H 0.082 0.918
37 H 0.187 0.813
38 H 0.091 0.909
39 H 0.100 0.900
40 H 0.103 0.897
41 H 0.218 0.782
42 H 0.225 0.775
43 H 0.065 0.935
44 H 0.095 0.905
45 H 0.078 0.922
46 H 0.117 0.883
47 H 0.098 0.902
48 H 0.084 0.916
49 H 0.096 0.904
50 H 0.084 0.916
51 H 0.085 0.915
52 H 0.101 0.899
Dipole moment (debyes) X Y Z Total
from point charges 10.028 -12.196 -5.396 16.686
hybrid contribution 0.846 0.891 0.107 1.234
sum 10.875 -11.305 -5.289 16.554
Atomic orbital electron populations
1.21670 0.92495 1.02595 1.02923 1.21272 0.97802 1.00325 0.92383
0.89080 1.20301 0.77099 0.86706 0.83514 1.90635 1.51693 1.16203
1.84818 1.48538 1.68000 1.09994 1.08013 1.22108 0.99323 0.96981
0.81925 1.20823 0.99535 0.96178 0.88287 0.88566 1.20093 0.91843
1.03823 0.93213 1.22345 0.99039 1.12883 0.85519 1.45350 1.11821
1.50883 1.07189 1.21959 0.83412 1.01824 0.99478 1.24974 0.81906
1.03948 0.99700 1.70215 1.29907 1.18917 0.99461 1.21143 0.95026
0.87813 0.88590 0.91950 1.42298 1.21159 1.01310 1.69684 1.19228
1.02947 0.82435 1.32927 0.91220 1.24726 0.96821 0.97766 1.01001
1.69688 1.10183 1.32558 1.44704 0.86185 0.77315 0.85839 0.77213
1.20617 1.21851 1.21919 1.23412 1.01894 0.78311 1.01816 1.22433
0.95950 1.00543 0.92075 1.22985 0.96736 1.01772 0.96052 1.22808
0.97008 1.01608 0.96487 1.22949 0.96885 1.01041 0.96378 0.91409
0.92287 0.91591 0.90763 0.91777 0.81311 0.90921 0.89962 0.89677
0.78243 0.77531 0.93520 0.90479 0.92236 0.88258 0.90248 0.91560
0.90416 0.91607 0.91530 0.89931
Number of geometries 1
Number of calculations of the screened coulomb radii 1
The total number of SCF iterations 14.
Note: The number of geometries may not
correspond to the number of cycles due to rejected geometry changes.
In the following table subtotal= G_P + SS G_CDS.
Atom Chem. CM2 G_P Area Sigma k SS G_CDS Subtotal M
number symbol chg. (kcal) (Ang**2) cal/(Ang**2) (kcal) (kcal) value
1 C -0.14 -1.71 9.07 37.16 0.34 -1.38 16
2 C -0.10 -1.26 2.01 -91.77 -0.18 -1.45 16
3 H 0.09 1.16 7.12 -51.93 -0.37 0.79 16
4 C 0.53 8.56 6.92 -10.99 -0.08 8.48 16
5 O -0.55 -10.30 16.29 5.55 0.09 -10.21 16
6 N -0.61 -9.77 3.33 -170.52 -0.57 -10.33 16
7 C 0.12 1.85 6.20 -2.53 -0.02 1.84 16
8 C -0.03 -0.47 3.18 -89.62 -0.28 -0.76 16
9 H 0.10 1.84 8.00 -51.93 -0.42 1.43 16
10 C 0.04 0.70 6.19 -82.85 -0.51 0.18 16
11 C -0.07 -0.93 10.68 -16.56 -0.18 -1.11 16
12 N -0.47 -5.75 3.15 -125.99 -0.40 -6.15 16
13 C 0.08 0.54 11.18 59.85 0.67 1.21 16
14 C 0.18 2.55 12.56 -90.38 -1.14 1.42 16
15 N -0.47 -8.80 11.25 -10.87 -0.12 -8.93 16
16 C 0.19 3.62 3.05 -66.10 -0.20 3.42 16
17 H 0.06 1.25 8.02 -51.93 -0.42 0.83 16
18 N -0.70 -17.12 5.18 -51.65 -0.27 -17.39 16
19 C -0.24 -6.55 9.98 -15.06 -0.15 -6.70 16
20 H 0.07 1.82 7.83 -52.49 -0.41 1.40 16
21 C -0.01 -0.18 5.50 -17.43 -0.10 -0.27 16
22 N -0.92 -27.56 13.06 55.49 0.72 -26.83 16
23 Si 0.90 22.01 29.55 -169.99 -5.02 16.99 16
24 H -0.28 -6.61 7.11 56.52 0.40 -6.21 16
25 H -0.29 -7.13 7.11 56.52 0.40 -6.72 16
26 H -0.29 -7.23 7.11 56.52 0.40 -6.83 16
27 C 0.07 1.22 5.84 -3.06 -0.02 1.20 16
28 C -0.09 -1.07 2.76 -89.81 -0.25 -1.32 16
29 C -0.14 -1.52 7.62 -25.83 -0.20 -1.71 16
30 C -0.14 -1.31 8.09 -25.83 -0.21 -1.52 16
31 C -0.14 -1.30 6.66 -25.92 -0.17 -1.47 16
32 H 0.07 0.71 7.34 -51.93 -0.38 0.33 16
33 H 0.06 0.74 8.14 -51.93 -0.42 0.32 16
34 H 0.07 0.88 7.59 -51.93 -0.39 0.49 16
35 H 0.07 1.16 8.09 -51.93 -0.42 0.74 16
36 H 0.06 0.90 6.35 -51.93 -0.33 0.57 16
37 H 0.17 1.95 6.87 -52.49 -0.36 1.59 16
38 H 0.07 0.41 8.14 -51.93 -0.42 -0.02 16
39 H 0.08 0.49 8.14 -51.93 -0.42 0.07 16
40 H 0.08 0.46 8.14 -51.93 -0.42 0.04 16
41 H 0.20 2.11 8.06 -52.48 -0.42 1.69 16
42 H 0.39 10.39 7.71 -40.82 -0.31 10.07 16
43 H 0.25 7.38 8.31 -40.82 -0.34 7.04 16
44 H 0.08 1.41 7.47 -51.93 -0.39 1.02 16
45 H 0.06 0.96 8.01 -51.93 -0.42 0.54 16
46 H 0.10 1.30 7.67 -51.93 -0.40 0.90 16
47 H 0.08 0.86 7.99 -51.93 -0.41 0.45 16
48 H 0.07 0.76 8.14 -51.93 -0.42 0.34 16
49 H 0.08 0.67 8.14 -51.93 -0.42 0.25 16
50 H 0.07 0.55 8.14 -51.93 -0.42 0.13 16
51 H 0.07 0.59 8.14 -51.93 -0.42 0.17 16
52 H 0.08 0.77 8.10 -51.93 -0.42 0.35 16
LS Contribution 416.27 15.07 6.27 6.27
Total: -1.00 -34.02 416.27 -10.75 -44.78
By element:
Atomic # 1 Polarization: 20.54 SS G_CDS: -8.79 Total: 11.75 kcal
Atomic # 6 Polarization: 2.72 SS G_CDS: -2.67 Total: 0.05 kcal
Atomic # 7 Polarization: -69.00 SS G_CDS: -0.63 Total: -69.63 kcal
Atomic # 8 Polarization: -10.30 SS G_CDS: 0.09 Total: -10.21 kcal
Atomic # 14 Polarization: 22.01 SS G_CDS: -5.02 Total: 16.99 kcal
Total LS contribution 6.27 Total: 6.27 kcal
Total: -34.02 -10.75 -44.78 kcal
The number of atoms in the molecule is 52
The average number of expansion shells was 16.00
The maximum number of expansion shells was 16
The minimum number of expansion shells was 16
**** NOTA BENE ****
This is the net solvation energy for this exact molecular structure
(nuclear and electronic)! The standard-state solvation energy should be
obtained as the difference between the heat of formation plus delta-G solvation
for the relaxed solvated system and that for the relaxed gas-phase system.
ZINC001070505771.mol2 52
1SCF run
This is a breakdown of the solvation energy calculated
without geometric relaxation in solution:
(1) E-EN(sol) electronic-nuclear energy of solute 106.621 kcal
(2) G-P(sol) polarization free energy of solvation -34.024 kcal
(3) G-ENP(sol) elect.-nuc.-pol. free energy of system 72.597 kcal
(4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.752 kcal
(5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -44.776 kcal
(6) G-S(sol) free energy of system = (1) + (5) 61.845 kcal
1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED
Total computer time = 1.21 seconds