Wall clock time and date at job start Tue Mar 31 2020 20:06:27 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.52996 * 1 3 3 H 1.09003 * 109.47073 * 2 1 4 4 C 1.52996 * 109.47246 * 239.99252 * 2 1 3 5 5 C 1.52994 * 109.47633 * 300.00471 * 4 2 1 6 6 C 1.53004 * 109.47061 * 300.00041 * 5 4 2 7 7 C 1.53001 * 109.46950 * 59.99865 * 6 5 4 8 8 H 1.08991 * 109.47541 * 60.00111 * 7 6 5 9 9 C 1.53001 * 109.46769 * 180.02562 * 7 6 5 10 10 C 1.52999 * 109.47066 * 299.99809 * 7 6 5 11 11 N 1.46501 * 109.47195 * 180.02562 * 10 7 6 12 12 C 1.34776 * 120.00213 * 155.00424 * 11 10 7 13 13 O 1.21348 * 119.99842 * 359.97438 * 12 11 10 14 14 C 1.49445 * 119.99936 * 180.02562 * 12 11 10 15 15 O 1.21349 * 120.00214 * 353.50159 * 14 12 11 16 16 N 1.34779 * 120.00153 * 173.50346 * 14 12 11 17 17 C 1.46927 * 120.62861 * 351.29699 * 16 14 12 18 18 C 1.53194 * 108.77358 * 233.58601 * 17 16 14 19 19 C 1.53041 * 109.31129 * 305.36269 * 18 17 16 20 20 C 1.53038 * 109.53879 * 61.36809 * 19 18 17 21 21 H 1.09003 * 109.49937 * 58.61143 * 20 19 18 22 22 C 1.50693 * 109.50018 * 178.68314 * 20 19 18 23 23 H 1.08001 * 119.99908 * 6.73620 * 22 20 19 24 24 C 1.37884 * 120.00340 * 186.74016 * 22 20 19 25 25 N 1.31514 * 120.00017 * 185.95516 * 24 22 20 26 26 Si 1.86797 * 119.99551 * 5.95762 * 24 22 20 27 27 H 1.48492 * 109.47361 * 269.99669 * 26 24 22 28 28 H 1.48504 * 109.46905 * 149.99867 * 26 24 22 29 29 H 1.48503 * 109.47412 * 30.00180 * 26 24 22 30 30 C 1.46923 * 120.63294 * 171.27003 * 16 14 12 31 31 H 1.08999 * 109.47103 * 300.00433 * 1 2 3 32 32 H 1.09005 * 109.46829 * 60.00038 * 1 2 3 33 33 H 1.08994 * 109.47330 * 180.02562 * 1 2 3 34 34 H 1.08999 * 109.47417 * 180.02562 * 4 2 1 35 35 H 1.09005 * 109.46626 * 60.00586 * 4 2 1 36 36 H 1.08999 * 109.47447 * 179.97438 * 5 4 2 37 37 H 1.09004 * 109.47008 * 59.99836 * 5 4 2 38 38 H 1.08999 * 109.46624 * 299.99590 * 6 5 4 39 39 H 1.08994 * 109.47092 * 179.97438 * 6 5 4 40 40 H 1.08999 * 109.47078 * 179.97438 * 9 7 6 41 41 H 1.09003 * 109.46930 * 300.00020 * 9 7 6 42 42 H 1.09004 * 109.47408 * 59.99796 * 9 7 6 43 43 H 1.09004 * 109.47271 * 300.00302 * 10 7 6 44 44 H 0.97006 * 119.99659 * 334.99617 * 11 10 7 45 45 H 1.08999 * 109.58916 * 113.79637 * 17 16 14 46 46 H 1.09004 * 109.70538 * 353.44067 * 17 16 14 47 47 H 1.08998 * 109.49744 * 185.41288 * 18 17 16 48 48 H 1.08999 * 109.52352 * 65.33476 * 18 17 16 49 49 H 1.09002 * 109.46067 * 301.34513 * 19 18 17 50 50 H 1.09000 * 109.45571 * 181.38639 * 19 18 17 51 51 H 0.97003 * 119.99822 * 354.97653 * 25 24 22 52 52 H 1.09000 * 109.58623 * 6.63017 * 30 16 14 53 53 H 1.08999 * 109.58751 * 246.19965 * 30 16 14 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.8933 1.0277 0.0000 4 6 2.0400 -0.7214 -1.2491 5 6 1.5300 -2.1638 -1.2491 6 6 2.0399 -2.8850 0.0002 7 6 1.5298 -2.1637 1.2494 8 1 0.4399 -2.1636 1.2494 9 6 2.0402 -2.8845 2.4987 10 6 2.0399 -0.7212 1.2493 11 7 1.5520 -0.0307 2.4457 12 6 2.2289 1.0261 2.9371 13 8 3.2428 1.4053 2.3887 14 6 1.7314 1.7302 4.1578 15 8 0.6588 1.4277 4.6380 16 7 2.4741 2.7012 4.7255 17 6 3.8506 2.9501 4.2761 18 6 4.7819 2.9043 5.4917 19 6 4.2793 3.8871 6.5516 20 6 2.8752 3.4790 7.0033 21 1 2.9020 2.4654 7.4034 22 6 2.3941 4.4284 8.0701 23 1 2.9772 5.3032 8.3173 24 6 1.2096 4.1801 8.7308 25 7 0.7305 5.0657 9.5769 26 14 0.2856 2.5859 8.4242 27 1 -0.6681 2.7745 7.3018 28 1 -0.4612 2.2032 9.6494 29 1 1.2516 1.5125 8.0779 30 6 1.9223 3.5313 5.8049 31 1 -0.3633 0.5139 0.8899 32 1 -0.3633 0.5139 -0.8900 33 1 -0.3634 -1.0276 0.0005 34 1 3.1300 -0.7218 -1.2489 35 1 1.6767 -0.2075 -2.1391 36 1 1.8930 -2.6776 -2.1392 37 1 0.4399 -2.1638 -1.2492 38 1 3.1299 -2.8851 0.0002 39 1 1.6762 -3.9125 0.0005 40 1 1.6772 -2.3704 3.3886 41 1 1.6772 -3.9124 2.4989 42 1 3.1303 -2.8842 2.4986 43 1 3.1299 -0.7213 1.2494 44 1 0.7415 -0.3340 2.8842 45 1 3.9100 3.9322 3.8071 46 1 4.1474 2.1843 3.5594 47 1 5.7910 3.1819 5.1874 48 1 4.7912 1.8960 5.9056 49 1 4.2473 4.8916 6.1294 50 1 4.9541 3.8745 7.4075 51 1 1.2451 5.8595 9.7913 52 1 0.9452 3.1473 6.0980 53 1 1.8245 4.5608 5.4605 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 ZINC001330851486.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 20:06:27 Heat of formation + Delta-G solvation = -88.363480 kcal Electronic energy + Delta-G solvation = -31240.268528 eV Core-core repulsion = 27303.930593 eV Total energy + Delta-G solvation = -3936.337935 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 1.13 seconds Orbital eigenvalues (eV) -41.03915 -39.79536 -38.14483 -37.10018 -34.60823 -33.53283 -33.27858 -31.92495 -30.73683 -30.39468 -27.23225 -26.79440 -25.69653 -24.87425 -23.37691 -21.99279 -21.79142 -21.43442 -20.00962 -19.23778 -18.81513 -17.67741 -17.21891 -16.89019 -16.14272 -15.79924 -15.32275 -15.10299 -15.02475 -14.32220 -14.13037 -14.06045 -13.88432 -13.65248 -13.46127 -13.42042 -13.20091 -12.95528 -12.85402 -12.62196 -12.36403 -12.32047 -11.96298 -11.91674 -11.62103 -11.59855 -11.39891 -11.19034 -11.01244 -10.95006 -10.76439 -10.68755 -10.45153 -10.36087 -10.33714 -10.16062 -9.64990 -9.51163 -9.27958 -8.94372 -8.70418 -8.43048 -8.13871 -6.45046 -4.14059 1.61226 2.52552 3.40539 3.49746 3.52557 3.65008 3.73579 4.20454 4.33059 4.43985 4.48953 4.53623 4.70864 4.76164 4.95302 4.98999 5.05467 5.12103 5.23047 5.26819 5.36268 5.41376 5.45526 5.48662 5.54692 5.55686 5.63468 5.63856 5.67901 5.75731 5.80904 5.86566 5.90180 5.95660 6.00385 6.04440 6.10590 6.13703 6.17234 6.21732 6.38211 6.48085 6.54848 6.62206 6.75845 6.80873 7.11004 7.17361 7.41392 7.69004 7.87546 8.04826 8.28370 8.34279 8.85970 9.57036 10.93625 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.018762 B = 0.003100 C = 0.002917 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1491.985941 B = 9029.980090 C = 9597.883080 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.152 4.152 2 C -0.103 4.103 3 H 0.077 0.923 4 C -0.111 4.111 5 C -0.124 4.124 6 C -0.109 4.109 7 C -0.101 4.101 8 H 0.081 0.919 9 C -0.151 4.151 10 C 0.154 3.846 11 N -0.696 5.696 12 C 0.513 3.487 13 O -0.515 6.515 14 C 0.494 3.506 15 O -0.531 6.531 16 N -0.571 5.571 17 C 0.073 3.927 18 C -0.136 4.136 19 C -0.102 4.102 20 C 0.005 3.995 21 H 0.033 0.967 22 C -0.590 4.590 23 H 0.059 0.941 24 C 0.030 3.970 25 N -0.909 5.909 26 Si 1.047 2.953 27 H -0.276 1.276 28 H -0.289 1.289 29 H -0.278 1.278 30 C 0.102 3.898 31 H 0.056 0.944 32 H 0.054 0.946 33 H 0.061 0.939 34 H 0.064 0.936 35 H 0.065 0.935 36 H 0.062 0.938 37 H 0.066 0.934 38 H 0.065 0.935 39 H 0.067 0.933 40 H 0.059 0.941 41 H 0.058 0.942 42 H 0.056 0.944 43 H 0.092 0.908 44 H 0.411 0.589 45 H 0.063 0.937 46 H 0.131 0.869 47 H 0.065 0.935 48 H 0.058 0.942 49 H 0.052 0.948 50 H 0.055 0.945 51 H 0.268 0.732 52 H 0.082 0.918 53 H 0.069 0.931 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 3.445 -15.777 -22.338 27.564 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.209 4.209 2 C -0.123 4.123 3 H 0.095 0.905 4 C -0.149 4.149 5 C -0.161 4.161 6 C -0.147 4.147 7 C -0.120 4.120 8 H 0.099 0.901 9 C -0.208 4.208 10 C 0.050 3.950 11 N -0.348 5.348 12 C 0.295 3.705 13 O -0.392 6.392 14 C 0.276 3.724 15 O -0.408 6.408 16 N -0.303 5.303 17 C -0.049 4.049 18 C -0.174 4.174 19 C -0.140 4.140 20 C -0.014 4.014 21 H 0.051 0.949 22 C -0.629 4.629 23 H 0.078 0.922 24 C -0.175 4.175 25 N -0.543 5.543 26 Si 0.875 3.125 27 H -0.203 1.203 28 H -0.215 1.215 29 H -0.204 1.204 30 C -0.019 4.019 31 H 0.075 0.925 32 H 0.073 0.927 33 H 0.080 0.920 34 H 0.083 0.917 35 H 0.083 0.917 36 H 0.081 0.919 37 H 0.085 0.915 38 H 0.084 0.916 39 H 0.085 0.915 40 H 0.078 0.922 41 H 0.077 0.923 42 H 0.075 0.925 43 H 0.110 0.890 44 H 0.248 0.752 45 H 0.082 0.918 46 H 0.147 0.853 47 H 0.083 0.917 48 H 0.077 0.923 49 H 0.071 0.929 50 H 0.074 0.926 51 H 0.078 0.922 52 H 0.100 0.900 53 H 0.087 0.913 Dipole moment (debyes) X Y Z Total from point charges 2.214 -16.452 -22.068 27.615 hybrid contribution 1.038 0.975 -0.087 1.427 sum 3.252 -15.476 -22.155 27.220 Atomic orbital electron populations 1.21935 0.94942 1.01989 1.02018 1.21491 0.96458 1.00175 0.94175 0.90494 1.21525 1.00808 0.95424 0.97171 1.21617 1.01211 0.95627 0.97688 1.21543 1.01074 0.98827 0.93287 1.21448 1.01316 0.93307 0.95909 0.90107 1.21865 1.01505 1.00224 0.97227 1.21126 0.99892 0.90601 0.83407 1.45891 1.29207 1.30648 1.29073 1.20490 0.83734 0.80903 0.85383 1.90474 1.23919 1.64583 1.60234 1.19994 0.85650 0.81921 0.84880 1.90747 1.22584 1.64261 1.63256 1.48404 1.14221 1.32242 1.35443 1.22743 0.80059 1.04275 0.97821 1.22005 0.98098 1.00705 0.96621 1.21483 0.98142 0.97545 0.96835 1.19485 0.91153 0.98290 0.92521 0.94941 1.22003 1.09996 1.10508 1.20371 0.92229 1.25822 0.95533 1.02429 0.93735 1.70819 1.38474 1.10676 1.34358 0.83933 0.76619 0.75716 0.76186 1.20257 1.21530 1.20419 1.22758 0.99662 0.92945 0.86491 0.92461 0.92677 0.92044 0.91717 0.91665 0.91904 0.91506 0.91624 0.91453 0.92197 0.92323 0.92476 0.89045 0.75202 0.91826 0.85252 0.91673 0.92293 0.92890 0.92577 0.92213 0.89962 0.91270 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 13. 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.15 -1.25 8.34 37.16 0.31 -0.94 16 2 C -0.10 -0.91 2.88 -90.62 -0.26 -1.17 16 3 H 0.08 0.84 8.14 -51.93 -0.42 0.42 16 4 C -0.11 -0.78 5.45 -26.74 -0.15 -0.93 16 5 C -0.12 -0.73 5.65 -26.73 -0.15 -0.88 16 6 C -0.11 -0.66 5.20 -26.73 -0.14 -0.80 16 7 C -0.10 -0.73 2.02 -90.62 -0.18 -0.92 16 8 H 0.08 0.57 6.43 -51.93 -0.33 0.23 16 9 C -0.15 -1.17 9.03 37.16 0.34 -0.83 16 10 C 0.15 1.44 1.59 -67.93 -0.11 1.33 16 11 N -0.70 -8.09 3.61 -54.23 -0.20 -8.29 16 12 C 0.51 7.68 7.71 -11.59 -0.09 7.59 16 13 O -0.51 -8.03 12.84 5.50 0.07 -7.96 16 14 C 0.49 9.17 7.73 -11.59 -0.09 9.08 16 15 O -0.53 -11.07 15.93 5.50 0.09 -10.98 16 16 N -0.57 -11.08 2.97 -173.27 -0.51 -11.60 16 17 C 0.07 1.20 5.98 -3.71 -0.02 1.18 16 18 C -0.14 -2.25 6.08 -26.61 -0.16 -2.41 16 19 C -0.10 -2.04 5.11 -26.69 -0.14 -2.18 16 20 C 0.01 0.13 1.82 -91.65 -0.17 -0.04 16 21 H 0.03 0.82 6.53 -51.93 -0.34 0.48 16 22 C -0.59 -16.49 8.90 -34.44 -0.31 -16.80 16 23 H 0.06 1.76 7.98 -52.49 -0.42 1.35 16 24 C 0.03 0.87 5.43 -17.82 -0.10 0.78 16 25 N -0.91 -28.14 13.70 55.43 0.76 -27.38 16 26 Si 1.05 26.89 25.42 -169.99 -4.32 22.57 16 27 H -0.28 -7.00 6.82 56.52 0.39 -6.62 16 28 H -0.29 -7.35 7.11 56.52 0.40 -6.94 16 29 H -0.28 -6.89 6.33 56.52 0.36 -6.54 16 30 C 0.10 2.30 5.33 -3.71 -0.02 2.28 16 31 H 0.06 0.55 7.61 -51.93 -0.39 0.16 16 32 H 0.05 0.43 8.14 -51.93 -0.42 0.01 16 33 H 0.06 0.44 4.81 -51.93 -0.25 0.19 16 34 H 0.06 0.47 8.14 -51.93 -0.42 0.04 16 35 H 0.06 0.43 8.14 -51.93 -0.42 0.01 16 36 H 0.06 0.33 8.14 -51.93 -0.42 -0.10 16 37 H 0.07 0.39 6.43 -51.93 -0.33 0.05 16 38 H 0.07 0.41 8.14 -51.93 -0.42 -0.01 16 39 H 0.07 0.36 8.14 -51.93 -0.42 -0.06 16 40 H 0.06 0.53 7.44 -51.93 -0.39 0.15 16 41 H 0.06 0.39 8.14 -51.93 -0.42 -0.03 16 42 H 0.06 0.46 8.14 -51.93 -0.42 0.03 16 43 H 0.09 0.97 7.58 -51.93 -0.39 0.57 16 44 H 0.41 4.59 7.88 -40.82 -0.32 4.27 16 45 H 0.06 1.01 8.14 -51.93 -0.42 0.58 16 46 H 0.13 2.06 4.32 -54.78 -0.24 1.82 16 47 H 0.06 0.91 8.14 -51.93 -0.42 0.48 16 48 H 0.06 0.99 8.14 -51.93 -0.42 0.57 16 49 H 0.05 1.05 8.14 -51.93 -0.42 0.62 16 50 H 0.06 1.12 8.14 -51.93 -0.42 0.69 16 51 H 0.27 8.18 8.72 -40.82 -0.36 7.83 16 52 H 0.08 2.03 3.55 -51.93 -0.18 1.84 16 53 H 0.07 1.55 8.14 -51.93 -0.42 1.13 16 LS Contribution 390.38 15.07 5.88 5.88 Total: -1.00 -31.39 390.38 -8.81 -40.20 By element: Atomic # 1 Polarization: 12.36 SS G_CDS: -9.15 Total: 3.21 kcal Atomic # 6 Polarization: -4.23 SS G_CDS: -1.43 Total: -5.66 kcal Atomic # 7 Polarization: -47.31 SS G_CDS: 0.05 Total: -47.27 kcal Atomic # 8 Polarization: -19.10 SS G_CDS: 0.16 Total: -18.94 kcal Atomic # 14 Polarization: 26.89 SS G_CDS: -4.32 Total: 22.57 kcal Total LS contribution 5.88 Total: 5.88 kcal Total: -31.39 -8.81 -40.20 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. ZINC001330851486.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 -48.162 kcal (2) G-P(sol) polarization free energy of solvation -31.392 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -79.554 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.809 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -40.201 kcal (6) G-S(sol) free energy of system = (1) + (5) -88.363 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.13 seconds