Wall clock time and date at job start Tue Mar 31 2020 19:25:36 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.52999 * 1 3 3 C 1.52999 * 109.47529 * 2 1 4 4 C 1.53000 * 109.47388 * 120.00583 * 2 1 3 5 5 C 1.52999 * 109.47530 * 174.99635 * 4 2 1 6 6 H 1.09000 * 109.12470 * 55.52156 * 5 4 2 7 7 C 1.52919 * 109.30221 * 174.99662 * 5 4 2 8 8 C 1.53339 * 108.67500 * 171.06647 * 7 5 4 9 9 H 1.08997 * 109.67567 * 172.69909 * 8 7 5 10 10 C 1.50706 * 109.67391 * 52.20656 * 8 7 5 11 11 O 1.21291 * 119.99537 * 261.57583 * 10 8 7 12 12 N 1.34767 * 120.00272 * 81.58125 * 10 8 7 13 13 C 1.37103 * 119.99958 * 356.28399 * 12 10 8 14 14 H 1.08006 * 119.99904 * 235.75582 * 13 12 10 15 15 C 1.38947 * 120.00229 * 55.75146 * 13 12 10 16 16 N 1.31415 * 120.00336 * 202.89695 * 15 13 12 17 17 Si 1.86802 * 119.99613 * 22.90207 * 15 13 12 18 18 H 1.48508 * 109.47185 * 180.34884 * 17 15 13 19 19 H 1.48496 * 109.47429 * 60.34922 * 17 15 13 20 20 H 1.48504 * 109.47198 * 300.34761 * 17 15 13 21 21 C 1.46500 * 120.00371 * 176.28798 * 12 10 8 22 22 C 1.54268 * 110.03106 * 237.90820 * 21 12 10 23 23 C 1.54968 * 104.15583 * 217.13237 * 22 21 12 24 24 C 1.54886 * 102.77692 * 37.88798 * 23 22 21 25 25 C 1.53869 * 110.03753 * 120.68784 * 21 12 10 26 26 C 1.53071 * 108.42721 * 292.47006 * 8 7 5 27 27 C 1.50176 * 109.88350 * 48.91695 * 26 8 7 28 28 O 1.21279 * 118.24547 * 162.38043 * 27 26 8 29 29 N 1.33353 * 123.50787 * 342.38165 * 27 26 8 30 30 H 1.09004 * 109.46685 * 299.99269 * 1 2 3 31 31 H 1.08993 * 109.47078 * 59.99807 * 1 2 3 32 32 H 1.09002 * 109.46904 * 179.97438 * 1 2 3 33 33 H 1.09000 * 109.46965 * 240.00045 * 2 1 3 34 34 H 1.08994 * 109.47328 * 179.97438 * 3 2 1 35 35 H 1.08993 * 109.47147 * 300.00735 * 3 2 1 36 36 H 1.09004 * 109.46540 * 60.00642 * 3 2 1 37 37 H 1.09005 * 109.47510 * 54.99327 * 4 2 1 38 38 H 1.09007 * 109.46901 * 294.99955 * 4 2 1 39 39 H 1.08999 * 109.60042 * 51.30906 * 7 5 4 40 40 H 1.08998 * 109.60564 * 290.82356 * 7 5 4 41 41 H 0.96998 * 119.99953 * 0.02562 * 16 15 13 42 42 H 1.09005 * 110.01492 * 359.29143 * 21 12 10 43 43 H 1.09001 * 110.48759 * 335.77733 * 22 21 12 44 44 H 1.09000 * 110.48734 * 98.42693 * 22 21 12 45 45 H 1.08993 * 110.72022 * 279.57049 * 23 22 21 46 46 H 1.09001 * 110.91903 * 156.29438 * 23 22 21 47 47 H 1.09005 * 110.48209 * 203.40145 * 24 23 22 48 48 H 1.09005 * 110.60079 * 80.81829 * 24 23 22 49 49 H 1.09000 * 110.03543 * 238.57562 * 25 21 12 50 50 H 1.09003 * 110.02951 * 359.94694 * 25 21 12 51 51 H 1.08997 * 109.40156 * 169.02891 * 26 8 7 52 52 H 1.08995 * 109.40205 * 288.77488 * 26 8 7 53 53 H 0.97005 * 117.94623 * 182.24993 * 29 27 26 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 6 2.0401 1.4425 0.0000 4 6 2.0401 -0.7214 1.2491 5 6 3.5649 -0.8312 1.1878 6 1 3.9899 0.1645 1.0620 7 6 4.0807 -1.4396 2.4926 8 6 5.5736 -1.7507 2.3328 9 1 5.9879 -2.0636 3.2911 10 6 6.2982 -0.5231 1.8436 11 8 6.5560 -0.4006 0.6647 12 7 6.6596 0.4385 2.7159 13 6 6.3018 0.3367 4.0355 14 1 5.7471 1.1341 4.5077 15 6 6.6527 -0.7924 4.7653 16 7 5.9699 -1.1188 5.8397 17 14 8.0951 -1.8464 4.2194 18 1 8.2562 -2.9881 5.1554 19 1 9.3342 -1.0279 4.2192 20 1 7.8416 -2.3612 2.8497 21 6 7.4348 1.5927 2.2542 22 6 8.7595 1.6876 3.0390 23 6 8.9789 3.2098 3.2295 24 6 7.5434 3.7339 3.4817 25 6 6.6621 2.8933 2.5351 26 6 5.7259 -2.8803 1.3111 27 6 4.9206 -2.5695 0.0822 28 8 5.1732 -3.1750 -0.9378 29 7 3.9427 -1.6634 0.0525 30 1 -0.3633 0.5138 0.8901 31 1 -0.3633 0.5138 -0.8899 32 1 -0.3633 -1.0277 -0.0005 33 1 1.8933 -0.5138 -0.8900 34 1 3.1300 1.4424 0.0005 35 1 1.6767 1.9563 -0.8899 36 1 1.6768 1.9562 0.8901 37 1 1.6054 -1.7199 1.2951 38 1 1.7517 -0.1585 2.1370 39 1 3.5364 -2.3592 2.7075 40 1 3.9402 -0.7298 3.3077 41 1 5.2207 -0.5717 6.1230 42 1 7.6387 1.5025 1.1872 43 1 8.6669 1.1919 4.0053 44 1 9.5763 1.2529 2.4628 45 1 9.3997 3.6546 2.3278 46 1 9.6186 3.4062 4.0900 47 1 7.4737 4.7919 3.2288 48 1 7.2523 3.5707 4.5194 49 1 6.4915 3.4347 1.6045 50 1 5.7098 2.6657 3.0141 51 1 6.7763 -2.9833 1.0390 52 1 5.3733 -3.8142 1.7487 53 1 3.4522 -1.5509 -0.7768 RHF calculation, no. of doubly occupied orbitals= 65 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 ZINC000357615708.mol2 53 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Tue Mar 31 2020 19:25:36 Heat of formation + Delta-G solvation = -93.549204 kcal Electronic energy + Delta-G solvation = -32753.752537 eV Core-core repulsion = 28817.189732 eV Total energy + Delta-G solvation = -3936.562805 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 0.96 seconds Orbital eigenvalues (eV) -40.56048 -38.89526 -37.27224 -36.94457 -34.80201 -34.30688 -32.40702 -31.78786 -29.99474 -28.16226 -27.48065 -27.20277 -26.35722 -24.38242 -23.62334 -22.01444 -20.61862 -20.18393 -19.75656 -19.32404 -18.63064 -17.11035 -16.74051 -16.07500 -15.74024 -15.41402 -15.21782 -14.86666 -14.60524 -14.33512 -14.09097 -13.83069 -13.61605 -13.43616 -13.36833 -13.06077 -12.77111 -12.53459 -12.39484 -12.21889 -12.07973 -11.84709 -11.58813 -11.39443 -11.36091 -11.04706 -10.97701 -10.80365 -10.68150 -10.65032 -10.63145 -10.57099 -10.38642 -10.27039 -10.20540 -10.08819 -9.91234 -9.64690 -9.17553 -8.88980 -8.57596 -8.48106 -7.76919 -6.81831 -4.35455 3.10034 3.25632 3.28981 3.38172 3.65829 3.82829 4.07355 4.33287 4.69677 4.79684 4.82784 4.87608 5.16550 5.17723 5.19072 5.25241 5.34049 5.36034 5.43714 5.55134 5.59012 5.65161 5.74169 5.79315 5.86998 5.89216 5.92569 5.98152 5.99748 6.08774 6.09878 6.12965 6.15207 6.17097 6.25010 6.31517 6.33331 6.35206 6.37814 6.38990 6.52430 6.53922 6.58359 6.64888 6.75270 6.87874 7.01528 7.28589 7.69947 7.89124 7.94402 8.48721 8.64441 8.78845 9.29946 9.71815 10.43700 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.010696 B = 0.005442 C = 0.004584 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2617.294911 B = 5144.088592 C = 6106.700696 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.145 4.145 2 C -0.098 4.098 3 C -0.151 4.151 4 C -0.125 4.125 5 C 0.140 3.860 6 H 0.078 0.922 7 C -0.106 4.106 8 C -0.105 4.105 9 H 0.113 0.887 10 C 0.498 3.502 11 O -0.586 6.586 12 N -0.460 5.460 13 C -0.420 4.420 14 H 0.082 0.918 15 C 0.038 3.962 16 N -0.864 5.864 17 Si 1.014 2.986 18 H -0.278 1.278 19 H -0.277 1.277 20 H -0.269 1.269 21 C 0.126 3.874 22 C -0.132 4.132 23 C -0.120 4.120 24 C -0.116 4.116 25 C -0.122 4.122 26 C -0.122 4.122 27 C 0.545 3.455 28 O -0.578 6.578 29 N -0.706 5.706 30 H 0.055 0.945 31 H 0.052 0.948 32 H 0.053 0.947 33 H 0.067 0.933 34 H 0.061 0.939 35 H 0.051 0.949 36 H 0.055 0.945 37 H 0.068 0.932 38 H 0.078 0.922 39 H 0.069 0.931 40 H 0.099 0.901 41 H 0.278 0.722 42 H 0.080 0.920 43 H 0.087 0.913 44 H 0.052 0.948 45 H 0.058 0.942 46 H 0.059 0.941 47 H 0.052 0.948 48 H 0.065 0.935 49 H 0.051 0.949 50 H 0.071 0.929 51 H 0.099 0.901 52 H 0.096 0.904 53 H 0.389 0.611 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -5.962 4.441 -3.265 8.120 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.202 4.202 2 C -0.117 4.117 3 C -0.208 4.208 4 C -0.163 4.163 5 C 0.035 3.965 6 H 0.096 0.904 7 C -0.144 4.144 8 C -0.126 4.126 9 H 0.131 0.869 10 C 0.288 3.712 11 O -0.470 6.470 12 N -0.181 5.181 13 C -0.547 4.547 14 H 0.099 0.901 15 C -0.173 4.173 16 N -0.495 5.495 17 Si 0.840 3.160 18 H -0.203 1.203 19 H -0.203 1.203 20 H -0.195 1.195 21 C 0.023 3.977 22 C -0.170 4.170 23 C -0.158 4.158 24 C -0.154 4.154 25 C -0.161 4.161 26 C -0.162 4.162 27 C 0.334 3.666 28 O -0.458 6.458 29 N -0.357 5.357 30 H 0.074 0.926 31 H 0.071 0.929 32 H 0.072 0.928 33 H 0.085 0.915 34 H 0.080 0.920 35 H 0.070 0.930 36 H 0.074 0.926 37 H 0.087 0.913 38 H 0.097 0.903 39 H 0.088 0.912 40 H 0.118 0.882 41 H 0.089 0.911 42 H 0.098 0.902 43 H 0.105 0.895 44 H 0.071 0.929 45 H 0.077 0.923 46 H 0.078 0.922 47 H 0.071 0.929 48 H 0.084 0.916 49 H 0.070 0.930 50 H 0.090 0.910 51 H 0.117 0.883 52 H 0.115 0.885 53 H 0.225 0.775 Dipole moment (debyes) X Y Z Total from point charges -4.726 2.844 -3.721 6.654 hybrid contribution -0.924 1.611 -0.305 1.882 sum -5.650 4.455 -4.026 8.245 Atomic orbital electron populations 1.21823 0.95224 1.01558 1.01620 1.21120 0.95637 0.96623 0.98328 1.21907 1.01804 0.95618 1.01473 1.22049 0.94698 1.01103 0.98454 1.21050 0.94052 0.94975 0.86461 0.90407 1.21686 0.93771 1.02168 0.96812 1.21509 0.96828 0.93664 1.00621 0.86889 1.20439 0.79891 0.85485 0.85431 1.90581 1.58829 1.80551 1.17027 1.43771 1.47104 1.18918 1.08333 1.20086 1.39579 1.08801 0.86247 0.90091 1.26384 0.99887 0.96687 0.94294 1.70859 1.17683 1.36461 1.24520 0.84622 0.76982 0.77402 0.77013 1.20314 1.20319 1.19467 1.21692 0.90644 0.87092 0.98319 1.22737 0.96515 0.95618 1.02147 1.22410 0.96711 0.96236 1.00430 1.22249 0.95194 0.97874 1.00066 1.22330 0.99733 0.95619 0.98373 1.21576 1.01812 0.99176 0.93616 1.19843 0.78855 0.77847 0.90082 1.90657 1.68336 1.51985 1.34855 1.45421 1.36139 1.39177 1.14962 0.92635 0.92901 0.92806 0.91515 0.92019 0.93013 0.92585 0.91347 0.90347 0.91242 0.88207 0.91122 0.90218 0.89483 0.92871 0.92336 0.92191 0.92875 0.91594 0.92998 0.91048 0.88281 0.88549 0.77499 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 15. 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.36 9.19 37.16 0.34 -1.02 16 2 C -0.10 -1.09 2.84 -90.62 -0.26 -1.35 16 3 C -0.15 -1.74 8.95 37.16 0.33 -1.41 16 4 C -0.12 -1.53 4.33 -26.73 -0.12 -1.65 16 5 C 0.14 2.20 2.28 -68.38 -0.16 2.04 16 6 H 0.08 1.39 5.53 -51.93 -0.29 1.10 16 7 C -0.11 -1.83 4.82 -26.60 -0.13 -1.96 16 8 C -0.11 -2.11 1.29 -91.56 -0.12 -2.23 16 9 H 0.11 2.42 2.98 -51.93 -0.15 2.27 16 10 C 0.50 11.00 3.94 -10.99 -0.04 10.96 16 11 O -0.59 -14.22 14.14 5.55 0.08 -14.14 16 12 N -0.46 -9.90 1.79 -106.54 -0.19 -10.09 16 13 C -0.42 -9.43 6.39 -14.93 -0.10 -9.53 16 14 H 0.08 1.95 7.40 -52.48 -0.39 1.56 16 15 C 0.04 0.89 3.90 -17.47 -0.07 0.82 16 16 N -0.86 -21.25 13.69 55.50 0.76 -20.49 16 17 Si 1.01 22.30 24.93 -169.99 -4.24 18.07 16 18 H -0.28 -6.16 7.11 56.52 0.40 -5.76 16 19 H -0.28 -6.02 6.92 56.52 0.39 -5.62 16 20 H -0.27 -5.79 4.84 56.52 0.27 -5.52 16 21 C 0.13 2.50 3.29 -66.81 -0.22 2.28 16 22 C -0.13 -2.44 6.37 -25.03 -0.16 -2.60 16 23 C -0.12 -1.82 6.91 -24.70 -0.17 -1.99 16 24 C -0.12 -1.77 6.83 -25.08 -0.17 -1.94 16 25 C -0.12 -2.11 6.27 -25.62 -0.16 -2.27 16 26 C -0.12 -2.38 4.88 -28.10 -0.14 -2.51 16 27 C 0.55 11.03 6.40 -11.87 -0.08 10.96 16 28 O -0.58 -13.37 17.78 5.56 0.10 -13.27 16 29 N -0.71 -12.11 4.96 -54.86 -0.27 -12.39 16 30 H 0.05 0.50 8.14 -51.93 -0.42 0.08 16 31 H 0.05 0.46 8.14 -51.93 -0.42 0.04 16 32 H 0.05 0.48 8.14 -51.93 -0.42 0.06 16 33 H 0.07 0.77 6.46 -51.93 -0.34 0.43 16 34 H 0.06 0.83 6.46 -51.93 -0.34 0.50 16 35 H 0.05 0.54 8.14 -51.93 -0.42 0.12 16 36 H 0.06 0.62 8.14 -51.93 -0.42 0.19 16 37 H 0.07 0.80 8.14 -51.93 -0.42 0.38 16 38 H 0.08 0.92 8.14 -51.93 -0.42 0.50 16 39 H 0.07 1.08 8.11 -51.93 -0.42 0.66 16 40 H 0.10 1.82 7.37 -51.93 -0.38 1.44 16 41 H 0.28 6.72 8.73 -40.82 -0.36 6.36 16 42 H 0.08 1.68 6.95 -51.93 -0.36 1.32 16 43 H 0.09 1.80 5.93 -51.93 -0.31 1.49 16 44 H 0.05 0.96 8.14 -51.93 -0.42 0.53 16 45 H 0.06 0.81 8.14 -51.93 -0.42 0.38 16 46 H 0.06 0.84 8.14 -51.93 -0.42 0.42 16 47 H 0.05 0.70 8.14 -51.93 -0.42 0.28 16 48 H 0.07 1.06 8.14 -51.93 -0.42 0.63 16 49 H 0.05 0.83 8.14 -51.93 -0.42 0.41 16 50 H 0.07 1.32 6.93 -51.93 -0.36 0.96 16 51 H 0.10 2.07 7.53 -51.93 -0.39 1.67 16 52 H 0.10 1.61 8.14 -51.93 -0.42 1.18 16 53 H 0.39 6.17 6.90 -40.82 -0.28 5.89 16 LS Contribution 386.31 15.07 5.82 5.82 Total: -1.00 -27.39 386.31 -8.56 -35.95 By element: Atomic # 1 Polarization: 23.17 SS G_CDS: -9.22 Total: 13.95 kcal Atomic # 6 Polarization: -2.00 SS G_CDS: -1.40 Total: -3.41 kcal Atomic # 7 Polarization: -43.27 SS G_CDS: 0.30 Total: -42.97 kcal Atomic # 8 Polarization: -27.59 SS G_CDS: 0.18 Total: -27.41 kcal Atomic # 14 Polarization: 22.30 SS G_CDS: -4.24 Total: 18.07 kcal Total LS contribution 5.82 Total: 5.82 kcal Total: -27.39 -8.56 -35.95 kcal The number of atoms in the molecule is 53 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. ZINC000357615708.mol2 53 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 -57.601 kcal (2) G-P(sol) polarization free energy of solvation -27.388 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -84.989 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.561 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -35.948 kcal (6) G-S(sol) free energy of system = (1) + (5) -93.549 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.97 seconds