Wall clock time and date at job start Wed Apr 1 2020 00:24:55 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.08995 * 109.51372 * 2 1 4 4 C 1.53295 * 109.51655 * 120.10801 * 2 1 3 5 5 N 1.47149 * 108.38941 * 172.31209 * 4 2 1 6 6 C 1.34782 * 120.98569 * 129.39961 * 5 4 2 7 7 O 1.21585 * 119.99705 * 185.21244 * 6 5 4 8 8 N 1.34767 * 119.99977 * 5.21035 * 6 5 4 9 9 C 1.46505 * 120.00055 * 175.02718 * 8 6 5 10 10 C 1.52995 * 109.47064 * 179.97438 * 9 8 6 11 11 C 1.52995 * 109.47514 * 299.99999 * 10 9 8 12 12 C 1.53002 * 109.47260 * 60.00266 * 10 9 8 13 13 C 1.50705 * 109.46939 * 179.97438 * 10 9 8 14 14 H 1.07994 * 119.99875 * 0.02562 * 13 10 9 15 15 C 1.37882 * 119.99847 * 179.97438 * 13 10 9 16 16 N 1.31514 * 120.00170 * 0.02562 * 15 13 10 17 17 Si 1.86799 * 119.99812 * 180.02562 * 15 13 10 18 18 H 1.48505 * 109.46988 * 90.00218 * 17 15 13 19 19 H 1.48500 * 109.46987 * 209.99979 * 17 15 13 20 20 H 1.48502 * 109.47086 * 330.00261 * 17 15 13 21 21 C 1.47167 * 118.02565 * 309.66606 * 5 4 2 22 22 C 1.52721 * 108.65506 * 50.02537 * 21 5 4 23 23 C 1.53065 * 109.32744 * 68.36541 * 22 21 5 24 24 C 1.53116 * 109.17076 * 176.67618 * 23 22 21 25 25 O 1.42905 * 109.40315 * 302.54062 * 24 23 22 26 26 C 1.42900 * 114.13625 * 61.18033 * 25 24 23 27 27 C 1.53077 * 109.54830 * 188.01853 * 22 21 5 28 28 O 1.42834 * 109.53803 * 308.24741 * 22 21 5 29 29 H 1.09003 * 109.46637 * 299.89134 * 1 2 3 30 30 H 1.08990 * 109.47307 * 59.89378 * 1 2 3 31 31 H 1.08997 * 109.46950 * 179.89685 * 1 2 3 32 32 H 1.09006 * 109.67751 * 52.58810 * 4 2 1 33 33 H 1.08997 * 109.68306 * 292.03813 * 4 2 1 34 34 H 0.97002 * 119.99976 * 355.02690 * 8 6 5 35 35 H 1.08997 * 109.46959 * 60.00221 * 9 8 6 36 36 H 1.08999 * 109.46798 * 300.00115 * 9 8 6 37 37 H 1.09001 * 109.47330 * 60.00061 * 11 10 9 38 38 H 1.09010 * 109.47077 * 299.99882 * 11 10 9 39 39 H 1.08997 * 109.47741 * 180.02562 * 11 10 9 40 40 H 1.09001 * 109.46781 * 299.99447 * 12 10 9 41 41 H 1.08999 * 109.47291 * 59.99404 * 12 10 9 42 42 H 1.09001 * 109.47284 * 179.97438 * 12 10 9 43 43 H 0.96998 * 120.00148 * 179.97438 * 16 15 13 44 44 H 1.09002 * 109.67804 * 169.81882 * 21 5 4 45 45 H 1.08999 * 109.67809 * 290.35397 * 21 5 4 46 46 H 1.09006 * 109.52518 * 296.58309 * 23 22 21 47 47 H 1.09002 * 109.52765 * 56.76748 * 23 22 21 48 48 H 1.08993 * 109.48449 * 182.55336 * 24 23 22 49 49 H 1.08999 * 109.48005 * 62.51243 * 24 23 22 50 50 H 1.08993 * 109.48954 * 178.82035 * 26 25 24 51 51 H 1.09002 * 109.49566 * 58.78552 * 26 25 24 52 52 H 1.08999 * 109.52069 * 303.34571 * 27 22 21 53 53 H 1.08992 * 109.52472 * 63.53962 * 27 22 21 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.8941 1.0273 0.0000 4 6 2.0422 -0.7248 1.2499 5 7 3.5016 -0.8740 1.1348 6 6 4.3187 -0.4995 2.1391 7 8 5.5246 -0.5419 1.9896 8 7 3.7991 -0.0782 3.3091 9 6 4.6840 0.2248 4.4367 10 6 3.8469 0.6794 5.6340 11 6 3.0548 1.9326 5.2562 12 6 2.8778 -0.4358 6.0316 13 6 4.7571 0.9905 6.7942 14 1 5.8256 0.8758 6.6879 15 6 4.2256 1.4224 7.9908 16 7 2.9243 1.5625 8.1200 17 14 5.3538 1.8087 9.4286 18 1 5.7524 3.2382 9.3748 19 1 4.6404 1.5414 10.7033 20 1 6.5655 0.9537 9.3506 21 6 4.0427 -1.4386 -0.1118 22 6 3.4301 -0.6899 -1.2936 23 6 3.9572 0.7470 -1.3112 24 6 3.4062 1.4696 -2.5435 25 8 3.7827 0.7536 -3.7216 26 6 3.2751 -0.5808 -3.7816 27 6 3.8212 -1.3863 -2.5995 28 8 2.0076 -0.6770 -1.1653 29 1 -0.3633 0.5122 0.8910 30 1 -0.3633 0.5154 -0.8889 31 1 -0.3633 -1.0276 -0.0019 32 1 1.5770 -1.7081 1.3205 33 1 1.8019 -0.1400 2.1378 34 1 2.8393 0.0244 3.4046 35 1 5.3740 1.0196 4.1537 36 1 5.2486 -0.6680 4.7056 37 1 2.3965 1.7072 4.4172 38 1 3.7453 2.7272 4.9730 39 1 2.4588 2.2568 6.1093 40 1 2.2194 -0.6612 5.1926 41 1 3.4420 -1.3287 6.3007 42 1 2.2817 -0.1122 6.8849 43 1 2.5503 1.8660 8.9620 44 1 5.1264 -1.3224 -0.1273 45 1 3.7856 -2.4958 -0.1779 46 1 3.6319 1.2655 -0.4092 47 1 5.0464 0.7350 -1.3527 48 1 3.8145 2.4793 -2.5857 49 1 2.3192 1.5191 -2.4801 50 1 3.5886 -1.0475 -4.7153 51 1 2.1863 -0.5583 -3.7342 52 1 4.9073 -1.4452 -2.6697 53 1 3.3990 -2.3909 -2.6172 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 ZINC000590948383.mol2 53 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 00:24:55 Heat of formation + Delta-G solvation = -116.566160 kcal Electronic energy + Delta-G solvation = -32378.423042 eV Core-core repulsion = 28250.406095 eV Total energy + Delta-G solvation = -4128.016947 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.97 seconds Orbital eigenvalues (eV) -41.17514 -39.43809 -37.96528 -36.46519 -34.68895 -34.33007 -33.17439 -31.42902 -31.12318 -30.79040 -27.85896 -26.45925 -25.76564 -25.53480 -24.71504 -24.07348 -23.02206 -21.87149 -20.11844 -19.45327 -18.92953 -18.03524 -16.94820 -16.49968 -16.39506 -16.05781 -15.75541 -15.31802 -15.17862 -15.07622 -14.79507 -14.48548 -14.23057 -13.96159 -13.64029 -13.11319 -13.01176 -12.86335 -12.79712 -12.58376 -12.44872 -12.36461 -12.31371 -12.04445 -11.81866 -11.62569 -11.42798 -11.30010 -11.22033 -10.89637 -10.72800 -10.71110 -10.59463 -10.42299 -10.22640 -10.10135 -9.85787 -9.70132 -9.66729 -9.33731 -8.82490 -8.67478 -8.42269 -8.16651 -6.04699 -4.14575 2.74090 3.26848 3.29393 3.33909 3.38210 3.49527 3.66005 4.17484 4.32017 4.39263 4.45325 4.47683 4.62054 4.73601 4.79318 4.88211 5.01041 5.07382 5.16750 5.22293 5.26893 5.35023 5.37880 5.43667 5.44377 5.48873 5.54570 5.63114 5.66810 5.75785 5.82514 5.92703 5.94837 6.13676 6.32644 6.41503 6.46502 6.75116 6.77411 6.86325 6.95494 6.96783 6.99201 7.03697 7.16999 7.34317 7.52209 7.67185 7.97624 8.08790 8.30826 8.37507 8.84781 8.94099 9.51341 10.96833 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.022120 B = 0.002645 C = 0.002551 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1265.504862 B =10584.987069 C =10974.660648 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.143 4.143 2 C 0.071 3.929 3 H 0.067 0.933 4 C 0.083 3.917 5 N -0.637 5.637 6 C 0.703 3.297 7 O -0.594 6.594 8 N -0.732 5.732 9 C 0.149 3.851 10 C 0.057 3.943 11 C -0.129 4.129 12 C -0.129 4.129 13 C -0.505 4.505 14 H 0.058 0.942 15 C 0.066 3.934 16 N -0.886 5.886 17 Si 0.894 3.106 18 H -0.281 1.281 19 H -0.287 1.287 20 H -0.274 1.274 21 C 0.112 3.888 22 C 0.106 3.894 23 C -0.177 4.177 24 C 0.067 3.933 25 O -0.384 6.384 26 C 0.065 3.935 27 C -0.138 4.138 28 O -0.366 6.366 29 H 0.076 0.924 30 H 0.066 0.934 31 H 0.066 0.934 32 H 0.076 0.924 33 H 0.088 0.912 34 H 0.392 0.608 35 H 0.053 0.947 36 H 0.052 0.948 37 H 0.003 0.997 38 H 0.027 0.973 39 H 0.086 0.914 40 H 0.004 0.996 41 H 0.027 0.973 42 H 0.086 0.914 43 H 0.260 0.740 44 H 0.102 0.898 45 H 0.069 0.931 46 H 0.079 0.921 47 H 0.084 0.916 48 H 0.092 0.908 49 H 0.053 0.947 50 H 0.093 0.907 51 H 0.058 0.942 52 H 0.088 0.912 53 H 0.084 0.916 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -6.398 -4.199 -20.639 22.012 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.201 4.201 2 C 0.013 3.987 3 H 0.085 0.915 4 C -0.041 4.041 5 N -0.376 5.376 6 C 0.405 3.595 7 O -0.473 6.473 8 N -0.388 5.388 9 C 0.026 3.974 10 C 0.056 3.944 11 C -0.187 4.187 12 C -0.187 4.187 13 C -0.544 4.544 14 H 0.076 0.924 15 C -0.135 4.135 16 N -0.525 5.525 17 Si 0.722 3.278 18 H -0.207 1.207 19 H -0.213 1.213 20 H -0.199 1.199 21 C -0.012 4.012 22 C 0.066 3.934 23 C -0.215 4.215 24 C -0.010 4.010 25 O -0.303 6.303 26 C -0.011 4.011 27 C -0.176 4.176 28 O -0.285 6.285 29 H 0.095 0.905 30 H 0.085 0.915 31 H 0.085 0.915 32 H 0.094 0.906 33 H 0.107 0.893 34 H 0.230 0.770 35 H 0.071 0.929 36 H 0.070 0.930 37 H 0.022 0.978 38 H 0.046 0.954 39 H 0.104 0.896 40 H 0.023 0.977 41 H 0.046 0.954 42 H 0.104 0.896 43 H 0.070 0.930 44 H 0.120 0.880 45 H 0.087 0.913 46 H 0.098 0.902 47 H 0.103 0.897 48 H 0.111 0.889 49 H 0.072 0.928 50 H 0.111 0.889 51 H 0.077 0.923 52 H 0.107 0.893 53 H 0.102 0.898 Dipole moment (debyes) X Y Z Total from point charges -5.621 -4.187 -20.800 21.949 hybrid contribution 0.071 0.222 1.449 1.468 sum -5.551 -3.965 -19.351 20.518 Atomic orbital electron populations 1.21739 0.92074 1.03299 1.02975 1.22499 0.95374 0.95834 0.85011 0.91525 1.22335 0.81131 1.01687 0.98929 1.47245 1.05573 1.68254 1.16532 1.15939 0.84964 0.78255 0.80370 1.90875 1.12637 1.61901 1.81847 1.45669 1.11054 1.70552 1.11525 1.21615 0.92230 0.99213 0.84315 1.18357 0.93046 0.92345 0.90661 1.21818 0.99022 0.98520 0.99308 1.21822 0.98911 0.97717 1.00268 1.20983 0.92985 1.44708 0.95731 0.92389 1.24875 0.95058 0.95114 0.98438 1.69798 1.10283 1.46889 1.25542 0.86566 0.80820 0.78704 0.81756 1.20701 1.21262 1.19886 1.21925 1.00794 0.96631 0.81859 1.21688 0.81555 0.95037 0.95109 1.22717 1.02752 0.96382 0.99676 1.22495 0.99105 0.95264 0.84089 1.87911 1.75559 1.25764 1.41039 1.22501 0.98398 0.82495 0.97742 1.22142 1.02377 1.00066 0.93020 1.87925 1.20137 1.77201 1.43270 0.90497 0.91526 0.91486 0.90621 0.89347 0.77050 0.92855 0.92974 0.97814 0.95422 0.89562 0.97747 0.95412 0.89573 0.93000 0.87962 0.91258 0.90214 0.89737 0.88940 0.92830 0.88865 0.92350 0.89321 0.89755 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 12. 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.09 9.57 37.16 0.36 -0.73 16 2 C 0.07 0.70 3.02 -26.58 -0.08 0.62 16 3 H 0.07 0.69 6.21 -51.93 -0.32 0.37 16 4 C 0.08 0.93 5.86 -3.54 -0.02 0.91 16 5 N -0.64 -8.87 2.63 -178.35 -0.47 -9.34 16 6 C 0.70 12.69 8.31 -86.92 -0.72 11.97 16 7 O -0.59 -12.83 15.82 5.29 0.08 -12.75 16 8 N -0.73 -13.27 4.87 -65.23 -0.32 -13.59 16 9 C 0.15 3.20 4.51 -4.04 -0.02 3.18 16 10 C 0.06 1.30 0.50 -155.66 -0.08 1.22 16 11 C -0.13 -2.86 8.08 37.16 0.30 -2.56 16 12 C -0.13 -2.89 8.08 37.16 0.30 -2.59 16 13 C -0.51 -13.16 8.01 -34.44 -0.28 -13.44 16 14 H 0.06 1.53 7.74 -52.49 -0.41 1.12 16 15 C 0.07 1.73 5.15 -17.82 -0.09 1.64 16 16 N -0.89 -23.36 9.35 55.43 0.52 -22.84 16 17 Si 0.89 20.61 29.54 -169.99 -5.02 15.59 16 18 H -0.28 -6.48 7.11 56.52 0.40 -6.08 16 19 H -0.29 -6.59 7.11 56.52 0.40 -6.19 16 20 H -0.27 -6.34 7.11 56.52 0.40 -5.93 16 21 C 0.11 1.37 5.28 -3.82 -0.02 1.35 16 22 C 0.11 1.11 0.68 -90.69 -0.06 1.05 16 23 C -0.18 -1.96 4.45 -26.64 -0.12 -2.08 16 24 C 0.07 0.67 6.70 37.22 0.25 0.92 16 25 O -0.38 -4.18 10.76 -35.23 -0.38 -4.56 16 26 C 0.07 0.56 6.70 37.22 0.25 0.81 16 27 C -0.14 -1.18 5.14 -26.63 -0.14 -1.32 16 28 O -0.37 -3.83 8.86 -35.23 -0.31 -4.14 16 29 H 0.08 0.51 8.14 -51.93 -0.42 0.09 16 30 H 0.07 0.48 8.14 -51.93 -0.42 0.06 16 31 H 0.07 0.49 8.14 -51.93 -0.42 0.07 16 32 H 0.08 0.79 8.14 -51.93 -0.42 0.36 16 33 H 0.09 0.97 5.77 -51.93 -0.30 0.68 16 34 H 0.39 6.42 3.72 -40.82 -0.15 6.27 16 35 H 0.05 1.21 8.14 -51.93 -0.42 0.79 16 36 H 0.05 1.18 8.14 -51.93 -0.42 0.76 16 37 H 0.00 0.06 6.86 -51.93 -0.36 -0.30 16 38 H 0.03 0.61 8.14 -51.92 -0.42 0.19 16 39 H 0.09 2.12 6.07 -51.93 -0.32 1.80 16 40 H 0.00 0.07 6.86 -51.93 -0.36 -0.29 16 41 H 0.03 0.61 8.14 -51.93 -0.42 0.19 16 42 H 0.09 2.12 6.07 -51.93 -0.32 1.80 16 43 H 0.26 6.62 8.31 -40.82 -0.34 6.28 16 44 H 0.10 1.41 7.05 -51.93 -0.37 1.04 16 45 H 0.07 0.78 8.13 -51.93 -0.42 0.36 16 46 H 0.08 0.97 6.07 -51.93 -0.32 0.66 16 47 H 0.08 1.01 8.14 -51.93 -0.42 0.59 16 48 H 0.09 0.84 8.14 -51.93 -0.42 0.42 16 49 H 0.05 0.52 7.83 -51.93 -0.41 0.12 16 50 H 0.09 0.63 8.14 -51.93 -0.42 0.21 16 51 H 0.06 0.51 7.83 -51.93 -0.41 0.11 16 52 H 0.09 0.73 8.14 -51.93 -0.42 0.31 16 53 H 0.08 0.63 8.14 -51.93 -0.42 0.20 16 LS Contribution 393.56 15.07 5.93 5.93 Total: -1.00 -29.51 393.56 -9.21 -38.72 By element: Atomic # 1 Polarization: 15.10 SS G_CDS: -9.07 Total: 6.03 kcal Atomic # 6 Polarization: 1.12 SS G_CDS: -0.17 Total: 0.95 kcal Atomic # 7 Polarization: -45.50 SS G_CDS: -0.27 Total: -45.77 kcal Atomic # 8 Polarization: -20.85 SS G_CDS: -0.61 Total: -21.46 kcal Atomic # 14 Polarization: 20.61 SS G_CDS: -5.02 Total: 15.59 kcal Total LS contribution 5.93 Total: 5.93 kcal Total: -29.51 -9.21 -38.72 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. ZINC000590948383.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 -77.847 kcal (2) G-P(sol) polarization free energy of solvation -29.511 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -107.359 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.207 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -38.719 kcal (6) G-S(sol) free energy of system = (1) + (5) -116.566 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.97 seconds