Wall clock time and date at job start Tue Mar 31 2020 06:45:50 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 N 2 2 C 1.31507 * 1 3 3 Si 1.86806 * 119.99804 * 2 1 4 4 H 1.48499 * 109.46909 * 269.99827 * 3 2 1 5 5 H 1.48501 * 109.47151 * 29.99820 * 3 2 1 6 6 H 1.48496 * 109.47051 * 149.99779 * 3 2 1 7 7 C 1.37877 * 120.00404 * 179.97438 * 2 1 3 8 8 H 1.08000 * 119.99507 * 359.97438 * 7 2 1 9 9 C 1.50700 * 119.99983 * 180.02562 * 7 2 1 10 10 C 1.53006 * 109.46680 * 179.97438 * 9 7 2 11 11 O 1.42906 * 109.46919 * 179.97438 * 10 9 7 12 12 C 1.35897 * 116.99968 * 179.97438 * 11 10 9 13 13 C 1.38804 * 120.02072 * 0.02562 * 12 11 10 14 14 C 1.38173 * 120.08432 * 180.02562 * 13 12 11 15 15 C 1.38164 * 120.12067 * 0.02562 * 14 13 12 16 16 C 1.38880 * 120.04398 * 0.02562 * 15 14 13 17 17 N 1.39942 * 120.04500 * 179.72174 * 16 15 14 18 18 C 1.34770 * 120.00559 * 144.92339 * 17 16 15 19 19 O 1.21590 * 119.99809 * 4.74898 * 18 17 16 20 20 N 1.34771 * 120.00471 * 184.75512 * 18 17 16 21 21 C 1.46929 * 120.60576 * 4.91175 * 20 18 17 22 22 C 1.52774 * 109.00893 * 233.74918 * 21 20 18 23 23 C 1.52991 * 109.35842 * 305.81381 * 22 21 20 24 24 C 1.52997 * 109.89763 * 60.90054 * 23 22 21 25 25 C 1.46918 * 120.61075 * 185.18827 * 20 18 17 26 26 O 1.43310 * 109.64798 * 181.46838 * 23 22 21 27 27 C 1.44206 * 107.26666 * 120.35359 * 26 23 22 28 28 C 1.54920 * 103.08077 * 21.39606 * 27 26 23 29 29 O 1.43704 * 109.40749 * 300.77825 * 23 22 21 30 30 C 1.38584 * 120.01997 * 180.02562 * 12 11 10 31 31 H 0.97001 * 120.00273 * 0.02562 * 1 2 3 32 32 H 1.08992 * 109.47570 * 59.99708 * 9 7 2 33 33 H 1.09000 * 109.47340 * 299.99461 * 9 7 2 34 34 H 1.08999 * 109.47354 * 60.00326 * 10 9 7 35 35 H 1.08993 * 109.46995 * 299.99460 * 10 9 7 36 36 H 1.08006 * 119.95758 * 0.02562 * 13 12 11 37 37 H 1.07998 * 119.93818 * 179.97438 * 14 13 12 38 38 H 1.08007 * 119.97720 * 179.97438 * 15 14 13 39 39 H 0.97003 * 119.99793 * 324.91640 * 17 16 15 40 40 H 1.09004 * 109.58217 * 353.58052 * 21 20 18 41 41 H 1.08992 * 109.58847 * 114.00370 * 21 20 18 42 42 H 1.09002 * 109.56885 * 185.78357 * 22 21 20 43 43 H 1.08998 * 109.40240 * 65.73449 * 22 21 20 44 44 H 1.09001 * 109.50157 * 179.01304 * 24 23 22 45 45 H 1.09001 * 109.52883 * 58.92286 * 24 23 22 46 46 H 1.09002 * 109.58542 * 245.98209 * 25 20 18 47 47 H 1.08997 * 109.58389 * 6.26193 * 25 20 18 48 48 H 1.09002 * 110.69157 * 263.00169 * 27 26 23 49 49 H 1.08992 * 110.84641 * 139.87182 * 27 26 23 50 50 H 1.09002 * 110.97890 * 206.36399 * 28 27 26 51 51 H 1.09001 * 111.00947 * 82.59898 * 28 27 26 52 52 H 1.08001 * 120.06048 * 359.97438 * 30 12 11 Note: An asterisk (*) indicates this parameter was or will be optimized. Cartesian coordinates (angstroms) NO. Atom X Y Z 1 7 0.0000 0.0000 0.0000 2 6 1.3151 0.0000 0.0000 3 14 2.2490 1.6178 0.0000 4 1 2.4965 2.0464 -1.4001 5 1 1.4464 2.6527 0.7000 6 1 3.5465 1.4403 0.7001 7 6 2.0045 -1.1940 0.0005 8 1 1.4645 -2.1293 0.0014 9 6 3.5115 -1.1939 0.0011 10 6 4.0215 -2.6364 0.0023 11 8 5.4506 -2.6364 0.0022 12 6 6.0676 -3.8472 0.0037 13 6 5.3121 -5.0116 0.0047 14 6 5.9380 -6.2434 0.0068 15 6 7.3175 -6.3201 0.0084 16 6 8.0785 -5.1584 0.0086 17 7 9.4758 -5.2361 0.0161 18 6 10.2052 -4.3220 -0.6538 19 8 9.6547 -3.3812 -1.1926 20 7 11.5448 -4.4508 -0.7267 21 6 12.2135 -5.6366 -0.1740 22 6 13.3041 -5.1860 0.7963 23 6 14.2392 -4.2041 0.0876 24 6 13.4532 -2.9732 -0.3682 25 6 12.3637 -3.4052 -1.3548 26 8 15.2848 -3.8036 0.9821 27 6 16.5362 -4.1870 0.3766 28 6 16.1740 -4.3187 -1.1238 29 8 14.8269 -4.8454 -1.0562 30 6 7.4515 -3.9202 0.0043 31 1 -0.4850 0.8400 -0.0004 32 1 3.8745 -0.6798 0.8910 33 1 3.8752 -0.6803 -0.8889 34 1 3.6586 -3.1509 -0.8875 35 1 3.6580 -3.1496 0.8925 36 1 4.2335 -4.9547 0.0038 37 1 5.3479 -7.1479 0.0079 38 1 7.8044 -7.2842 0.0105 39 1 9.9193 -5.9478 0.5037 40 1 11.4866 -6.2523 0.3558 41 1 12.6630 -6.2128 -0.9827 42 1 13.8735 -6.0521 1.1333 43 1 12.8450 -4.6959 1.6548 44 1 14.1285 -2.2705 -0.8564 45 1 12.9925 -2.4945 0.4959 46 1 12.8258 -3.7979 -2.2606 47 1 11.7355 -2.5507 -1.6061 48 1 16.8804 -5.1415 0.7750 49 1 17.2912 -3.4149 0.5244 50 1 16.8418 -5.0183 -1.6266 51 1 16.1885 -3.3463 -1.6162 52 1 8.0410 -3.0152 0.0032 RHF calculation, no. of doubly occupied orbitals= 71 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 ZINC000047736295.mol2 52 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 06:45:50 Heat of formation + Delta-G solvation = -110.857699 kcal Electronic energy + Delta-G solvation = -34139.495118 eV Core-core repulsion = 29491.214193 eV Total energy + Delta-G solvation = -4648.280925 eV No. of doubly occupied orbitals = 71 Molecular weight (most abundant/longest-lived isotopes) = 376.188 amu Computer time = 1.24 seconds Orbital eigenvalues (eV) -41.63235 -40.75300 -39.49982 -37.70184 -36.75147 -35.57092 -34.73695 -33.27981 -32.64312 -31.79283 -31.05723 -29.78948 -29.36520 -28.79078 -25.37744 -25.02354 -24.23614 -23.67150 -22.87317 -22.06626 -21.34907 -20.39793 -19.79227 -19.02674 -18.05191 -17.77366 -17.07743 -16.77442 -16.34797 -16.19397 -16.06856 -15.67144 -15.47889 -15.35529 -15.12527 -14.94790 -14.44380 -14.35309 -14.11991 -14.08823 -14.03962 -13.56624 -13.37029 -13.13216 -12.99810 -12.84315 -12.67003 -12.49747 -12.49230 -12.14337 -11.86777 -11.82552 -11.73807 -11.48366 -11.47659 -10.90343 -10.85689 -10.80477 -10.36359 -10.18790 -10.10632 -9.95770 -9.90114 -9.74644 -9.29483 -9.12006 -8.58123 -8.21563 -7.81328 -6.09542 -4.10066 1.42764 1.66106 2.19639 2.72723 2.93413 3.27979 3.33095 3.36912 3.39908 3.46682 3.88409 3.92100 3.99431 4.04664 4.12426 4.30697 4.42049 4.43642 4.46478 4.67425 4.70640 4.83884 4.86914 4.95552 5.01727 5.04651 5.16098 5.21194 5.21610 5.24946 5.32566 5.38273 5.41298 5.46105 5.47903 5.65836 5.69519 5.88308 5.97560 5.99166 6.03371 6.26572 6.39282 6.47826 6.48804 6.57403 6.65992 6.73331 6.97408 7.18088 7.29552 7.54022 7.67352 8.00437 8.26176 8.43051 8.91278 9.51169 11.02929 Molecular weight = 376.19amu Principal moments of inertia in cm(-1) A = 0.014719 B = 0.001630 C = 0.001501 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1901.909608 B =17178.948284 C =18653.457274 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 N -0.898 5.898 2 C 0.072 3.928 3 Si 0.880 3.120 4 H -0.275 1.275 5 H -0.286 1.286 6 H -0.263 1.263 7 C -0.526 4.526 8 H 0.075 0.925 9 C -0.010 4.010 10 C 0.072 3.928 11 O -0.316 6.316 12 C 0.162 3.838 13 C -0.225 4.225 14 C -0.061 4.061 15 C -0.195 4.195 16 C 0.192 3.808 17 N -0.684 5.684 18 C 0.706 3.294 19 O -0.560 6.560 20 N -0.627 5.627 21 C 0.118 3.882 22 C -0.136 4.136 23 C 0.260 3.740 24 C -0.153 4.153 25 C 0.137 3.863 26 O -0.376 6.376 27 C 0.020 3.980 28 C 0.015 3.985 29 O -0.366 6.366 30 C -0.161 4.161 31 H 0.265 0.735 32 H 0.028 0.972 33 H 0.024 0.976 34 H 0.057 0.943 35 H 0.057 0.943 36 H 0.133 0.867 37 H 0.126 0.874 38 H 0.120 0.880 39 H 0.404 0.596 40 H 0.081 0.919 41 H 0.082 0.918 42 H 0.088 0.912 43 H 0.089 0.911 44 H 0.084 0.916 45 H 0.086 0.914 46 H 0.075 0.925 47 H 0.102 0.898 48 H 0.077 0.923 49 H 0.100 0.900 50 H 0.110 0.890 51 H 0.077 0.923 52 H 0.146 0.854 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 37.197 -16.334 0.521 40.628 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 N -0.537 5.537 2 C -0.130 4.130 3 Si 0.706 3.294 4 H -0.200 1.200 5 H -0.212 1.212 6 H -0.187 1.187 7 C -0.564 4.564 8 H 0.093 0.907 9 C -0.047 4.047 10 C -0.002 4.002 11 O -0.229 6.229 12 C 0.114 3.886 13 C -0.244 4.244 14 C -0.080 4.080 15 C -0.216 4.216 16 C 0.097 3.903 17 N -0.333 5.333 18 C 0.408 3.592 19 O -0.435 6.435 20 N -0.365 5.365 21 C -0.006 4.006 22 C -0.175 4.175 23 C 0.183 3.817 24 C -0.192 4.192 25 C 0.014 3.986 26 O -0.297 6.297 27 C -0.056 4.056 28 C -0.062 4.062 29 O -0.287 6.287 30 C -0.182 4.182 31 H 0.075 0.925 32 H 0.046 0.954 33 H 0.042 0.958 34 H 0.075 0.925 35 H 0.076 0.924 36 H 0.150 0.850 37 H 0.144 0.856 38 H 0.138 0.862 39 H 0.242 0.758 40 H 0.100 0.900 41 H 0.100 0.900 42 H 0.107 0.893 43 H 0.107 0.893 44 H 0.103 0.897 45 H 0.105 0.895 46 H 0.093 0.907 47 H 0.120 0.880 48 H 0.095 0.905 49 H 0.118 0.882 50 H 0.128 0.872 51 H 0.095 0.905 52 H 0.163 0.837 Dipole moment (debyes) X Y Z Total from point charges 36.612 -15.355 0.151 39.702 hybrid contribution 0.669 0.313 0.223 0.772 sum 37.281 -15.042 0.373 40.203 Atomic orbital electron populations 1.69890 1.06262 1.27168 1.50342 1.24922 0.94749 0.99183 0.94116 0.86768 0.80489 0.83286 0.78814 1.20041 1.21209 1.18709 1.20850 0.88360 0.94931 1.52227 0.90719 1.19591 0.96850 0.90329 0.97952 1.23087 0.74654 1.01912 1.00521 1.86141 1.23933 1.25838 1.86998 1.18932 0.91710 0.85223 0.92745 1.21402 1.00851 0.92255 1.09900 1.20646 0.93482 0.97095 0.96794 1.20604 0.93971 0.96300 1.10747 1.17057 0.80724 0.93306 0.99170 1.42690 1.04802 1.32150 1.53697 1.15808 0.81160 0.82741 0.79463 1.90996 1.69482 1.32440 1.50599 1.47279 1.04253 1.22343 1.62655 1.21483 0.94549 0.87451 0.97121 1.22177 0.96583 0.99867 0.98866 1.23241 0.87662 0.91132 0.79674 1.22543 0.98176 0.97079 1.01386 1.21396 0.91314 0.91437 0.94418 1.88926 1.14705 1.81718 1.44383 1.24202 0.87464 1.01943 0.92042 1.24401 0.83981 1.00959 0.96841 1.89270 1.26039 1.66090 1.47275 1.20435 0.91907 0.98124 1.07731 0.92493 0.95405 0.95779 0.92470 0.92441 0.84961 0.85629 0.86186 0.75768 0.90047 0.89976 0.89316 0.89276 0.89744 0.89509 0.90679 0.87971 0.90527 0.88154 0.87215 0.90470 0.83674 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 N -0.90 -27.69 13.96 55.43 0.77 -26.92 16 2 C 0.07 2.13 5.47 -17.82 -0.10 2.03 16 3 Si 0.88 22.17 27.47 -169.99 -4.67 17.50 16 4 H -0.27 -6.82 7.11 56.52 0.40 -6.41 16 5 H -0.29 -7.07 7.11 56.52 0.40 -6.67 16 6 H -0.26 -6.47 6.54 56.52 0.37 -6.10 16 7 C -0.53 -15.15 9.66 -34.44 -0.33 -15.48 16 8 H 0.07 2.26 7.99 -52.49 -0.42 1.84 16 9 C -0.01 -0.25 4.41 -27.88 -0.12 -0.37 16 10 C 0.07 1.56 4.68 37.16 0.17 1.73 16 11 O -0.32 -6.69 10.13 -19.42 -0.20 -6.89 16 12 C 0.16 2.97 6.68 -39.24 -0.26 2.71 16 13 C -0.22 -3.56 9.04 -39.40 -0.36 -3.92 16 14 C -0.06 -0.79 10.03 -39.63 -0.40 -1.19 16 15 C -0.20 -2.33 9.96 -39.37 -0.39 -2.72 16 16 C 0.19 2.65 6.30 -83.71 -0.53 2.12 16 17 N -0.68 -7.43 5.24 -14.24 -0.07 -7.50 16 18 C 0.71 8.93 8.14 -86.93 -0.71 8.22 16 19 O -0.56 -9.57 13.21 5.29 0.07 -9.50 16 20 N -0.63 -5.75 2.98 -178.76 -0.53 -6.28 16 21 C 0.12 0.70 6.17 -3.92 -0.02 0.67 16 22 C -0.14 -0.79 5.49 -26.85 -0.15 -0.94 16 23 C 0.26 1.76 1.03 -26.74 -0.03 1.73 16 24 C -0.15 -1.11 5.45 -26.63 -0.15 -1.25 16 25 C 0.14 1.18 6.19 -3.71 -0.02 1.16 16 26 O -0.38 -2.75 10.87 -59.04 -0.64 -3.39 16 27 C 0.02 0.08 8.08 38.17 0.31 0.39 16 28 C 0.01 0.05 7.91 38.17 0.30 0.36 16 29 O -0.37 -2.36 9.75 -59.04 -0.58 -2.94 16 30 C -0.16 -2.77 8.64 -39.25 -0.34 -3.11 16 31 H 0.27 7.64 8.31 -40.82 -0.34 7.30 16 32 H 0.03 0.70 7.20 -51.93 -0.37 0.33 16 33 H 0.02 0.61 7.57 -51.93 -0.39 0.21 16 34 H 0.06 1.19 7.66 -51.93 -0.40 0.79 16 35 H 0.06 1.18 7.66 -51.93 -0.40 0.78 16 36 H 0.13 2.05 6.29 -52.48 -0.33 1.72 16 37 H 0.13 1.25 8.06 -52.49 -0.42 0.83 16 38 H 0.12 1.05 8.06 -52.48 -0.42 0.62 16 39 H 0.40 2.86 6.30 -40.82 -0.26 2.60 16 40 H 0.08 0.33 5.57 -51.93 -0.29 0.04 16 41 H 0.08 0.47 7.82 -51.93 -0.41 0.07 16 42 H 0.09 0.41 8.14 -51.93 -0.42 -0.01 16 43 H 0.09 0.56 8.14 -51.93 -0.42 0.13 16 44 H 0.08 0.52 8.14 -51.93 -0.42 0.10 16 45 H 0.09 0.71 8.14 -51.93 -0.42 0.29 16 46 H 0.07 0.60 7.82 -51.93 -0.41 0.20 16 47 H 0.10 1.06 7.03 -51.93 -0.37 0.69 16 48 H 0.08 0.24 8.14 -51.93 -0.42 -0.19 16 49 H 0.10 0.23 8.14 -51.93 -0.42 -0.19 16 50 H 0.11 0.19 8.14 -51.93 -0.42 -0.23 16 51 H 0.08 0.24 8.14 -51.93 -0.42 -0.19 16 52 H 0.15 2.75 6.22 -52.49 -0.33 2.42 16 LS Contribution 412.36 15.07 6.21 6.21 Total: -1.00 -36.10 412.36 -10.51 -46.61 By element: Atomic # 1 Polarization: 8.72 SS G_CDS: -7.76 Total: 0.97 kcal Atomic # 6 Polarization: -4.75 SS G_CDS: -3.12 Total: -7.86 kcal Atomic # 7 Polarization: -40.87 SS G_CDS: 0.17 Total: -40.70 kcal Atomic # 8 Polarization: -21.38 SS G_CDS: -1.34 Total: -22.72 kcal Atomic # 14 Polarization: 22.17 SS G_CDS: -4.67 Total: 17.50 kcal Total LS contribution 6.21 Total: 6.21 kcal Total: -36.10 -10.51 -46.61 kcal The number of atoms in the molecule is 52 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. ZINC000047736295.mol2 52 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 -64.252 kcal (2) G-P(sol) polarization free energy of solvation -36.099 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -100.351 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.507 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -46.606 kcal (6) G-S(sol) free energy of system = (1) + (5) -110.858 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.24 seconds