Wall clock time and date at job start Tue Mar 31 2020 18:00:18                  
 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.52999 *                                 1
    3       3  H            1.08999 *   109.47305 *                   2   1
    4       4  C            1.52993 *   109.47275 *   120.00423 *     2   1   3
    5       5  O            1.42903 *   109.47133 *    54.99736 *     4   2   1
    6       6  C            1.42906 *   113.99959 *   180.02562 *     5   4   2
    7       7  H            1.09004 *   110.43212 *   328.90410 *     6   5   4
    8       8  C            1.54811 *   110.48603 *   206.22034 *     6   5   4
    9       9  C            1.54997 *   102.08899 *   220.44002 *     8   6   5
   10      10  O            1.44382 *   103.49183 *    37.02498 *     9   8   6
   11      11  C            1.43747 *   106.97323 *   319.93922 *    10   9   8
   12      12  N            1.46506 *   109.46935 *   240.00403 *     2   1   3
   13      13  C            1.34772 *   120.00149 *    84.99972 *    12   2   1
   14      14  O            1.21586 *   120.00137 *   359.97438 *    13  12   2
   15      15  N            1.34776 *   120.00002 *   180.02562 *    13  12   2
   16      16  C            1.39246 *   120.00162 *   175.71136 *    15  13  12
   17      17  C            1.39481 *   119.83362 *   215.02326 *    16  15  13
   18      18  C            1.36375 *   120.16744 *   179.76196 *    17  16  15
   19      19  C            1.40867 *   119.82699 *     0.48658 *    18  17  16
   20      20  C            1.41155 *   120.15997 *   179.78792 *    19  18  17
   21      21  H            1.08001 *   119.99832 *   214.75897 *    20  19  18
   22      22  C            1.39896 *   120.00522 *    34.76534 *    20  19  18
   23      23  N            1.30881 *   119.99622 *   355.55802 *    22  20  19
   24      24  Si           1.86799 *   120.00225 *   175.55594 *    22  20  19
   25      25  H            1.48496 *   109.47197 *   359.97438 *    24  22  20
   26      26  H            1.48503 *   109.46553 *   119.99962 *    24  22  20
   27      27  H            1.48495 *   109.47239 *   239.99636 *    24  22  20
   28      28  C            1.40866 *   119.67248 *   359.75952 *    19  18  17
   29      29  C            1.36368 *   119.83113 *     0.02562 *    28  19  18
   30      30  H            1.09008 *   109.46895 *   299.99998 *     1   2   3
   31      31  H            1.09000 *   109.47226 *    59.99152 *     1   2   3
   32      32  H            1.08999 *   109.47305 *   180.02562 *     1   2   3
   33      33  H            1.09001 *   109.47286 *   294.99777 *     4   2   1
   34      34  H            1.09002 *   109.47017 *   174.99608 *     4   2   1
   35      35  H            1.08994 *   110.90611 *   338.63663 *     8   6   5
   36      36  H            1.08997 *   110.90623 *   102.24168 *     8   6   5
   37      37  H            1.09000 *   110.61818 *   278.51873 *     9   8   6
   38      38  H            1.09001 *   110.75184 *   155.58833 *     9   8   6
   39      39  H            1.09002 *   109.87807 *   267.08899 *    11  10   9
   40      40  H            1.08994 *   109.88049 *   145.93901 *    11  10   9
   41      41  H            0.97003 *   119.99736 *   265.00167 *    12   2   1
   42      42  H            0.97001 *   120.00082 *   355.71186 *    15  13  12
   43      43  H            1.08002 *   119.91282 *   359.97438 *    17  16  15
   44      44  H            1.08003 *   120.08516 *   180.27346 *    18  17  16
   45      45  H            0.97001 *   120.00153 *   180.02562 *    23  22  20
   46      46  H            1.08002 *   120.08444 *   180.02562 *    28  19  18
   47      47  H            1.08006 *   119.91880 *   179.97438 *    29  28  19

 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.8933    1.0276    0.0000
     4         6        2.0400   -0.7213    1.2491
     5         8        1.4702   -2.0304    1.3093
     6         6        1.8836   -2.7927    2.4452
     7         1        2.8957   -2.5163    2.7410
     8         6        1.7999   -4.3082    2.1403
     9         6        1.2351   -4.8755    3.4675
    10         8        0.2998   -3.8625    3.8959
    11         6        0.9043   -2.5884    3.6172
    12         7        2.0183   -0.6906   -1.1963
    13         6        2.1470   -0.0182   -2.3572
    14         8        1.8574    1.1614   -2.4119
    15         7        2.5957   -0.6536   -3.4578
    16         6        2.6443    0.0165   -4.6775
    17         6        2.3832   -0.6754   -5.8601
    18         6        2.4241   -0.0248   -7.0580
    19         6        2.7411    1.3471   -7.0997
    20         6        2.7898    2.0264   -8.3361
    21         1        3.5103    2.8160   -8.4904
    22         6        1.9048    1.6770   -9.3617
    23         7        1.0909    0.6639   -9.2067
    24        14        1.8846    2.6565  -10.9521
    25         1        2.8929    3.7441  -10.8765
    26         1        2.2082    1.7591  -12.0902
    27         1        0.5364    3.2453  -11.1538
    28         6        3.0086    2.0378   -5.9015
    29         6        2.9595    1.3748   -4.7109
    30         1       -0.3633    0.5139    0.8901
    31         1       -0.3634    0.5140   -0.8899
    32         1       -0.3634   -1.0276    0.0005
    33         1        1.7517   -0.1584    2.1369
    34         1        3.1263   -0.7996    1.2054
    35         1        1.1172   -4.5013    1.3129
    36         1        2.7880   -4.7180    1.9313
    37         1        2.0310   -4.9984    4.2020
    38         1        0.7232   -5.8225    3.2959
    39         1        1.4453   -2.2350    4.4950
    40         1        0.1365   -1.8672    3.3373
    41         1        2.2489   -1.6318   -1.1528
    42         1        2.8849   -1.5777   -3.4003
    43         1        2.1447   -1.7282   -5.8268
    44         1        2.2174   -0.5608   -7.9726
    45         1        0.4774    0.4213   -9.9179
    46         1        3.2519    3.0898   -5.9257
    47         1        3.1648    1.9028   -3.7914

          RHF calculation, no. of doubly occupied orbitals= 63


 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                                       
 ZINC000122690281.mol2 47                                                        
                                                                                 

     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 18:00:18                  

 Heat of formation + Delta-G solvation     =     -96.409009 kcal
 Electronic energy + Delta-G solvation     =  -27950.709633 eV
 Core-core repulsion                       =   23905.501202 eV
 Total energy + Delta-G solvation          =   -4045.208431 eV
 No. of doubly occupied orbitals           =      63
 Molecular weight
  (most abundant/longest-lived isotopes)   =     334.173     amu
 Computer time                             =       2.02     seconds

               Orbital eigenvalues (eV)

 -40.24576 -39.83934 -38.07688 -37.25605 -35.22051 -34.23148 -32.92663 -31.83354
 -30.66114 -29.48785 -29.11381 -28.51681 -26.69274 -24.86882 -23.28537 -22.13381
 -21.78608 -21.14877 -20.55629 -20.04601 -18.06304 -17.71729 -17.20986 -17.04501
 -16.65286 -16.13975 -15.54715 -15.35860 -15.11116 -14.90496 -14.54799 -14.44829
 -13.89539 -13.71422 -13.41719 -13.08818 -12.95108 -12.65701 -12.59744 -12.51810
 -12.19645 -12.06377 -11.89809 -11.64363 -11.61113 -11.30929 -11.17906 -11.13273
 -11.05370 -10.90997 -10.68675 -10.35034 -10.02621  -9.86860  -9.48509  -9.40741
  -9.20230  -8.88480  -8.81600  -7.48112  -6.77425  -6.55912  -4.28035   2.71238
   2.99451   3.11457   3.13066   3.14420   3.18745   3.22848   3.46947   3.75832
   4.00809   4.02655   4.36930   4.46050   4.55610   4.63983   4.67419   4.76178
   4.81799   4.88938   4.98801   5.07274   5.21067   5.24428   5.31659   5.39938
   5.46660   5.54486   5.68643   5.72834   5.97485   6.25061   6.31830   6.33380
   6.49841   6.65302   6.76486   6.82839   6.90706   7.10877   7.21966   7.33109
   7.47904   7.71233   7.92645   8.03745   8.15804   8.36928   8.53922   8.70069
   8.98965   8.99313   9.15961  10.35388

          Molecular weight =  334.17amu

           Principal moments of inertia in cm(-1)
          A =    0.031179   B =    0.001845   C =    0.001805

           Principal moments of inertia in units of 10**(-40)*gram-cm**2
          A =  897.822465   B =15169.245727   C =15506.202129

              Net atomic charges, atomic populations, and dipole contributions
                                   calculated with CM2 

         Atom NO.   Type          Charge        No. of electrons
           1         C            -0.161           4.161
           2         C             0.168           3.832
           3         H             0.100           0.900
           4         C             0.051           3.949
           5         O            -0.366           6.366
           6         C             0.056           3.944
           7         H             0.077           0.923
           8         C            -0.153           4.153
           9         C             0.046           3.954
          10         O            -0.374           6.374
          11         C             0.024           3.976
          12         N            -0.745           5.745
          13         C             0.695           3.305
          14         O            -0.579           6.579
          15         N            -0.668           5.668
          16         C             0.026           3.974
          17         C            -0.124           4.124
          18         C            -0.140           4.140
          19         C             0.096           3.904
          20         C            -0.462           4.462
          21         H             0.075           0.925
          22         C             0.069           3.931
          23         N            -0.814           5.814
          24         Si            0.914           3.086
          25         H            -0.263           1.263
          26         H            -0.275           1.275
          27         H            -0.275           1.275
          28         C            -0.203           4.203
          29         C            -0.087           4.087
          30         H             0.059           0.941
          31         H             0.064           0.936
          32         H             0.064           0.936
          33         H             0.065           0.935
          34         H             0.058           0.942
          35         H             0.093           0.907
          36         H             0.087           0.913
          37         H             0.062           0.938
          38         H             0.097           0.903
          39         H             0.062           0.938
          40         H             0.094           0.906
          41         H             0.395           0.605
          42         H             0.391           0.609
          43         H             0.090           0.910
          44         H             0.132           0.868
          45         H             0.276           0.724
          46         H             0.095           0.905
          47         H             0.109           0.891

For the charges calculated by CM2:

 Dipole moment (debyes)    X         Y         Z       Total
 from point charges     2.329   -10.357    22.653    25.017


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.219           4.219
           2         C             0.063           3.937
           3         H             0.118           0.882
           4         C            -0.026           4.026
           5         O            -0.285           6.285
           6         C            -0.003           4.003
           7         H             0.095           0.905
           8         C            -0.190           4.190
           9         C            -0.030           4.030
          10         O            -0.292           6.292
          11         C            -0.054           4.054
          12         N            -0.402           5.402
          13         C             0.396           3.604
          14         O            -0.456           6.456
          15         N            -0.314           5.314
          16         C            -0.067           4.067
          17         C            -0.145           4.145
          18         C            -0.160           4.160
          19         C             0.094           3.906
          20         C            -0.493           4.493
          21         H             0.093           0.907
          22         C            -0.146           4.146
          23         N            -0.442           5.442
          24         Si            0.740           3.260
          25         H            -0.186           1.186
          26         H            -0.201           1.201
          27         H            -0.201           1.201
          28         C            -0.224           4.224
          29         C            -0.108           4.108
          30         H             0.078           0.922
          31         H             0.083           0.917
          32         H             0.083           0.917
          33         H             0.083           0.917
          34         H             0.076           0.924
          35         H             0.111           0.889
          36         H             0.106           0.894
          37         H             0.080           0.920
          38         H             0.115           0.885
          39         H             0.081           0.919
          40         H             0.112           0.888
          41         H             0.229           0.771
          42         H             0.224           0.776
          43         H             0.108           0.892
          44         H             0.149           0.851
          45         H             0.087           0.913
          46         H             0.113           0.887
          47         H             0.126           0.874

 Dipole moment (debyes)     X         Y         Z       Total
 from point charges      1.529    -9.474    22.445    24.411
 hybrid contribution     0.396    -0.373    -0.512     0.748
 sum                     1.925    -9.848    21.932    24.119

          Atomic orbital electron populations

   1.22111   0.94862   1.02693   1.02270   1.20768   0.94903   0.96321   0.81728
   0.88214   1.22980   0.97459   0.84013   0.98173   1.88017   1.69997   1.26355
   1.44138   1.23529   0.97474   0.92398   0.86878   0.90528   1.23222   1.03207
   0.92995   0.99618   1.23554   0.92774   0.92861   0.93858   1.88919   1.41358
   1.14822   1.84146   1.23762   0.97411   0.87474   0.96708   1.45126   1.72610
   1.13547   1.08916   1.16034   0.78601   0.85427   0.80329   1.90901   1.56719
   1.14125   1.83895   1.42121   1.66909   1.15235   1.07091   1.17080   1.14016
   0.93133   0.82477   1.20298   1.03757   0.97974   0.92435   1.21175   1.03888
   0.94674   0.96285   1.17721   0.89255   0.91344   0.92290   1.19508   1.16275
   1.16886   0.96609   0.90655   1.25867   0.94413   0.94647   0.99651   1.69777
   1.17766   1.23503   1.33121   0.86163   0.78891   0.80165   0.80823   1.18632
   1.20072   1.20057   1.20802   1.11168   0.98221   0.92169   1.20422   0.99748
   0.93054   0.97536   0.92238   0.91724   0.91722   0.91705   0.92354   0.88878
   0.89434   0.92008   0.88452   0.91944   0.88765   0.77106   0.77592   0.89195
   0.85063   0.91304   0.88679   0.87353

   Number of geometries                                      1
   Number of calculations of the screened coulomb radii      1
   The total number of SCF iterations                       16.
    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    -1.87     9.29      37.16       0.35    -1.53     16
   2       C     0.17     1.87     2.98     -67.93      -0.20     1.67     16
   3       H     0.10     1.23     7.58     -51.93      -0.39     0.84     16
   4       C     0.05     0.39     5.54      37.15       0.21     0.60     16
   5       O    -0.37    -3.08     8.81     -35.23      -0.31    -3.39     16
   6       C     0.06     0.37     3.51     -25.21      -0.09     0.28     16
   7       H     0.08     0.38     7.86     -51.93      -0.41    -0.02     16
   8       C    -0.15    -0.91     6.90     -24.72      -0.17    -1.08     16
   9       C     0.05     0.28     7.72      38.21       0.29     0.57     16
  10       O    -0.37    -3.43    11.52     -35.23      -0.41    -3.83     16
  11       C     0.02     0.17     7.52      37.73       0.28     0.45     16
  12       N    -0.74    -9.51     5.06     -58.68      -0.30    -9.81     16
  13       C     0.70    12.41     8.36     -86.92      -0.73    11.69     16
  14       O    -0.58   -12.77    13.42       5.29       0.07   -12.70     16
  15       N    -0.67   -11.98     5.31     -14.10      -0.07   -12.05     16
  16       C     0.03     0.58     6.31     -83.57      -0.53     0.05     16
  17       C    -0.12    -2.88     9.90     -39.79      -0.39    -3.27     16
  18       C    -0.14    -3.66     8.61     -39.26      -0.34    -4.00     16
  19       C     0.10     2.63     5.55    -106.80      -0.59     2.04     16
  20       C    -0.46   -12.53     9.02     -37.86      -0.34   -12.87     16
  21       H     0.08     1.97     7.85     -52.49      -0.41     1.56     16
  22       C     0.07     1.77     5.42     -17.34      -0.09     1.68     16
  23       N    -0.81   -21.57    10.59      55.55       0.59   -20.98     16
  24       Si    0.91    19.36    29.56    -169.99      -5.03    14.33     16
  25       H    -0.26    -5.48     7.11      56.52       0.40    -5.08     16
  26       H    -0.28    -5.67     7.11      56.52       0.40    -5.27     16
  27       H    -0.28    -5.79     7.11      56.52       0.40    -5.39     16
  28       C    -0.20    -5.26     9.74     -39.26      -0.38    -5.64     16
  29       C    -0.09    -2.14     8.61     -39.79      -0.34    -2.48     16
  30       H     0.06     0.57     8.14     -51.92      -0.42     0.15     16
  31       H     0.06     0.94     8.14     -51.93      -0.42     0.52     16
  32       H     0.06     0.76     7.56     -51.93      -0.39     0.37     16
  33       H     0.06     0.40     8.14     -51.93      -0.42    -0.02     16
  34       H     0.06     0.41     8.10     -51.93      -0.42    -0.01     16
  35       H     0.09     0.66     8.09     -51.93      -0.42     0.24     16
  36       H     0.09     0.38     8.14     -51.93      -0.42    -0.05     16
  37       H     0.06     0.28     8.14     -51.93      -0.42    -0.14     16
  38       H     0.10     0.48     8.14     -51.93      -0.42     0.06     16
  39       H     0.06     0.35     8.14     -51.93      -0.42    -0.07     16
  40       H     0.09     0.72     7.97     -51.93      -0.41     0.31     16
  41       H     0.40     4.13     8.40     -40.82      -0.34     3.79     16
  42       H     0.39     5.52     8.83     -40.82      -0.36     5.15     16
  43       H     0.09     1.83     8.06     -52.49      -0.42     1.41     16
  44       H     0.13     3.54     6.13     -52.48      -0.32     3.22     16
  45       H     0.28     6.69     8.30     -40.82      -0.34     6.35     16
  46       H     0.10     2.33     8.06     -52.49      -0.42     1.91     16
  47       H     0.11     2.63     6.22     -52.48      -0.33     2.30     16

 LS Contribution                 386.57      15.07       5.83     5.83

 Total:         -1.00   -32.50   386.57                 -9.85   -42.35

 By element:

 Atomic #  1  Polarization:   19.25  SS G_CDS:  -7.15  Total:    12.10 kcal
 Atomic #  6  Polarization:   -8.77  SS G_CDS:  -3.07  Total:   -11.84 kcal
 Atomic #  7  Polarization:  -43.06  SS G_CDS:   0.22  Total:   -42.85 kcal
 Atomic #  8  Polarization:  -19.27  SS G_CDS:  -0.65  Total:   -19.92 kcal
 Atomic # 14  Polarization:   19.36  SS G_CDS:  -5.03  Total:    14.33 kcal
 Total LS contribution                           5.83  Total:     5.83 kcal

 Total:                      -32.50             -9.85           -42.35 kcal

   The number of atoms in the molecule is     47
   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.

                                                                                 
 ZINC000122690281.mol2 47                                                        

 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                -54.054 kcal
 (2)  G-P(sol) polarization free energy of solvation               -32.504 kcal
 (3)  G-ENP(sol) elect.-nuc.-pol. free energy of system            -86.558 kcal
 (4)  G-CDS(sol) cavity-dispersion-solvent structure free energy    -9.851 kcal
 (5)  G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4)             -42.355 kcal
 (6)  G-S(sol) free energy of system = (1) + (5)                   -96.409 kcal

     1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED

 Total computer time =       2.02 seconds