Wall clock time and date at job start Tue Mar 31 2020 21:18:16 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.31004 * 1 3 3 C 1.50708 * 119.99226 * 2 1 4 4 C 1.53003 * 109.46870 * 125.00265 * 3 2 1 5 5 C 1.50702 * 109.53670 * 299.99585 * 4 3 2 6 6 O 1.21283 * 119.99562 * 239.82074 * 5 4 3 7 7 N 1.34772 * 120.00533 * 59.82683 * 5 4 3 8 8 C 1.46498 * 120.00148 * 179.97438 * 7 5 4 9 9 H 1.09002 * 109.52726 * 325.08489 * 8 7 5 10 10 C 1.53036 * 109.49480 * 85.16248 * 8 7 5 11 11 C 1.53043 * 109.53585 * 181.31612 * 10 8 7 12 12 C 1.53190 * 109.31461 * 298.63391 * 11 10 8 13 13 N 1.46931 * 108.77598 * 54.63484 * 12 11 10 14 14 C 1.36936 * 120.63056 * 126.41335 * 13 12 11 15 15 H 1.07997 * 119.99636 * 209.03877 * 14 13 12 16 16 C 1.38808 * 120.00109 * 29.03639 * 14 13 12 17 17 N 1.31625 * 119.99625 * 15.05995 * 16 14 13 18 18 Si 1.86798 * 120.00303 * 195.05562 * 16 14 13 19 19 H 1.48501 * 109.47082 * 0.02562 * 18 16 14 20 20 H 1.48490 * 109.47451 * 120.00213 * 18 16 14 21 21 H 1.48504 * 109.47126 * 240.00622 * 18 16 14 22 22 C 1.46930 * 118.73743 * 306.39537 * 13 12 11 23 23 H 1.09000 * 109.58345 * 293.67936 * 22 13 12 24 24 C 1.52998 * 109.58550 * 173.39399 * 22 13 12 25 25 C 1.53172 * 109.61973 * 60.21374 * 4 3 2 26 26 C 1.53095 * 109.16580 * 176.89064 * 25 4 3 27 27 O 1.42909 * 109.41247 * 302.40097 * 26 25 4 28 28 C 1.42899 * 114.10049 * 61.15249 * 27 26 25 29 29 C 1.53092 * 109.41448 * 298.83892 * 28 27 26 30 30 H 1.07993 * 120.00141 * 0.02562 * 1 2 3 31 31 H 1.07998 * 119.99018 * 179.97438 * 1 2 3 32 32 H 1.08000 * 120.00476 * 179.97438 * 2 1 3 33 33 H 1.08993 * 109.47181 * 5.00003 * 3 2 1 34 34 H 1.09000 * 109.46720 * 244.99529 * 3 2 1 35 35 H 0.96998 * 120.00054 * 359.97438 * 7 5 4 36 36 H 1.08997 * 109.46521 * 61.29940 * 10 8 7 37 37 H 1.09000 * 109.45292 * 301.32541 * 10 8 7 38 38 H 1.09000 * 109.49845 * 178.68384 * 11 10 8 39 39 H 1.08995 * 109.52835 * 58.60263 * 11 10 8 40 40 H 1.09003 * 109.59078 * 294.84978 * 12 11 10 41 41 H 1.09001 * 109.70863 * 174.49336 * 12 11 10 42 42 H 0.97010 * 119.99482 * 179.97438 * 17 16 14 43 43 H 1.08999 * 109.46974 * 59.99824 * 24 22 13 44 44 H 1.08996 * 109.46857 * 180.02562 * 24 22 13 45 45 H 1.09005 * 109.46875 * 299.99761 * 24 22 13 46 46 H 1.09000 * 109.52817 * 296.80276 * 25 4 3 47 47 H 1.09006 * 109.52295 * 56.99085 * 25 4 3 48 48 H 1.08993 * 109.48454 * 182.41546 * 26 25 4 49 49 H 1.09000 * 109.48766 * 62.37775 * 26 25 4 50 50 H 1.09002 * 109.47677 * 178.85970 * 28 27 26 51 51 H 1.08994 * 109.48657 * 58.83235 * 28 27 26 52 52 H 1.08994 * 109.52632 * 297.70076 * 29 28 27 53 53 H 1.08997 * 109.57023 * 177.54582 * 29 28 27 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.3100 0.0000 0.0000 3 6 2.0634 1.3053 0.0000 4 6 3.0350 1.3333 1.1816 5 6 2.2653 1.1985 2.4702 6 8 2.3255 2.0745 3.3069 7 7 1.5102 0.1046 2.6930 8 6 0.7625 -0.0265 3.9460 9 1 1.3345 0.4188 4.7600 10 6 -0.5830 0.6908 3.8152 11 6 -1.3511 0.5824 5.1346 12 6 -1.6183 -0.8937 5.4455 13 7 -0.3406 -1.6191 5.4297 14 6 0.0412 -2.3831 6.5001 15 1 1.0898 -2.5232 6.7167 16 6 -0.9197 -2.9774 7.3065 17 7 -2.1659 -3.0625 6.8917 18 14 -0.4441 -3.6564 8.9805 19 1 1.0067 -3.4375 9.2095 20 1 -1.2219 -2.9600 10.0364 21 1 -0.7389 -5.1112 9.0265 22 6 0.5193 -1.5095 4.2434 23 1 0.0281 -1.9794 3.3914 24 6 1.8555 -2.2045 4.5126 25 6 4.0308 0.1762 1.0562 26 6 5.0448 0.2596 2.2002 27 8 5.7033 1.5273 2.1598 28 6 4.8275 2.6451 2.3198 29 6 3.8049 2.6575 1.1806 30 1 -0.5400 0.9352 -0.0004 31 1 -0.5398 -0.9354 -0.0004 32 1 1.8501 -0.9353 0.0004 33 1 1.3584 2.1317 0.0896 34 1 2.6212 1.4027 -0.9314 35 1 1.4619 -0.5958 2.0238 36 1 -1.1648 0.2291 3.0175 37 1 -0.4129 1.7413 3.5794 38 1 -2.2991 1.1134 5.0483 39 1 -0.7597 1.0220 5.9376 40 1 -2.2849 -1.3113 4.6909 41 1 -2.0780 -0.9846 6.4296 42 1 -2.8373 -3.4782 7.4551 43 1 1.6792 -3.2566 4.7364 44 1 2.4920 -2.1217 3.6317 45 1 2.3473 -1.7299 5.3617 46 1 3.4963 -0.7721 1.1116 47 1 4.5522 0.2462 0.1016 48 1 5.7815 -0.5363 2.0917 49 1 4.5277 0.1489 3.1533 50 1 5.4077 3.5676 2.2975 51 1 4.3074 2.5658 3.2744 52 1 4.3216 2.7777 0.2285 53 1 3.1086 3.4838 1.3240 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 ZINC001486186395.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 21:18:16 Heat of formation + Delta-G solvation = -50.770387 kcal Electronic energy + Delta-G solvation = -32340.281902 eV Core-core repulsion = 28405.574126 eV Total energy + Delta-G solvation = -3934.707776 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 1.28 seconds Orbital eigenvalues (eV) -41.09672 -38.64832 -37.38734 -36.18900 -34.31710 -32.87200 -32.54527 -31.69447 -31.10555 -29.19056 -27.84543 -26.28254 -25.05123 -24.00531 -23.79855 -22.54335 -21.27383 -20.74078 -19.97535 -19.68972 -18.89191 -16.75389 -16.53734 -16.19468 -16.04017 -15.73793 -15.14795 -15.10049 -14.84282 -14.53091 -14.17180 -14.01248 -13.75623 -13.68364 -13.40300 -13.16259 -13.00060 -12.69179 -12.53027 -12.36192 -12.27125 -11.78477 -11.69947 -11.61755 -11.54034 -11.32793 -11.13065 -10.93574 -10.78152 -10.67349 -10.61321 -10.48328 -10.29758 -10.13564 -10.02943 -9.86687 -9.64948 -9.56220 -9.35533 -9.15437 -8.59832 -8.45134 -6.69028 -5.81063 -3.25406 2.05479 2.82122 3.35406 3.38791 3.44436 3.46719 3.83926 4.28946 4.45322 4.52268 4.57557 4.80521 4.87558 4.96180 5.00494 5.16077 5.28295 5.33068 5.37642 5.43624 5.46320 5.50658 5.52257 5.58188 5.63771 5.69671 5.73265 5.79086 5.79694 5.88344 5.93112 5.96183 6.03948 6.12130 6.20195 6.31642 6.38112 6.42365 6.52474 6.54373 6.64174 6.73592 6.74878 6.85783 6.96646 7.06346 7.15236 7.39146 7.66630 7.79336 7.81056 8.41405 8.48981 9.18550 9.81280 10.37376 11.27388 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.015813 B = 0.003660 C = 0.003360 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1770.260465 B = 7648.271203 C = 8331.821804 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.168 4.168 2 C -0.178 4.178 3 C -0.077 4.077 4 C -0.076 4.076 5 C 0.521 3.479 6 O -0.543 6.543 7 N -0.711 5.711 8 C 0.130 3.870 9 H 0.084 0.916 10 C -0.113 4.113 11 C -0.141 4.141 12 C 0.169 3.831 13 N -0.582 5.582 14 C -0.251 4.251 15 H 0.071 0.929 16 C 0.004 3.996 17 N -0.902 5.902 18 Si 0.901 3.099 19 H -0.277 1.277 20 H -0.289 1.289 21 H -0.290 1.290 22 C 0.169 3.831 23 H 0.025 0.975 24 C -0.137 4.137 25 C -0.142 4.142 26 C 0.057 3.943 27 O -0.385 6.385 28 C 0.049 3.951 29 C -0.127 4.127 30 H 0.096 0.904 31 H 0.102 0.898 32 H 0.114 0.886 33 H 0.073 0.927 34 H 0.082 0.918 35 H 0.395 0.605 36 H 0.053 0.947 37 H 0.047 0.953 38 H 0.050 0.950 39 H 0.055 0.945 40 H 0.018 0.982 41 H 0.089 0.911 42 H 0.253 0.747 43 H 0.064 0.936 44 H 0.036 0.964 45 H 0.054 0.946 46 H 0.078 0.922 47 H 0.082 0.918 48 H 0.095 0.905 49 H 0.053 0.947 50 H 0.089 0.911 51 H 0.079 0.921 52 H 0.074 0.926 53 H 0.079 0.921 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 9.100 7.667 -15.821 19.797 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.205 4.205 2 C -0.196 4.196 3 C -0.115 4.115 4 C -0.079 4.079 5 C 0.307 3.693 6 O -0.421 6.421 7 N -0.364 5.364 8 C 0.024 3.976 9 H 0.102 0.898 10 C -0.152 4.152 11 C -0.180 4.180 12 C 0.044 3.956 13 N -0.305 5.305 14 C -0.380 4.380 15 H 0.088 0.912 16 C -0.195 4.195 17 N -0.547 5.547 18 Si 0.731 3.269 19 H -0.202 1.202 20 H -0.216 1.216 21 H -0.217 1.217 22 C 0.061 3.939 23 H 0.043 0.957 24 C -0.194 4.194 25 C -0.180 4.180 26 C -0.020 4.020 27 O -0.304 6.304 28 C -0.028 4.028 29 C -0.165 4.165 30 H 0.115 0.885 31 H 0.120 0.880 32 H 0.132 0.868 33 H 0.091 0.909 34 H 0.101 0.899 35 H 0.233 0.767 36 H 0.072 0.928 37 H 0.066 0.934 38 H 0.069 0.931 39 H 0.074 0.926 40 H 0.036 0.964 41 H 0.107 0.893 42 H 0.063 0.937 43 H 0.083 0.917 44 H 0.055 0.945 45 H 0.073 0.927 46 H 0.097 0.903 47 H 0.101 0.899 48 H 0.113 0.887 49 H 0.071 0.929 50 H 0.107 0.893 51 H 0.097 0.903 52 H 0.092 0.908 53 H 0.098 0.902 Dipole moment (debyes) X Y Z Total from point charges 9.194 8.198 -15.124 19.506 hybrid contribution -0.392 -0.765 0.493 0.991 sum 8.802 7.433 -14.632 18.623 Atomic orbital electron populations 1.23963 0.95844 1.02367 0.98334 1.22903 0.96220 0.98343 1.02153 1.20751 0.97474 0.96181 0.97068 1.20913 0.95446 0.97684 0.93889 1.20280 0.79180 0.82091 0.87790 1.90635 1.64328 1.40391 1.46701 1.46439 1.47525 1.21432 1.20961 1.21761 0.91745 0.99013 0.85040 0.89828 1.21726 0.95946 0.99060 0.98481 1.22185 0.98491 1.01151 0.96211 1.20845 0.87103 0.85756 1.01899 1.44037 1.19032 1.53747 1.13675 1.19212 0.92335 1.23590 1.02876 0.91175 1.25064 0.96167 0.98435 0.99788 1.70034 1.01628 1.46905 1.36112 0.86243 0.78000 0.78282 0.84381 1.20234 1.21571 1.21748 1.20538 0.91966 0.90281 0.91123 0.95739 1.21823 0.96121 1.01719 0.99749 1.22300 0.97066 0.97736 1.00931 1.22581 0.95729 0.86338 0.97306 1.87883 1.34489 1.14043 1.93962 1.22838 0.91911 0.88966 0.99066 1.22004 0.98386 0.96012 1.00095 0.88520 0.87968 0.86833 0.90877 0.89934 0.76718 0.92840 0.93385 0.93123 0.92600 0.96408 0.89251 0.93722 0.91725 0.94510 0.92672 0.90309 0.89930 0.88669 0.92910 0.89286 0.90317 0.90780 0.90209 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 17. 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.17 -1.57 13.91 29.01 0.40 -1.17 16 2 C -0.18 -1.55 5.82 -36.02 -0.21 -1.76 16 3 C -0.08 -0.68 4.73 -27.88 -0.13 -0.81 16 4 C -0.08 -0.80 0.51 -155.48 -0.08 -0.88 16 5 C 0.52 7.38 5.08 -10.99 -0.06 7.33 16 6 O -0.54 -9.63 13.50 5.55 0.08 -9.55 16 7 N -0.71 -10.01 3.80 -53.69 -0.20 -10.21 16 8 C 0.13 2.30 1.99 -67.81 -0.14 2.17 16 9 H 0.08 1.73 7.58 -51.93 -0.39 1.34 16 10 C -0.11 -2.01 5.46 -26.69 -0.15 -2.16 16 11 C -0.14 -2.76 6.08 -26.61 -0.16 -2.92 16 12 C 0.17 3.97 5.29 -3.71 -0.02 3.95 16 13 N -0.58 -13.93 2.99 -164.77 -0.49 -14.43 16 14 C -0.25 -6.59 8.60 -15.06 -0.13 -6.72 16 15 H 0.07 1.82 6.21 -52.49 -0.33 1.49 16 16 C 0.00 0.11 4.96 -17.43 -0.09 0.02 16 17 N -0.90 -25.53 11.14 55.49 0.62 -24.91 16 18 Si 0.90 21.39 29.58 -169.99 -5.03 16.37 16 19 H -0.28 -6.46 7.11 56.52 0.40 -6.06 16 20 H -0.29 -6.93 7.11 56.52 0.40 -6.53 16 21 H -0.29 -6.95 7.11 56.52 0.40 -6.55 16 22 C 0.17 3.32 3.03 -67.60 -0.21 3.11 16 23 H 0.03 0.50 8.14 -51.93 -0.42 0.07 16 24 C -0.14 -2.48 7.83 37.16 0.29 -2.19 16 25 C -0.14 -1.34 4.52 -26.59 -0.12 -1.46 16 26 C 0.06 0.65 6.62 37.21 0.25 0.89 16 27 O -0.38 -4.95 10.77 -35.23 -0.38 -5.33 16 28 C 0.05 0.59 6.55 37.20 0.24 0.83 16 29 C -0.13 -1.34 4.56 -26.59 -0.12 -1.46 16 30 H 0.10 0.92 7.84 -52.49 -0.41 0.51 16 31 H 0.10 0.92 8.06 -52.49 -0.42 0.50 16 32 H 0.11 0.87 6.77 -52.49 -0.36 0.51 16 33 H 0.07 0.69 7.79 -51.93 -0.40 0.28 16 34 H 0.08 0.59 8.14 -51.93 -0.42 0.17 16 35 H 0.40 4.85 5.00 -40.82 -0.20 4.64 16 36 H 0.05 0.91 8.14 -51.93 -0.42 0.49 16 37 H 0.05 0.85 8.14 -51.93 -0.42 0.43 16 38 H 0.05 0.90 8.14 -51.93 -0.42 0.48 16 39 H 0.06 1.15 8.14 -51.93 -0.42 0.73 16 40 H 0.02 0.43 8.14 -51.93 -0.42 0.01 16 41 H 0.09 2.35 6.06 -51.93 -0.31 2.04 16 42 H 0.25 6.94 8.31 -40.82 -0.34 6.60 16 43 H 0.06 1.25 7.62 -51.93 -0.40 0.86 16 44 H 0.04 0.54 7.80 -51.93 -0.41 0.13 16 45 H 0.05 1.06 7.45 -51.93 -0.39 0.67 16 46 H 0.08 0.67 6.11 -51.93 -0.32 0.35 16 47 H 0.08 0.67 8.14 -51.93 -0.42 0.24 16 48 H 0.10 0.90 8.14 -51.93 -0.42 0.47 16 49 H 0.05 0.70 7.38 -51.93 -0.38 0.32 16 50 H 0.09 0.93 8.14 -51.93 -0.42 0.51 16 51 H 0.08 1.16 5.80 -51.93 -0.30 0.86 16 52 H 0.07 0.66 8.14 -51.93 -0.42 0.23 16 53 H 0.08 0.87 8.03 -51.93 -0.42 0.46 16 LS Contribution 392.00 15.07 5.91 5.91 Total: -1.00 -29.98 392.00 -9.14 -39.12 By element: Atomic # 1 Polarization: 15.47 SS G_CDS: -9.22 Total: 6.25 kcal Atomic # 6 Polarization: -2.80 SS G_CDS: -0.42 Total: -3.22 kcal Atomic # 7 Polarization: -49.47 SS G_CDS: -0.08 Total: -49.55 kcal Atomic # 8 Polarization: -14.58 SS G_CDS: -0.30 Total: -14.88 kcal Atomic # 14 Polarization: 21.39 SS G_CDS: -5.03 Total: 16.37 kcal Total LS contribution 5.91 Total: 5.91 kcal Total: -29.98 -9.14 -39.12 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. ZINC001486186395.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 -11.646 kcal (2) G-P(sol) polarization free energy of solvation -29.980 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -41.626 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.144 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -39.124 kcal (6) G-S(sol) free energy of system = (1) + (5) -50.770 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.28 seconds