Wall clock time and date at job start Wed Apr 1 2020 04:57:13
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.53000 * 1
3 3 H 1.08997 * 110.63716 * 2 1
4 4 C 1.54325 * 110.87403 * 236.75172 * 2 1 3
5 5 N 1.47290 * 107.30599 * 219.38314 * 4 2 1
6 6 C 1.36938 * 125.79079 * 178.92205 * 5 4 2
7 7 H 1.08005 * 119.99558 * 180.02562 * 6 5 4
8 8 C 1.38811 * 120.00146 * 0.02562 * 6 5 4
9 9 N 1.31620 * 120.00086 * 359.97438 * 8 6 5
10 10 Si 1.86800 * 119.99779 * 180.02562 * 8 6 5
11 11 H 1.48496 * 109.47294 * 89.99594 * 10 8 6
12 12 H 1.48503 * 109.46797 * 209.99511 * 10 8 6
13 13 H 1.48500 * 109.47035 * 329.99310 * 10 8 6
14 14 C 1.47425 * 108.55742 * 358.92448 * 5 4 2
15 15 C 1.54908 * 104.92961 * 24.15407 * 14 5 4
16 16 H 1.09008 * 111.00494 * 81.01141 * 15 14 5
17 17 C 1.52999 * 111.00368 * 204.85960 * 15 14 5
18 18 C 1.52997 * 109.47275 * 293.54646 * 17 15 14
19 19 N 1.46502 * 109.46940 * 179.97438 * 18 17 15
20 20 C 1.34780 * 119.99919 * 180.02562 * 19 18 17
21 21 O 1.21285 * 120.00027 * 0.02562 * 20 19 18
22 22 C 1.50701 * 120.00258 * 180.02562 * 20 19 18
23 23 H 1.09001 * 109.67191 * 59.47879 * 22 20 19
24 24 C 1.53398 * 109.68371 * 179.97438 * 22 20 19
25 25 C 1.53111 * 109.44812 * 169.37801 * 24 22 20
26 26 C 1.53354 * 109.53468 * 290.00340 * 25 24 22
27 27 C 1.50862 * 108.41867 * 49.97671 * 26 25 24
28 28 N 1.33572 * 128.03527 * 163.24799 * 27 26 25
29 29 C 1.30751 * 109.14750 * 180.02562 * 28 27 26
30 30 N 1.35312 * 108.76665 * 0.02562 * 29 28 27
31 31 H 0.96998 * 126.50110 * 179.97438 * 30 29 28
32 32 C 1.35030 * 123.98338 * 343.30562 * 27 26 25
33 33 H 1.08999 * 109.47394 * 296.85917 * 1 2 3
34 34 H 1.09001 * 109.47080 * 56.85884 * 1 2 3
35 35 H 1.08999 * 109.46882 * 176.85698 * 1 2 3
36 36 H 1.09004 * 109.88239 * 338.79110 * 4 2 1
37 37 H 1.08995 * 109.88239 * 99.80032 * 4 2 1
38 38 H 0.96999 * 120.00175 * 180.02562 * 9 8 6
39 39 H 1.09002 * 110.36696 * 265.26054 * 14 5 4
40 40 H 1.08995 * 110.36608 * 143.04891 * 14 5 4
41 41 H 1.08994 * 109.47238 * 53.55041 * 17 15 14
42 42 H 1.09004 * 109.46626 * 173.54596 * 17 15 14
43 43 H 1.08996 * 109.47431 * 299.99824 * 18 17 15
44 44 H 1.08998 * 109.46915 * 60.00403 * 18 17 15
45 45 H 0.96993 * 120.00063 * 0.02562 * 19 18 17
46 46 H 1.09000 * 109.47131 * 289.37370 * 24 22 20
47 47 H 1.08999 * 109.47331 * 49.38441 * 24 22 20
48 48 H 1.09004 * 109.46046 * 50.01970 * 25 24 22
49 49 H 1.08998 * 109.46662 * 169.99160 * 25 24 22
50 50 H 1.09005 * 109.67229 * 169.71105 * 26 25 24
51 51 H 1.08999 * 109.67504 * 290.24053 * 26 25 24
52 52 H 1.07997 * 125.61748 * 179.97438 * 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.5300 0.0000 0.0000
3 1 1.9142 1.0200 0.0000
4 6 2.0799 -0.7906 -1.2059
5 7 3.2516 -1.5489 -0.7352
6 6 4.0155 -2.3951 -1.4939
7 1 4.8575 -2.9060 -1.0508
8 6 3.7075 -2.5961 -2.8325
9 7 2.6810 -1.9738 -3.3723
10 14 4.7492 -3.7510 -3.8672
11 1 4.1913 -5.1249 -3.7890
12 1 4.7412 -3.2930 -5.2798
13 1 6.1432 -3.7529 -3.3553
14 6 3.4652 -1.2446 0.6914
15 6 2.0789 -0.7933 1.2151
16 1 1.4468 -1.6529 1.4384
17 6 2.2236 0.1112 2.4406
18 6 2.7977 -0.6969 3.6060
19 7 2.9368 0.1694 4.7792
20 6 3.4267 -0.3299 5.9313
21 8 3.7528 -1.4962 5.9967
22 6 3.5692 0.5611 7.1384
23 1 2.5947 0.9643 7.4140
24 6 4.1447 -0.2445 8.3101
25 6 4.5346 0.7099 9.4421
26 6 5.7579 1.5330 9.0204
27 6 5.4889 2.1228 7.6581
28 7 6.1284 3.1291 7.0561
29 6 5.5912 3.3376 5.8824
30 7 4.5800 2.4565 5.7038
31 1 4.0082 2.3848 4.9236
32 6 4.5170 1.6898 6.8268
33 1 -0.3634 0.4643 0.9168
34 1 -0.3633 0.5618 -0.8605
35 1 -0.3633 -1.0261 -0.0563
36 1 1.3191 -1.4773 -1.5770
37 1 2.3750 -0.1017 -1.9974
38 1 2.4656 -2.1146 -4.3075
39 1 4.1929 -0.4409 0.8030
40 1 3.8001 -2.1352 1.2232
41 1 2.8950 0.9370 2.2055
42 1 1.2461 0.5055 2.7185
43 1 2.1263 -1.5227 3.8413
44 1 3.7752 -1.0911 3.3281
45 1 2.6764 1.1023 4.7269
46 1 3.3945 -0.9484 8.6703
47 1 5.0263 -0.7919 7.9765
48 1 3.7018 1.3801 9.6551
49 1 4.7751 0.1342 10.3359
50 1 5.9271 2.3340 9.7402
51 1 6.6358 0.8887 8.9747
52 1 5.9028 4.0898 5.1728
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
ZINC001103425652.mol2 52
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 04:57:13
Heat of formation + Delta-G solvation = 22.825811 kcal
Electronic energy + Delta-G solvation = -30650.552470 eV
Core-core repulsion = 26650.414153 eV
Total energy + Delta-G solvation = -4000.138316 eV
No. of doubly occupied orbitals = 66
Molecular weight
(most abundant/longest-lived isotopes) = 346.211 amu
Computer time = 1.20 seconds
Orbital eigenvalues (eV)
-42.73797 -40.01217 -38.27854 -36.79296 -35.39182 -33.73727 -33.16561 -31.37183
-30.72146 -30.33892 -29.24996 -27.02522 -25.76470 -24.82455 -23.84369 -22.58961
-22.30503 -21.76244 -21.04106 -19.79738 -19.31622 -19.25334 -18.34191 -17.25832
-16.89170 -16.13073 -15.92963 -15.51724 -15.08869 -15.01075 -14.70049 -14.47595
-14.20487 -14.06227 -13.85864 -13.62789 -13.38878 -13.12351 -12.92196 -12.62382
-12.42083 -12.29838 -12.02761 -11.78096 -11.57987 -11.48461 -11.28232 -11.27013
-11.23296 -11.15342 -10.81171 -10.79139 -10.51177 -10.19755 -10.10377 -10.01689
-9.94187 -9.88350 -9.69650 -9.54164 -9.32098 -8.61266 -8.50166 -6.83334
-5.71849 -3.04964 1.59614 2.30633 2.67925 2.76824 2.91909 3.39751
3.44958 3.46619 4.09206 4.16872 4.35976 4.42300 4.49417 4.64313
4.76969 4.84624 4.92646 4.97793 5.00160 5.11383 5.16539 5.27343
5.41280 5.42635 5.53771 5.55488 5.60053 5.61979 5.66192 5.69339
5.76461 5.81262 5.88205 5.88773 5.99573 6.02634 6.12064 6.16971
6.19197 6.26768 6.39761 6.45016 6.57990 6.60606 6.73583 6.94599
6.97349 7.11748 7.56179 7.61547 7.74027 7.82512 7.88931 8.25304
9.21695 9.82186 10.48564 11.40031
Molecular weight = 346.21amu
Principal moments of inertia in cm(-1)
A = 0.015742 B = 0.002137 C = 0.002092
Principal moments of inertia in units of 10**(-40)*gram-cm**2
A = 1778.257910 B =13100.809109 C =13377.916413
Net atomic charges, atomic populations, and dipole contributions
calculated with CM2
Atom NO. Type Charge No. of electrons
1 C -0.133 4.133
2 C -0.114 4.114
3 H 0.063 0.937
4 C 0.183 3.817
5 N -0.607 5.607
6 C -0.227 4.227
7 H 0.077 0.923
8 C -0.009 4.009
9 N -0.923 5.923
10 Si 0.911 3.089
11 H -0.289 1.289
12 H -0.291 1.291
13 H -0.281 1.281
14 C 0.140 3.860
15 C -0.095 4.095
16 H 0.077 0.923
17 C -0.116 4.116
18 C 0.116 3.884
19 N -0.733 5.733
20 C 0.514 3.486
21 O -0.523 6.523
22 C -0.027 4.027
23 H 0.109 0.891
24 C -0.097 4.097
25 C -0.117 4.117
26 C -0.037 4.037
27 C 0.021 3.979
28 N -0.486 5.486
29 C 0.189 3.811
30 N -0.560 5.560
31 H 0.416 0.584
32 C -0.029 4.029
33 H 0.042 0.958
34 H 0.051 0.949
35 H 0.054 0.946
36 H 0.048 0.952
37 H 0.040 0.960
38 H 0.252 0.748
39 H 0.015 0.985
40 H 0.037 0.963
41 H 0.066 0.934
42 H 0.065 0.935
43 H 0.071 0.929
44 H 0.075 0.925
45 H 0.389 0.611
46 H 0.079 0.921
47 H 0.081 0.919
48 H 0.067 0.933
49 H 0.077 0.923
50 H 0.079 0.921
51 H 0.078 0.922
52 H 0.211 0.789
For the charges calculated by CM2:
Dipole moment (debyes) X Y Z Total
from point charges -2.421 11.529 23.533 26.317
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.191 4.191
2 C -0.133 4.133
3 H 0.081 0.919
4 C 0.058 3.942
5 N -0.333 5.333
6 C -0.355 4.355
7 H 0.095 0.905
8 C -0.205 4.205
9 N -0.570 5.570
10 Si 0.742 3.258
11 H -0.216 1.216
12 H -0.217 1.217
13 H -0.207 1.207
14 C 0.013 3.987
15 C -0.115 4.115
16 H 0.096 0.904
17 C -0.154 4.154
18 C -0.007 4.007
19 N -0.386 5.386
20 C 0.299 3.701
21 O -0.397 6.397
22 C -0.049 4.049
23 H 0.127 0.873
24 C -0.135 4.135
25 C -0.155 4.155
26 C -0.075 4.075
27 C -0.110 4.110
28 N -0.204 5.204
29 C -0.094 4.094
30 N -0.163 5.163
31 H 0.256 0.744
32 C -0.145 4.145
33 H 0.061 0.939
34 H 0.070 0.930
35 H 0.073 0.927
36 H 0.066 0.934
37 H 0.058 0.942
38 H 0.062 0.938
39 H 0.033 0.967
40 H 0.055 0.945
41 H 0.084 0.916
42 H 0.084 0.916
43 H 0.089 0.911
44 H 0.093 0.907
45 H 0.223 0.777
46 H 0.098 0.902
47 H 0.099 0.901
48 H 0.085 0.915
49 H 0.096 0.904
50 H 0.097 0.903
51 H 0.096 0.904
52 H 0.228 0.772
Dipole moment (debyes) X Y Z Total
from point charges -1.473 11.021 23.681 26.162
hybrid contribution -0.588 -0.415 -1.762 1.904
sum -2.061 10.606 21.919 24.437
Atomic orbital electron populations
1.21677 0.95521 1.01614 1.00253 1.22478 0.94959 0.98857 0.97017
0.91866 1.21246 0.89068 0.93107 0.90771 1.44522 1.32425 1.51588
1.04729 1.19066 1.07440 1.17722 0.91301 0.90497 1.24797 0.98480
1.00366 0.96892 1.69960 1.28891 1.45677 1.12513 0.86125 0.79585
0.80403 0.79716 1.21591 1.21681 1.20702 1.21595 0.93198 0.97030
0.86843 1.22328 0.98989 0.97579 0.92618 0.90449 1.21752 1.00718
0.98686 0.94287 1.21567 1.02316 0.93288 0.83539 1.46228 1.66278
1.13798 1.12269 1.20527 0.76986 0.88580 0.83981 1.90746 1.48089
1.17243 1.83657 1.19934 0.98523 0.91679 0.94790 0.87335 1.21439
1.02441 0.97023 0.92609 1.21633 0.97880 0.97970 0.98058 1.20000
0.98430 0.98562 0.90540 1.20530 0.95467 0.96909 0.98050 1.70030
1.27555 1.16981 1.05880 1.25132 0.89679 0.97442 0.97132 1.43068
1.29121 1.27923 1.16220 0.74447 1.20437 1.04348 0.99663 0.90049
0.93877 0.92953 0.92713 0.93418 0.94210 0.93804 0.96721 0.94490
0.91571 0.91587 0.91078 0.90711 0.77728 0.90188 0.90072 0.91467
0.90411 0.90252 0.90353 0.77220
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.13 -1.96 9.17 37.16 0.34 -1.62 16
2 C -0.11 -1.96 3.01 -88.79 -0.27 -2.23 16
3 H 0.06 1.03 8.04 -51.93 -0.42 0.62 16
4 C 0.18 4.34 4.98 -2.92 -0.01 4.32 16
5 N -0.61 -15.44 3.40 -169.86 -0.58 -16.02 16
6 C -0.23 -6.39 10.28 -15.06 -0.15 -6.54 16
7 H 0.08 2.10 7.83 -52.48 -0.41 1.69 16
8 C -0.01 -0.26 5.49 -17.43 -0.10 -0.36 16
9 N -0.92 -27.15 10.33 55.49 0.57 -26.58 16
10 Si 0.91 22.32 29.55 -169.99 -5.02 17.30 16
11 H -0.29 -7.05 7.11 56.52 0.40 -6.65 16
12 H -0.29 -7.15 7.11 56.52 0.40 -6.75 16
13 H -0.28 -6.74 7.11 56.52 0.40 -6.34 16
14 C 0.14 2.80 5.89 -2.53 -0.01 2.79 16
15 C -0.10 -1.51 2.52 -88.47 -0.22 -1.73 16
16 H 0.08 1.28 8.05 -51.92 -0.42 0.86 16
17 C -0.12 -1.30 4.84 -26.73 -0.13 -1.43 16
18 C 0.12 1.34 5.52 -4.04 -0.02 1.32 16
19 N -0.73 -6.75 5.10 -60.29 -0.31 -7.06 16
20 C 0.51 5.70 7.13 -10.98 -0.08 5.62 16
21 O -0.52 -8.00 15.82 5.55 0.09 -7.91 16
22 C -0.03 -0.22 2.90 -92.74 -0.27 -0.49 16
23 H 0.11 0.66 8.14 -51.93 -0.42 0.23 16
24 C -0.10 -0.78 4.68 -26.47 -0.12 -0.90 16
25 C -0.12 -0.74 5.90 -26.49 -0.16 -0.89 16
26 C -0.04 -0.30 6.30 -27.62 -0.17 -0.47 16
27 C 0.02 0.22 6.98 -82.90 -0.58 -0.36 16
28 N -0.49 -5.54 11.25 -10.74 -0.12 -5.66 16
29 C 0.19 1.61 13.05 -89.92 -1.17 0.43 16
30 N -0.56 -4.27 5.90 -10.29 -0.06 -4.33 16
31 H 0.42 2.21 6.85 -40.82 -0.28 1.93 16
32 C -0.03 -0.26 6.32 -81.86 -0.52 -0.78 16
33 H 0.04 0.52 8.14 -51.93 -0.42 0.09 16
34 H 0.05 0.76 8.14 -51.93 -0.42 0.33 16
35 H 0.05 0.85 8.07 -51.93 -0.42 0.43 16
36 H 0.05 1.26 7.03 -51.93 -0.36 0.90 16
37 H 0.04 1.06 7.51 -51.93 -0.39 0.67 16
38 H 0.25 7.20 8.31 -40.82 -0.34 6.86 16
39 H 0.02 0.30 8.00 -51.93 -0.42 -0.12 16
40 H 0.04 0.75 8.06 -51.93 -0.42 0.33 16
41 H 0.07 0.70 8.10 -51.93 -0.42 0.28 16
42 H 0.07 0.63 8.14 -51.93 -0.42 0.20 16
43 H 0.07 0.89 8.14 -51.93 -0.42 0.47 16
44 H 0.07 1.03 7.45 -51.93 -0.39 0.65 16
45 H 0.39 2.40 6.58 -40.82 -0.27 2.14 16
46 H 0.08 0.59 8.14 -51.93 -0.42 0.17 16
47 H 0.08 0.82 7.73 -51.93 -0.40 0.42 16
48 H 0.07 0.38 8.12 -51.93 -0.42 -0.05 16
49 H 0.08 0.37 8.14 -51.93 -0.42 -0.05 16
50 H 0.08 0.59 8.14 -51.93 -0.42 0.16 16
51 H 0.08 0.64 8.14 -51.93 -0.42 0.21 16
52 H 0.21 1.01 8.06 -52.49 -0.42 0.58 16
LS Contribution 404.72 15.07 6.10 6.10
Total: -1.00 -35.43 404.72 -11.78 -47.21
By element:
Atomic # 1 Polarization: 9.07 SS G_CDS: -8.80 Total: 0.27 kcal
Atomic # 6 Polarization: 0.33 SS G_CDS: -3.65 Total: -3.32 kcal
Atomic # 7 Polarization: -59.15 SS G_CDS: -0.49 Total: -59.65 kcal
Atomic # 8 Polarization: -8.00 SS G_CDS: 0.09 Total: -7.91 kcal
Atomic # 14 Polarization: 22.32 SS G_CDS: -5.02 Total: 17.30 kcal
Total LS contribution 6.10 Total: 6.10 kcal
Total: -35.43 -11.78 -47.21 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.
ZINC001103425652.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 70.038 kcal
(2) G-P(sol) polarization free energy of solvation -35.434 kcal
(3) G-ENP(sol) elect.-nuc.-pol. free energy of system 34.604 kcal
(4) G-CDS(sol) cavity-dispersion-solvent structure free energy -11.778 kcal
(5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -47.212 kcal
(6) G-S(sol) free energy of system = (1) + (5) 22.826 kcal
1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED
Total computer time = 1.20 seconds