Wall clock time and date at job start Wed Apr 1 2020 00:25:37 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.52997 * 1 3 3 H 1.08994 * 109.89884 * 2 1 4 4 O 1.43886 * 110.07005 * 238.80541 * 2 1 3 5 5 C 1.44425 * 106.99785 * 214.15242 * 4 2 1 6 6 C 1.55118 * 103.56050 * 39.73915 * 5 4 2 7 7 C 1.53588 * 109.90663 * 121.07986 * 2 1 3 8 8 C 1.53107 * 110.50269 * 3.41670 * 7 2 1 9 9 N 1.47284 * 108.14827 * 288.09552 * 8 7 2 10 10 C 1.34776 * 121.04020 * 130.33964 * 9 8 7 11 11 O 1.21589 * 120.00097 * 4.81268 * 10 9 8 12 12 N 1.34779 * 119.99755 * 184.81185 * 10 9 8 13 13 C 1.46495 * 119.99681 * 185.03836 * 12 10 9 14 14 C 1.53005 * 109.46906 * 179.97438 * 13 12 10 15 15 C 1.52998 * 109.46893 * 59.99948 * 14 13 12 16 16 C 1.52997 * 109.47118 * 299.99716 * 14 13 12 17 17 C 1.50701 * 109.47155 * 179.97438 * 14 13 12 18 18 H 1.08002 * 119.99993 * 1.66307 * 17 14 13 19 19 C 1.37875 * 120.00133 * 181.66642 * 17 14 13 20 20 N 1.31518 * 120.00131 * 195.14291 * 19 17 14 21 21 Si 1.86800 * 119.99948 * 15.13566 * 19 17 14 22 22 H 1.48503 * 109.46976 * 83.85953 * 21 19 17 23 23 H 1.48503 * 109.47308 * 203.85620 * 21 19 17 24 24 H 1.48498 * 109.47523 * 323.86055 * 21 19 17 25 25 C 1.47161 * 117.91263 * 310.07142 * 9 8 7 26 26 C 1.53280 * 108.50643 * 51.30674 * 25 9 8 27 27 H 1.08996 * 109.51252 * 67.13111 * 26 25 9 28 28 C 1.53003 * 109.50650 * 187.26850 * 26 25 9 29 29 O 1.42228 * 110.61817 * 240.44638 * 7 2 1 30 30 H 1.08996 * 109.47440 * 306.36457 * 1 2 3 31 31 H 1.09000 * 109.47139 * 66.37189 * 1 2 3 32 32 H 1.09007 * 109.47271 * 186.36883 * 1 2 3 33 33 H 1.09004 * 110.60531 * 158.25054 * 5 4 2 34 34 H 1.09000 * 110.60786 * 281.06595 * 5 4 2 35 35 H 1.09001 * 110.93025 * 205.24562 * 6 5 4 36 36 H 1.08997 * 110.92948 * 81.58279 * 6 5 4 37 37 H 1.09001 * 109.85395 * 168.42331 * 8 7 2 38 38 H 1.08995 * 109.75168 * 47.71432 * 8 7 2 39 39 H 0.96997 * 119.99641 * 5.03220 * 12 10 9 40 40 H 1.09001 * 109.47977 * 59.99853 * 13 12 10 41 41 H 1.09007 * 109.47274 * 299.99522 * 13 12 10 42 42 H 1.09003 * 109.47322 * 299.99683 * 15 14 13 43 43 H 1.09001 * 109.46832 * 180.02562 * 15 14 13 44 44 H 1.08997 * 109.47545 * 59.99676 * 15 14 13 45 45 H 1.09000 * 109.47038 * 181.37937 * 16 14 13 46 46 H 1.08995 * 109.47218 * 301.37919 * 16 14 13 47 47 H 1.09001 * 109.47179 * 61.38342 * 16 14 13 48 48 H 0.96996 * 119.99836 * 179.97438 * 20 19 17 49 49 H 1.08995 * 109.62648 * 171.01274 * 25 9 8 50 50 H 1.08993 * 109.62829 * 291.40739 * 25 9 8 51 51 H 1.09000 * 109.46860 * 60.00280 * 28 26 25 52 52 H 1.09002 * 109.47401 * 179.97438 * 28 26 25 53 53 H 1.08989 * 109.47390 * 300.00227 * 28 26 25 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.9009 1.0249 0.0000 4 8 2.0237 -0.7000 -1.1561 5 6 3.2422 -1.3655 -0.7582 6 6 2.9237 -1.8852 0.6682 7 6 2.0529 -0.7455 1.2368 8 6 0.8895 -1.3306 2.0420 9 7 0.1848 -0.2194 2.7036 10 6 -1.1547 -0.0974 2.6190 11 8 -1.8098 -0.9583 2.0639 12 7 -1.7682 0.9790 3.1498 13 6 -3.2309 1.0593 3.1561 14 6 -3.6668 2.3715 3.8113 15 6 -3.0947 3.5504 3.0215 16 6 -3.1471 2.4207 5.2494 17 6 -5.1715 2.4544 3.8172 18 1 -5.7588 1.6396 3.4203 19 6 -5.7984 3.5705 4.3294 20 7 -7.0823 3.5321 4.6120 21 14 -4.8241 5.1372 4.6220 22 1 -4.7392 5.9108 3.3572 23 1 -5.5030 5.9561 5.6584 24 1 -3.4558 4.7905 5.0830 25 6 0.9886 0.7501 3.4650 26 6 2.1305 1.2447 2.5700 27 1 1.7207 1.8310 1.7476 28 6 3.0846 2.1134 3.3923 29 8 2.8408 0.1186 2.0464 30 1 -0.3634 0.6093 0.8275 31 1 -0.3633 0.4119 -0.9415 32 1 -0.3634 -1.0214 0.1140 33 1 3.4621 -2.1962 -1.4288 34 1 4.0741 -0.6615 -0.7356 35 1 3.8365 -2.0059 1.2515 36 1 2.3636 -2.8193 0.6263 37 1 1.2702 -2.0228 2.7930 38 1 0.2051 -1.8532 1.3738 39 1 -1.2411 1.7003 3.5276 40 1 -3.6386 0.2196 3.7191 41 1 -3.6013 1.0238 2.1315 42 1 -2.0063 3.4906 3.0168 43 1 -3.4048 4.4852 3.4885 44 1 -3.4648 3.5153 1.9969 45 1 -3.4783 3.3444 5.7240 46 1 -3.5351 1.5673 5.8054 47 1 -2.0577 2.3857 5.2435 48 1 -7.5232 4.3172 4.9726 49 1 0.3633 1.5926 3.7600 50 1 1.3978 0.2690 4.3533 51 1 3.4953 1.5253 4.2130 52 1 3.8966 2.4644 2.7555 53 1 2.5419 2.9688 3.7944 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 ZINC000331504950.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:25:37 Heat of formation + Delta-G solvation = -72.810569 kcal Electronic energy + Delta-G solvation = -33040.393340 eV Core-core repulsion = 28914.273778 eV Total energy + Delta-G solvation = -4126.119562 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 1.05 seconds Orbital eigenvalues (eV) -41.60335 -39.59008 -38.04016 -36.88832 -34.84924 -34.01068 -33.12536 -31.75340 -31.44175 -30.75826 -27.64098 -26.75037 -26.43361 -25.79995 -25.02794 -24.17655 -22.78779 -21.94447 -20.54229 -19.43871 -19.25354 -18.12169 -17.75680 -16.44665 -16.13760 -16.06944 -15.82789 -15.62910 -15.53013 -15.27565 -14.59865 -14.36483 -14.27071 -13.99441 -13.90653 -13.22413 -13.03066 -13.00356 -12.88350 -12.75907 -12.59528 -12.49915 -12.27065 -12.16327 -12.08182 -11.96726 -11.79180 -11.57358 -11.42194 -11.30911 -11.24614 -10.71741 -10.48139 -10.32651 -10.26940 -10.23324 -10.08859 -9.89179 -9.79375 -9.18233 -8.78182 -8.68522 -8.48481 -8.31535 -6.17858 -4.14672 2.60560 3.11619 3.19814 3.26850 3.35078 3.35936 3.62734 4.01846 4.11469 4.25378 4.32333 4.39572 4.50452 4.66231 4.70761 4.82375 4.90179 4.99290 5.12671 5.13537 5.24977 5.30114 5.33748 5.35077 5.38922 5.44955 5.48670 5.59333 5.62725 5.71704 5.74660 5.80718 5.92411 6.00070 6.03630 6.18722 6.28944 6.31799 6.42250 6.48513 6.57135 6.60922 6.66208 6.76831 6.87359 6.92026 7.02597 7.19213 7.39986 7.69512 8.12774 8.42785 8.78421 8.88487 9.46242 10.95327 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.019426 B = 0.003157 C = 0.002988 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1441.000068 B = 8867.979924 C = 9368.250257 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.159 4.159 2 C 0.058 3.942 3 H 0.071 0.929 4 O -0.369 6.369 5 C 0.047 3.953 6 C -0.140 4.140 7 C 0.077 3.923 8 C 0.113 3.887 9 N -0.640 5.640 10 C 0.707 3.293 11 O -0.585 6.585 12 N -0.738 5.738 13 C 0.159 3.841 14 C 0.030 3.970 15 C -0.148 4.148 16 C -0.143 4.143 17 C -0.514 4.514 18 H 0.077 0.923 19 C 0.076 3.924 20 N -0.889 5.889 21 Si 0.885 3.115 22 H -0.285 1.285 23 H -0.306 1.306 24 H -0.249 1.249 25 C 0.083 3.917 26 C 0.069 3.931 27 H 0.067 0.933 28 C -0.143 4.143 29 O -0.359 6.359 30 H 0.076 0.924 31 H 0.059 0.941 32 H 0.078 0.922 33 H 0.095 0.905 34 H 0.063 0.937 35 H 0.095 0.905 36 H 0.094 0.906 37 H 0.072 0.928 38 H 0.113 0.887 39 H 0.391 0.609 40 H 0.061 0.939 41 H 0.062 0.938 42 H 0.026 0.974 43 H 0.042 0.958 44 H 0.054 0.946 45 H 0.056 0.944 46 H 0.050 0.950 47 H 0.019 0.981 48 H 0.265 0.735 49 H 0.089 0.911 50 H 0.078 0.922 51 H 0.067 0.933 52 H 0.067 0.933 53 H 0.077 0.923 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 24.607 -5.896 -4.309 25.668 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.216 4.216 2 C 0.000 4.000 3 H 0.089 0.911 4 O -0.288 6.288 5 C -0.029 4.029 6 C -0.178 4.178 7 C 0.036 3.964 8 C -0.011 4.011 9 N -0.379 5.379 10 C 0.410 3.590 11 O -0.463 6.463 12 N -0.395 5.395 13 C 0.036 3.964 14 C 0.028 3.972 15 C -0.204 4.204 16 C -0.201 4.201 17 C -0.552 4.552 18 H 0.095 0.905 19 C -0.126 4.126 20 N -0.526 5.526 21 Si 0.718 3.282 22 H -0.213 1.213 23 H -0.234 1.234 24 H -0.175 1.175 25 C -0.041 4.041 26 C 0.011 3.989 27 H 0.085 0.915 28 C -0.201 4.201 29 O -0.279 6.279 30 H 0.095 0.905 31 H 0.078 0.922 32 H 0.097 0.903 33 H 0.114 0.886 34 H 0.081 0.919 35 H 0.113 0.887 36 H 0.113 0.887 37 H 0.090 0.910 38 H 0.131 0.869 39 H 0.227 0.773 40 H 0.079 0.921 41 H 0.080 0.920 42 H 0.045 0.955 43 H 0.060 0.940 44 H 0.073 0.927 45 H 0.075 0.925 46 H 0.069 0.931 47 H 0.038 0.962 48 H 0.075 0.925 49 H 0.107 0.893 50 H 0.096 0.904 51 H 0.086 0.914 52 H 0.086 0.914 53 H 0.096 0.904 Dipole moment (debyes) X Y Z Total from point charges 23.949 -6.472 -4.789 25.266 hybrid contribution 0.530 1.514 0.585 1.707 sum 24.479 -4.958 -4.204 25.327 Atomic orbital electron populations 1.22367 0.92061 1.04366 1.02791 1.23257 0.96293 0.95794 0.84652 0.91110 1.88962 1.39470 1.67641 1.32765 1.23562 0.88095 0.96622 0.94641 1.23064 1.00921 0.97107 0.96700 1.22867 0.91744 0.91025 0.90786 1.22137 0.93367 0.88749 0.96802 1.47444 1.05081 1.29335 1.56089 1.15761 0.80663 0.83041 0.79561 1.90917 1.63496 1.42619 1.49299 1.45716 1.09362 1.20327 1.64099 1.21403 0.74980 0.99050 1.00945 1.18719 0.97092 0.90010 0.91335 1.22368 0.97403 0.99554 1.01122 1.22308 0.96626 1.02058 0.99089 1.20702 0.87744 1.05031 1.41686 0.90475 1.24963 0.94741 0.98798 0.94130 1.69959 1.10431 1.29750 1.42449 0.87320 0.80166 0.83538 0.77188 1.21318 1.23435 1.17542 1.22333 0.93466 0.93281 0.95058 1.22551 0.92747 0.87599 0.96015 0.91507 1.21738 0.99418 0.98886 1.00058 1.87876 1.36340 1.36882 1.66757 0.90515 0.92185 0.90297 0.88639 0.91885 0.88651 0.88714 0.90990 0.86875 0.77256 0.92067 0.91967 0.95479 0.93986 0.92711 0.92487 0.93137 0.96250 0.92527 0.89330 0.90377 0.91391 0.91374 0.90443 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 14. 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.16 -2.16 7.57 37.16 0.28 -1.88 16 2 C 0.06 0.68 3.24 -26.43 -0.09 0.60 16 3 H 0.07 0.84 6.66 -51.93 -0.35 0.50 16 4 O -0.37 -4.49 11.20 -35.23 -0.39 -4.88 16 5 C 0.05 0.40 7.73 38.28 0.30 0.70 16 6 C -0.14 -1.12 6.44 -24.94 -0.16 -1.28 16 7 C 0.08 0.79 0.97 -89.59 -0.09 0.70 16 8 C 0.11 1.30 4.65 -3.56 -0.02 1.28 16 9 N -0.64 -8.07 1.99 -178.02 -0.35 -8.43 16 10 C 0.71 11.48 5.77 -86.92 -0.50 10.98 16 11 O -0.59 -11.53 14.70 5.29 0.08 -11.45 16 12 N -0.74 -11.94 4.87 -65.22 -0.32 -12.26 16 13 C 0.16 3.14 4.51 -4.04 -0.02 3.12 16 14 C 0.03 0.63 0.50 -155.66 -0.08 0.55 16 15 C -0.15 -2.92 7.49 37.16 0.28 -2.64 16 16 C -0.14 -2.84 7.81 37.16 0.29 -2.55 16 17 C -0.51 -13.61 8.62 -34.44 -0.30 -13.90 16 18 H 0.08 2.25 7.99 -52.49 -0.42 1.83 16 19 C 0.08 2.10 5.34 -17.82 -0.10 2.00 16 20 N -0.89 -26.34 13.96 55.43 0.77 -25.57 16 21 Si 0.88 20.72 22.94 -169.99 -3.90 16.82 16 22 H -0.29 -6.68 7.11 56.52 0.40 -6.28 16 23 H -0.31 -7.35 7.11 56.52 0.40 -6.95 16 24 H -0.25 -5.24 4.54 56.52 0.26 -4.99 16 25 C 0.08 0.82 5.92 -3.54 -0.02 0.80 16 26 C 0.07 0.63 3.11 -26.58 -0.08 0.55 16 27 H 0.07 0.65 6.55 -51.93 -0.34 0.31 16 28 C -0.14 -0.98 9.57 37.16 0.36 -0.62 16 29 O -0.36 -3.65 9.81 -35.23 -0.35 -4.00 16 30 H 0.08 1.11 4.62 -51.93 -0.24 0.87 16 31 H 0.06 0.83 8.14 -51.93 -0.42 0.40 16 32 H 0.08 1.16 5.22 -51.93 -0.27 0.89 16 33 H 0.10 0.63 8.14 -51.93 -0.42 0.20 16 34 H 0.06 0.50 8.14 -51.93 -0.42 0.08 16 35 H 0.09 0.66 8.07 -51.93 -0.42 0.24 16 36 H 0.09 0.67 7.84 -51.93 -0.41 0.27 16 37 H 0.07 0.75 8.14 -51.93 -0.42 0.33 16 38 H 0.11 1.47 5.03 -51.93 -0.26 1.21 16 39 H 0.39 5.51 3.71 -40.82 -0.15 5.36 16 40 H 0.06 1.29 8.14 -51.93 -0.42 0.87 16 41 H 0.06 1.31 8.14 -51.93 -0.42 0.89 16 42 H 0.03 0.44 6.86 -51.93 -0.36 0.08 16 43 H 0.04 0.88 2.15 -51.93 -0.11 0.77 16 44 H 0.05 1.09 8.14 -51.93 -0.42 0.67 16 45 H 0.06 1.19 4.22 -51.93 -0.22 0.97 16 46 H 0.05 1.01 8.14 -51.93 -0.42 0.59 16 47 H 0.02 0.31 6.90 -51.93 -0.36 -0.05 16 48 H 0.26 7.39 8.31 -40.82 -0.34 7.05 16 49 H 0.09 0.83 5.65 -51.93 -0.29 0.54 16 50 H 0.08 0.71 8.14 -51.93 -0.42 0.29 16 51 H 0.07 0.44 8.14 -51.93 -0.42 0.02 16 52 H 0.07 0.45 8.14 -51.93 -0.42 0.03 16 53 H 0.08 0.44 8.14 -51.94 -0.42 0.02 16 LS Contribution 374.94 15.07 5.65 5.65 Total: -1.00 -31.40 374.94 -7.30 -38.70 By element: Atomic # 1 Polarization: 15.56 SS G_CDS: -8.55 Total: 7.01 kcal Atomic # 6 Polarization: -1.65 SS G_CDS: 0.06 Total: -1.59 kcal Atomic # 7 Polarization: -46.36 SS G_CDS: 0.10 Total: -46.26 kcal Atomic # 8 Polarization: -19.67 SS G_CDS: -0.66 Total: -20.33 kcal Atomic # 14 Polarization: 20.72 SS G_CDS: -3.90 Total: 16.82 kcal Total LS contribution 5.65 Total: 5.65 kcal Total: -31.40 -7.30 -38.70 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. ZINC000331504950.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 -34.115 kcal (2) G-P(sol) polarization free energy of solvation -31.397 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -65.512 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -7.299 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -38.696 kcal (6) G-S(sol) free energy of system = (1) + (5) -72.811 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.05 seconds