Wall clock time and date at job start Tue Mar 31 2020 07:00:35 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 * SM5.42R - SM5.42R CALCULATIONS WILL BE PERFORMED * - CHARGE MODEL 2 WILL BE USED (CM2A) * - THE SOLVENT IS WATER ******************************************************************************* 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 O 1.42902 * 1 3 3 C 1.35989 * 117.00207 * 2 1 4 4 C 1.38570 * 120.07458 * 0.02562 * 3 2 1 5 5 C 1.38238 * 120.00282 * 179.97438 * 4 3 2 6 6 C 1.38258 * 120.14337 * 0.02562 * 5 4 3 7 7 C 1.50700 * 119.92400 * 179.97438 * 6 5 4 8 8 N 1.46903 * 109.47243 * 89.99728 * 7 6 5 9 9 C 1.46898 * 111.00266 * 179.97438 * 8 7 6 10 10 H 1.09000 * 109.47049 * 39.45060 * 9 8 7 11 11 C 1.53004 * 109.46941 * 279.45803 * 9 8 7 12 12 C 1.52994 * 109.47193 * 159.45517 * 9 8 7 13 13 N 1.46499 * 109.47398 * 302.09308 * 12 9 8 14 14 C 1.46761 * 124.35068 * 276.16783 * 13 12 9 15 15 C 1.54484 * 105.54943 * 162.91662 * 14 13 12 16 16 C 1.54858 * 101.66135 * 25.43011 * 15 14 13 17 17 C 1.34427 * 124.34594 * 96.13930 * 13 12 9 18 18 O 1.21273 * 124.94326 * 0.45240 * 17 13 12 19 19 C 1.38236 * 120.14712 * 0.22795 * 6 5 4 20 20 C 1.38569 * 120.00037 * 359.49444 * 19 6 5 21 21 O 1.35995 * 120.07197 * 180.25138 * 20 19 6 22 22 C 1.42899 * 117.00208 * 0.27568 * 21 20 19 23 23 C 1.53003 * 109.46895 * 179.97438 * 22 21 20 24 24 C 1.50699 * 109.46906 * 180.02562 * 23 22 21 25 25 H 1.08002 * 120.00156 * 0.02562 * 24 23 22 26 26 C 1.37873 * 119.99888 * 180.02562 * 24 23 22 27 27 N 1.31511 * 120.00486 * 0.02562 * 26 24 23 28 28 Si 1.86806 * 120.00128 * 180.02562 * 26 24 23 29 29 H 1.48506 * 109.46871 * 359.97438 * 28 26 24 30 30 H 1.48498 * 109.47528 * 119.99786 * 28 26 24 31 31 H 1.48501 * 109.47313 * 240.00088 * 28 26 24 32 32 H 1.09001 * 109.47303 * 299.99994 * 1 2 3 33 33 H 1.09002 * 109.46849 * 59.99844 * 1 2 3 34 34 H 1.08997 * 109.47666 * 180.02562 * 1 2 3 35 35 H 1.08001 * 119.99454 * 359.97438 * 4 3 2 36 36 H 1.07993 * 119.93186 * 179.97438 * 5 4 3 37 37 H 1.08999 * 109.47116 * 329.99750 * 7 6 5 38 38 H 1.08999 * 109.47059 * 209.99668 * 7 6 5 39 39 H 1.00903 * 111.00043 * 303.95917 * 8 7 6 40 40 H 1.09000 * 109.46643 * 299.99506 * 11 9 8 41 41 H 1.08992 * 109.47377 * 59.99598 * 11 9 8 42 42 H 1.08994 * 109.47109 * 180.02562 * 11 9 8 43 43 H 1.09007 * 109.47024 * 182.09468 * 12 9 8 44 44 H 1.08999 * 109.47201 * 62.09585 * 12 9 8 45 45 H 1.08998 * 110.24702 * 281.95599 * 14 13 12 46 46 H 1.09001 * 110.24713 * 43.87160 * 14 13 12 47 47 H 1.08992 * 110.97109 * 143.57887 * 15 14 13 48 48 H 1.09003 * 110.96674 * 267.34007 * 15 14 13 49 49 H 1.08998 * 110.46933 * 216.51659 * 16 15 14 50 50 H 1.09000 * 110.60360 * 93.93790 * 16 15 14 51 51 H 1.08004 * 119.99748 * 179.74637 * 19 6 5 52 52 H 1.08994 * 109.47044 * 300.00211 * 22 21 20 53 53 H 1.08999 * 109.46827 * 60.00514 * 22 21 20 54 54 H 1.09001 * 109.47421 * 300.00360 * 23 22 21 55 55 H 1.08999 * 109.46774 * 60.00713 * 23 22 21 56 56 H 0.96996 * 120.00765 * 180.02562 * 27 26 24 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 8 1.4290 0.0000 0.0000 3 6 2.0464 1.2116 0.0000 4 6 1.2933 2.3748 0.0005 5 6 1.9225 3.6057 0.0011 6 6 3.3031 3.6797 0.0005 7 6 3.9838 5.0242 0.0005 8 7 4.2055 5.4609 1.3855 9 6 4.8694 6.7707 1.4248 10 1 4.4606 7.4088 0.6414 11 6 6.3719 6.5877 1.2011 12 6 4.6324 7.4219 2.7887 13 7 3.1923 7.5515 3.0244 14 6 2.3676 8.6579 2.5248 15 6 0.9072 8.1644 2.6259 16 6 1.0098 7.1132 3.7585 17 6 2.4564 6.6718 3.7257 18 8 2.8893 5.6704 4.2553 19 6 4.0603 2.5233 -0.0053 20 6 3.4360 1.2862 0.0000 21 8 4.1795 0.1476 -0.0006 22 6 5.6003 0.3005 -0.0006 23 6 6.2618 -1.0792 -0.0006 24 6 7.7601 -0.9179 0.0000 25 1 8.1969 0.0699 0.0007 26 6 8.5733 -2.0312 -0.0005 27 7 8.0415 -3.2340 -0.0019 28 14 10.4306 -1.8314 0.0010 29 1 10.7729 -0.3863 0.0019 30 1 10.9969 -2.4754 1.2133 31 1 10.9985 -2.4733 -1.2117 32 1 -0.3634 0.5138 0.8900 33 1 -0.3633 0.5139 -0.8900 34 1 -0.3634 -1.0276 0.0005 35 1 0.2147 2.3196 0.0005 36 1 1.3344 4.5115 0.0011 37 1 3.3538 5.7505 -0.5129 38 1 4.9418 4.9464 -0.5136 39 1 3.3390 5.4787 1.9022 40 1 6.7808 5.9496 1.9845 41 1 6.5408 6.1237 0.2294 42 1 6.8646 7.5595 1.2306 43 1 5.0935 8.4096 2.8048 44 1 5.0740 6.8021 3.5690 45 1 2.5100 9.5438 3.1436 46 1 2.6203 8.8783 1.4877 47 1 0.2390 8.9783 2.9068 48 1 0.5861 7.7025 1.6923 49 1 0.3490 6.2701 3.5570 50 1 0.7699 7.5629 4.7220 51 1 5.1386 2.5839 -0.0098 52 1 5.9065 0.8507 -0.8903 53 1 5.9066 0.8498 0.8896 54 1 5.9555 -1.6290 0.8893 55 1 5.9555 -1.6288 -0.8907 56 1 8.6135 -4.0174 -0.0018 RHF calculation, no. of doubly occupied orbitals= 72 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=WATER ZINC000069922095.mol2 56 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 07:00:35 Heat of formation + Delta-G solvation = -111.436923 kcal Electronic energy + Delta-G solvation = -35060.264039 eV Core-core repulsion = 30547.791092 eV Total energy + Delta-G solvation = -4512.472948 eV No. of doubly occupied orbitals = 72 Molecular weight (most abundant/longest-lived isotopes) = 376.220 amu Computer time = 1.54 seconds Orbital eigenvalues (eV) -41.60874 -41.41064 -39.49654 -38.83910 -37.70259 -37.07210 -35.50122 -33.96045 -32.55259 -32.06185 -31.01008 -30.19468 -28.44051 -28.22039 -27.27184 -27.05584 -25.56502 -23.99468 -22.83110 -22.54232 -22.01452 -21.60205 -20.34786 -19.53470 -18.66133 -17.99900 -17.79786 -17.62178 -17.45780 -17.18417 -16.70730 -16.60386 -16.50883 -16.36327 -15.72490 -15.65851 -15.27520 -15.10956 -14.74301 -14.62241 -14.46241 -14.35770 -14.29796 -14.10962 -13.85901 -13.57893 -13.34739 -13.21082 -13.10925 -13.01337 -12.91936 -12.73966 -12.68660 -12.55290 -12.51765 -12.25465 -12.15813 -12.04140 -11.92036 -11.89586 -11.76623 -11.61626 -11.60484 -11.40457 -10.80468 -10.64641 -9.71993 -9.66596 -9.35823 -8.75032 -8.27899 -6.42175 0.28108 0.40725 1.38555 1.40735 1.54288 1.58427 2.05719 2.25333 2.39985 2.49264 2.63375 3.05814 3.29303 3.49913 3.66969 3.71546 3.81618 3.87738 4.02794 4.07634 4.10999 4.20300 4.22354 4.22525 4.23907 4.35452 4.38280 4.40754 4.41185 4.46387 4.54976 4.60189 4.67986 4.70414 4.72004 4.74208 4.77647 4.82359 4.91337 4.91704 5.03060 5.11507 5.16947 5.22155 5.23946 5.29241 5.36762 5.43422 5.54899 5.61775 5.66078 5.76356 6.03796 6.08162 6.18333 6.22402 6.39962 6.68690 6.72538 7.19692 7.27674 8.85817 Molecular weight = 376.22amu Principal moments of inertia in cm(-1) A = 0.007893 B = 0.002349 C = 0.001933 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 3546.773081 B =11916.122236 C =14480.777962 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.018 3.982 2 O -0.342 6.342 3 C 0.075 3.925 4 C -0.164 4.164 5 C -0.110 4.110 6 C -0.085 4.085 7 C 0.112 3.888 8 N -0.697 5.697 9 C 0.072 3.928 10 H 0.065 0.935 11 C -0.140 4.140 12 C 0.122 3.878 13 N -0.622 5.622 14 C 0.095 3.905 15 C -0.120 4.120 16 C -0.155 4.155 17 C 0.494 3.506 18 O -0.578 6.578 19 C -0.162 4.162 20 C 0.068 3.932 21 O -0.336 6.336 22 C 0.087 3.913 23 C 0.031 3.969 24 C -0.545 4.545 25 H 0.069 0.931 26 C 0.000 4.000 27 N -0.931 5.931 28 Si 0.979 3.021 29 H -0.258 1.258 30 H -0.276 1.276 31 H -0.273 1.273 32 H 0.068 0.932 33 H 0.072 0.928 34 H 0.110 0.890 35 H 0.167 0.833 36 H 0.157 0.843 37 H 0.062 0.938 38 H 0.076 0.924 39 H 0.339 0.661 40 H 0.043 0.957 41 H 0.052 0.948 42 H 0.083 0.917 43 H 0.115 0.885 44 H 0.062 0.938 45 H 0.090 0.910 46 H 0.093 0.907 47 H 0.118 0.882 48 H 0.087 0.913 49 H 0.101 0.899 50 H 0.116 0.884 51 H 0.125 0.875 52 H 0.059 0.941 53 H 0.051 0.949 54 H -0.006 1.006 55 H -0.002 1.002 56 H 0.267 0.733 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -19.548 33.536 1.473 38.845 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.076 4.076 2 O -0.256 6.256 3 C 0.029 3.971 4 C -0.182 4.182 5 C -0.128 4.128 6 C -0.086 4.086 7 C -0.016 4.016 8 N -0.331 5.331 9 C -0.038 4.038 10 H 0.082 0.918 11 C -0.199 4.199 12 C 0.000 4.000 13 N -0.356 5.356 14 C -0.027 4.027 15 C -0.158 4.158 16 C -0.195 4.195 17 C 0.284 3.716 18 O -0.458 6.458 19 C -0.182 4.182 20 C 0.023 3.977 21 O -0.251 6.251 22 C 0.013 3.987 23 C -0.006 4.006 24 C -0.582 4.582 25 H 0.087 0.913 26 C -0.198 4.198 27 N -0.571 5.571 28 Si 0.802 3.198 29 H -0.181 1.181 30 H -0.201 1.201 31 H -0.199 1.199 32 H 0.087 0.913 33 H 0.091 0.909 34 H 0.128 0.872 35 H 0.184 0.816 36 H 0.174 0.826 37 H 0.080 0.920 38 H 0.095 0.905 39 H 0.160 0.840 40 H 0.062 0.938 41 H 0.071 0.929 42 H 0.102 0.898 43 H 0.133 0.867 44 H 0.080 0.920 45 H 0.108 0.892 46 H 0.111 0.889 47 H 0.136 0.864 48 H 0.106 0.894 49 H 0.119 0.881 50 H 0.134 0.866 51 H 0.143 0.857 52 H 0.077 0.923 53 H 0.069 0.931 54 H 0.012 0.988 55 H 0.017 0.983 56 H 0.076 0.924 Dipole moment (debyes) X Y Z Total from point charges -18.602 32.124 1.807 37.166 hybrid contribution 0.445 -0.350 -0.030 0.566 sum -18.158 31.775 1.777 36.640 Atomic orbital electron populations 1.23358 0.76394 1.05345 1.02509 1.86187 1.23575 1.26663 1.89222 1.18628 0.88935 0.86373 1.03135 1.21302 1.01024 0.90328 1.05546 1.21004 0.93264 0.97851 1.00701 1.19688 0.93828 0.94344 1.00711 1.20800 0.99625 0.91295 0.89923 1.60444 1.30452 1.29256 1.12996 1.21960 0.95848 0.88838 0.97198 0.91757 1.21766 0.94051 1.02405 1.01639 1.21638 0.79743 1.02327 0.96325 1.48533 1.08669 1.25169 1.53184 1.22393 0.90132 0.88861 1.01275 1.22682 0.95858 0.98963 0.98292 1.22441 0.91976 1.01988 1.03069 1.21800 0.90910 0.81959 0.76964 1.90793 1.78535 1.28011 1.48433 1.21043 1.00504 0.90830 1.05776 1.18920 0.89540 0.86733 1.02553 1.86167 1.20292 1.29793 1.88801 1.22815 0.75933 1.00023 0.99971 1.18939 0.94891 0.90979 0.95779 1.21393 0.90937 0.93711 1.52185 0.91252 1.26091 1.05154 0.94532 0.94060 1.70138 1.37657 0.98666 1.50600 0.85195 0.77575 0.79012 0.78010 1.18139 1.20136 1.19878 0.91346 0.90911 0.87159 0.81559 0.82560 0.92030 0.90541 0.84018 0.93801 0.92874 0.89818 0.86655 0.91989 0.89181 0.88890 0.86369 0.89396 0.88080 0.86631 0.85675 0.92285 0.93059 0.98756 0.98342 0.92382 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.02 0.38 9.81 113.37 1.11 1.50 16 2 O -0.34 -12.24 10.28 -79.56 -0.82 -13.06 16 3 C 0.07 2.66 6.70 22.57 0.15 2.81 16 4 C -0.16 -4.16 9.04 22.35 0.20 -3.96 16 5 C -0.11 -2.53 9.69 22.27 0.22 -2.32 16 6 C -0.08 -2.35 5.39 -19.86 -0.11 -2.46 16 7 C 0.11 2.42 5.61 85.56 0.48 2.90 16 8 N -0.70 -16.13 6.29 -495.58 -3.11 -19.25 16 9 C 0.07 1.15 2.92 44.91 0.13 1.28 16 10 H 0.06 0.76 7.87 -2.39 -0.02 0.74 16 11 C -0.14 -2.35 9.36 71.98 0.67 -1.68 16 12 C 0.12 1.79 5.65 86.38 0.49 2.28 16 13 N -0.62 -8.23 2.58 -817.98 -2.11 -10.34 16 14 C 0.09 0.37 7.34 86.83 0.64 1.01 16 15 C -0.12 -0.27 7.20 31.73 0.23 -0.04 16 16 C -0.15 -1.47 6.87 30.67 0.21 -1.26 16 17 C 0.49 9.69 6.50 87.90 0.57 10.26 16 18 O -0.58 -18.54 16.94 -3.01 -0.05 -18.59 16 19 C -0.16 -5.84 8.69 22.35 0.19 -5.65 16 20 C 0.07 2.83 6.70 22.57 0.15 2.98 16 21 O -0.34 -16.76 9.92 -79.63 -0.79 -17.55 16 22 C 0.09 4.35 4.68 71.98 0.34 4.69 16 23 C 0.03 1.87 5.34 29.85 0.16 2.03 16 24 C -0.54 -32.83 9.11 25.60 0.23 -32.60 16 25 H 0.07 3.93 7.74 -2.91 -0.02 3.91 16 26 C 0.00 0.00 5.46 84.65 0.46 0.46 16 27 N -0.93 -59.37 13.29 21.45 0.28 -59.08 16 28 Si 0.98 45.64 29.54 68.60 2.03 47.66 16 29 H -0.26 -11.33 7.11 99.48 0.71 -10.63 16 30 H -0.28 -12.50 7.11 99.48 0.71 -11.79 16 31 H -0.27 -12.26 7.11 99.48 0.71 -11.55 16 32 H 0.07 1.20 7.66 -2.39 -0.02 1.18 16 33 H 0.07 1.20 7.66 -2.39 -0.02 1.18 16 34 H 0.11 2.00 8.14 -2.39 -0.02 1.98 16 35 H 0.17 2.88 6.30 -2.91 -0.02 2.86 16 36 H 0.16 2.32 8.06 -2.91 -0.02 2.29 16 37 H 0.06 0.98 8.01 -2.39 -0.02 0.96 16 38 H 0.08 1.73 7.09 -2.39 -0.02 1.71 16 39 H 0.34 8.37 5.66 -89.69 -0.51 7.87 16 40 H 0.04 0.94 8.14 -2.39 -0.02 0.92 16 41 H 0.05 0.97 7.19 -2.39 -0.02 0.95 16 42 H 0.08 0.93 8.14 -2.39 -0.02 0.91 16 43 H 0.12 0.92 8.14 -2.38 -0.02 0.90 16 44 H 0.06 1.35 7.79 -2.39 -0.02 1.33 16 45 H 0.09 0.09 8.14 -2.39 -0.02 0.08 16 46 H 0.09 0.13 8.14 -2.39 -0.02 0.11 16 47 H 0.12 -0.31 8.14 -2.39 -0.02 -0.32 16 48 H 0.09 0.14 8.14 -2.39 -0.02 0.12 16 49 H 0.10 1.00 8.14 -2.39 -0.02 0.98 16 50 H 0.12 0.81 8.14 -2.39 -0.02 0.79 16 51 H 0.13 4.65 6.30 -2.91 -0.02 4.64 16 52 H 0.06 2.75 7.65 -2.39 -0.02 2.74 16 53 H 0.05 2.50 7.68 -2.39 -0.02 2.48 16 54 H -0.01 -0.39 8.14 -2.39 -0.02 -0.41 16 55 H 0.00 -0.12 8.14 -2.39 -0.02 -0.14 16 56 H 0.27 15.77 8.31 -92.71 -0.77 15.00 16 Total: -1.00 -88.50 450.90 2.32 -86.18 By element: Atomic # 1 Polarization: 21.41 SS G_CDS: 0.36 Total: 21.78 kcal Atomic # 6 Polarization: -24.29 SS G_CDS: 6.53 Total: -17.75 kcal Atomic # 7 Polarization: -83.73 SS G_CDS: -4.94 Total: -88.67 kcal Atomic # 8 Polarization: -47.54 SS G_CDS: -1.66 Total: -49.20 kcal Atomic # 14 Polarization: 45.64 SS G_CDS: 2.03 Total: 47.66 kcal Total: -88.50 2.32 -86.18 kcal The number of atoms in the molecule is 56 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. ZINC000069922095.mol2 56 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 -25.260 kcal (2) G-P(sol) polarization free energy of solvation -88.499 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -113.759 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy 2.322 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -86.177 kcal (6) G-S(sol) free energy of system = (1) + (5) -111.437 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.54 seconds