Wall clock time and date at job start Wed Apr  1 2020 03:01:23                  
 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
 *  SM5.42R  - SM5.42R CALCULATIONS WILL BE PERFORMED
 *           - CHARGE MODEL 2 WILL BE USED (CM2A)
 *           - THE SOLVENT IS WATER               
 *******************************************************************************

   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.31622 *                                 1
    3       3  Si           1.86804 *   119.99736 *                   2   1
    4       4  H            1.48501 *   109.46996 *   274.97922 *     3   2   1
    5       5  H            1.48501 *   109.47087 *    34.98446 *     3   2   1
    6       6  H            1.48494 *   109.46720 *   154.98030 *     3   2   1
    7       7  C            1.38813 *   119.99977 *   180.02562 *     2   1   3
    8       8  H            1.07994 *   119.99764 *   188.94993 *     7   2   1
    9       9  N            1.36940 *   119.99972 *     8.94664 *     7   2   1
   10      10  C            1.47214 *   126.70049 *   194.06017 *     9   7   2
   11      11  C            1.53778 *   106.37898 *   211.30578 *    10   9   7
   12      12  H            1.09000 *   111.70757 *   219.69833 *    11  10   9
   13      13  C            1.54319 *   109.79574 *    95.21973 *    11  10   9
   14      14  C            1.55050 *   103.93821 *   221.68965 *    13  11  10
   15      15  C            1.55053 *   102.10916 *    40.05819 *    14  13  11
   16      16  C            1.54320 *   103.93618 *   319.94341 *    15  14  13
   17      17  N            1.46508 *   111.69296 *   262.79035 *    16  15  14
   18      18  C            1.34774 *   119.99807 *   298.88403 *    17  16  15
   19      19  O            1.21279 *   120.00510 *   359.97438 *    18  17  16
   20      20  C            1.50704 *   119.99882 *   179.97438 *    18  17  16
   21      21  H            1.09002 *   109.47357 *   299.99691 *    20  18  17
   22      22  C            1.53112 *   109.47542 *   179.97438 *    20  18  17
   23      23  C            1.53362 *   109.53085 *   170.00239 *    22  20  18
   24      24  C            1.50865 *   108.41046 *    49.98696 *    23  22  20
   25      25  N            1.33572 *   128.03360 *   163.26721 *    24  23  22
   26      26  C            1.30749 *   109.15081 *   180.08578 *    25  24  23
   27      27  N            1.35317 *   108.76292 *   359.72367 *    26  25  24
   28      28  H            0.97009 *   126.49750 *   180.26777 *    27  26  25
   29      29  C            1.35024 *   123.98800 *   343.29840 *    24  23  22
   30      30  C            1.50649 *   124.45581 *   359.89814 *    29  24  23
   31      31  C            1.47208 *   126.70251 *    14.08548 *     9   7   2
   32      32  H            0.97000 *   120.00347 *   359.97438 *     1   2   3
   33      33  H            1.08993 *   110.07391 *    92.07231 *    10   9   7
   34      34  H            1.09005 *   110.06989 *   330.53135 *    10   9   7
   35      35  H            1.08999 *   110.53616 *   340.30325 *    13  11  10
   36      36  H            1.08999 *   110.53542 *   103.07497 *    13  11  10
   37      37  H            1.09003 *   110.90733 *   158.26023 *    14  13  11
   38      38  H            1.09001 *   110.90339 *   281.82519 *    14  13  11
   39      39  H            1.09005 *   110.53759 *   201.32956 *    15  14  13
   40      40  H            1.09001 *   110.53566 *    78.55506 *    15  14  13
   41      41  H            0.96993 *   119.99500 *   118.88220 *    17  16  15
   42      42  H            1.08997 *   109.46698 *   290.02738 *    22  20  18
   43      43  H            1.09007 *   109.46046 *    49.99265 *    22  20  18
   44      44  H            1.09001 *   109.67100 *   169.71286 *    23  22  20
   45      45  H            1.08994 *   109.67433 *   290.25206 *    23  22  20
   46      46  H            1.08001 *   125.62180 *   179.97438 *    26  25  24
   47      47  H            1.09004 *   109.68460 *   103.26160 *    30  29  24
   48      48  H            1.09003 *   109.68856 *   223.82242 *    30  29  24
   49      49  H            1.08999 *   110.07530 *    29.46915 *    31   9   7
   50      50  H            1.09001 *   110.07691 *   267.92647 *    31   9   7

 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.3162    0.0000    0.0000
     3        14        2.2502    1.6178    0.0000
     4         1        2.3924    2.1072   -1.3948
     5         1        1.5042    2.6200    0.8027
     6         1        3.5963    1.4121    0.5921
     7         6        2.0103   -1.2022   -0.0005
     8         1        3.0804   -1.2079    0.1450
     9         7        1.3381   -2.3808   -0.1855
    10         6        1.8531   -3.7365    0.0677
    11         6        1.1648   -4.6593   -0.9518
    12         1        0.8897   -5.6161   -0.5081
    13         6        2.0721   -4.8428   -2.1865
    14         6        1.0924   -4.7603   -3.3855
    15         6        0.1373   -3.6202   -2.9472
    16         6       -0.0777   -3.8741   -1.4403
    17         7       -1.3195   -4.6090   -1.1867
    18         6       -2.5036   -4.0758   -1.5472
    19         8       -2.5413   -2.9875   -2.0811
    20         6       -3.7810   -4.8320   -1.2868
    21         1       -3.7514   -5.7868   -1.8118
    22         6       -4.9718   -4.0103   -1.7879
    23         6       -6.2397   -4.8724   -1.7550
    24         6       -6.3447   -5.5138   -0.3936
    25         7       -7.4155   -6.0803    0.1689
    26         6       -7.0944   -6.5175    1.3586
    27         7       -5.7892   -6.2420    1.5859
    28         1       -5.2850   -6.4570    2.3862
    29         6       -5.3187   -5.6036    0.4796
    30         6       -3.9299   -5.0810    0.2195
    31         6       -0.0405   -2.5397   -0.6768
    32         1       -0.4851    0.8400    0.0004
    33         1        1.6041   -4.0479    1.0821
    34         1        2.9333   -3.7592   -0.0767
    35         1        2.8107   -4.0431   -2.2406
    36         1        2.5622   -5.8160   -2.1596
    37         1        1.6202   -4.4949   -4.3015
    38         1        0.5493   -5.6972   -3.5092
    39         1       -0.8087   -3.6859   -3.4848
    40         1        0.6031   -2.6480   -3.1084
    41         1       -1.2893   -5.4792   -0.7594
    42         1       -5.1079   -3.1396   -1.1464
    43         1       -4.7822   -3.6821   -2.8099
    44         1       -7.1138   -4.2463   -1.9342
    45         1       -6.1779   -5.6452   -2.5211
    46         1       -7.7600   -7.0197    2.0450
    47         1       -3.7812   -4.1475    0.7623
    48         1       -3.1939   -5.8165    0.5446
    49         1       -0.2960   -1.7184   -1.3463
    50         1       -0.7357   -2.5692    0.1622

          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=WATER                                                                  
 ZINC001099111662.mol2 50                                                        
                                                                                 

     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 03:01:23                  

 Heat of formation + Delta-G solvation     =      11.087590 kcal
 Electronic energy + Delta-G solvation     =  -31517.846151 eV
 Core-core repulsion                       =   27544.510305 eV
 Total energy + Delta-G solvation          =   -3973.335845 eV
 No. of doubly occupied orbitals           =      65
 Molecular weight
  (most abundant/longest-lived isotopes)   =     344.195     amu
 Computer time                             =       1.00     seconds

               Orbital eigenvalues (eV)

 -43.20854 -41.61739 -39.84020 -37.09030 -36.17185 -34.75786 -34.02129 -33.48213
 -31.22997 -30.60237 -30.19761 -29.30001 -26.31682 -26.28786 -25.00364 -23.90597
 -22.53192 -22.36549 -22.10794 -20.71880 -20.24498 -19.76372 -18.85941 -18.76597
 -17.87585 -17.32811 -16.51181 -16.45961 -16.15870 -15.76962 -15.51525 -15.36040
 -15.04733 -14.95540 -14.84830 -14.71245 -14.59432 -14.27323 -14.19078 -13.66110
 -13.60127 -13.00149 -12.92087 -12.85837 -12.63968 -12.57089 -12.35138 -12.10645
 -12.03468 -11.97143 -11.92748 -11.89740 -11.64602 -11.59038 -11.52856 -11.41233
 -10.87025 -10.77951 -10.68099 -10.55713  -9.91808  -9.10431  -8.85668  -7.88299
  -5.25925   1.09398   1.41650   1.44505   1.53334   1.72298   2.14919   2.21701
   2.32105   2.48489   3.30112   3.44243   3.50564   3.66439   3.76824   3.79337
   3.82504   3.90976   3.94074   4.17185   4.23730   4.33272   4.43763   4.45663
   4.46748   4.53743   4.55433   4.63553   4.66348   4.72398   4.77886   4.81225
   4.93980   5.02129   5.04698   5.07116   5.10885   5.13421   5.14829   5.28937
   5.32139   5.33162   5.35716   5.44143   5.47891   5.59031   5.72578   5.79918
   5.90055   6.22095   6.60780   6.74174   7.04253   7.15842   7.46999   7.75664
   8.45548   9.22553

          Molecular weight =  344.20amu

           Principal moments of inertia in cm(-1)
          A =    0.012838   B =    0.003182   C =    0.002896

           Principal moments of inertia in units of 10**(-40)*gram-cm**2
          A = 2180.506507   B = 8796.013733   C = 9667.819360

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

         Atom NO.   Type          Charge        No. of electrons
           1         N            -0.945           5.945
           2         C            -0.066           4.066
           3         Si            0.974           3.026
           4         H            -0.290           1.290
           5         H            -0.297           1.297
           6         H            -0.260           1.260
           7         C            -0.240           4.240
           8         H             0.084           0.916
           9         N            -0.617           5.617
          10         C             0.128           3.872
          11         C            -0.093           4.093
          12         H             0.119           0.881
          13         C            -0.105           4.105
          14         C            -0.110           4.110
          15         C            -0.110           4.110
          16         C             0.157           3.843
          17         N            -0.698           5.698
          18         C             0.511           3.489
          19         O            -0.570           6.570
          20         C            -0.083           4.083
          21         H             0.131           0.869
          22         C            -0.094           4.094
          23         C            -0.026           4.026
          24         C            -0.020           4.020
          25         N            -0.530           5.530
          26         C             0.187           3.813
          27         N            -0.549           5.549
          28         H             0.437           0.563
          29         C            -0.028           4.028
          30         C            -0.027           4.027
          31         C             0.197           3.803
          32         H             0.243           0.757
          33         H             0.016           0.984
          34         H             0.055           0.945
          35         H             0.049           0.951
          36         H             0.099           0.901
          37         H             0.070           0.930
          38         H             0.089           0.911
          39         H             0.056           0.944
          40         H             0.032           0.968
          41         H             0.413           0.587
          42         H             0.045           0.955
          43         H             0.075           0.925
          44         H             0.075           0.925
          45         H             0.090           0.910
          46         H             0.247           0.753
          47         H             0.062           0.938
          48         H             0.109           0.891
          49         H             0.030           0.970
          50         H            -0.025           1.025

For the charges calculated by CM2:

 Dipole moment (debyes)    X         Y         Z       Total
 from point charges    -7.072   -21.531     0.266    22.664


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.592           5.592
           2         C            -0.259           4.259
           3         Si            0.802           3.198
           4         H            -0.217           1.217
           5         H            -0.224           1.224
           6         H            -0.184           1.184
           7         C            -0.369           4.369
           8         H             0.102           0.898
           9         N            -0.342           5.342
          10         C             0.003           3.997
          11         C            -0.113           4.113
          12         H             0.137           0.863
          13         C            -0.143           4.143
          14         C            -0.147           4.147
          15         C            -0.148           4.148
          16         C             0.069           3.931
          17         N            -0.350           5.350
          18         C             0.298           3.702
          19         O            -0.448           6.448
          20         C            -0.105           4.105
          21         H             0.148           0.852
          22         C            -0.133           4.133
          23         C            -0.064           4.064
          24         C            -0.149           4.149
          25         N            -0.250           5.250
          26         C            -0.097           4.097
          27         N            -0.153           5.153
          28         H             0.278           0.722
          29         C            -0.145           4.145
          30         C            -0.065           4.065
          31         C             0.070           3.930
          32         H             0.053           0.947
          33         H             0.034           0.966
          34         H             0.073           0.927
          35         H             0.068           0.932
          36         H             0.118           0.882
          37         H             0.089           0.911
          38         H             0.108           0.892
          39         H             0.075           0.925
          40         H             0.051           0.949
          41         H             0.251           0.749
          42         H             0.064           0.936
          43         H             0.093           0.907
          44         H             0.093           0.907
          45         H             0.108           0.892
          46         H             0.263           0.737
          47         H             0.081           0.919
          48         H             0.127           0.873
          49         H             0.048           0.952
          50         H            -0.007           1.007

 Dipole moment (debyes)     X         Y         Z       Total
 from point charges     -7.893   -20.598    -0.452    22.063
 hybrid contribution     0.952     1.525     1.294     2.215
 sum                    -6.941   -19.073     0.842    20.314

          Atomic orbital electron populations

   1.70211   1.06622   1.25181   1.57216   1.25757   0.96320   1.01415   1.02431
   0.85123   0.78142   0.79309   0.77186   1.21715   1.22370   1.18440   1.19396
   0.91367   0.83777   1.42348   0.89815   1.44703   1.12457   1.01445   1.75547
   1.21639   0.97587   0.84695   0.95774   1.22828   0.93424   0.98274   0.96782
   0.86299   1.22413   0.96855   1.02392   0.92634   1.22557   0.96434   0.97787
   0.97948   1.22495   0.99416   0.98523   0.94397   1.21674   0.84077   0.92322
   0.95016   1.45811   1.04618   1.22510   1.62097   1.21052   0.85807   0.85801
   0.77567   1.90769   1.87574   1.22592   1.43899   1.21242   0.91046   1.00914
   0.97298   0.85160   1.21337   0.92659   0.97230   1.02066   1.19896   0.97254
   0.97970   0.91288   1.20777   0.89783   1.05505   0.98843   1.69943   1.30116
   1.16485   1.08411   1.25805   0.84309   1.03530   0.96016   1.42903   1.11753
   1.44425   1.16193   0.72239   1.20797   0.97816   1.05691   0.90184   1.19926
   0.91085   1.01153   0.94364   1.20943   0.85010   0.93749   0.93308   0.94748
   0.96588   0.92654   0.93204   0.88208   0.91143   0.89207   0.92508   0.94864
   0.74909   0.93575   0.90676   0.90683   0.89165   0.73710   0.91923   0.87278
   0.95183   1.00715

   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       N    -0.94   -60.84    10.84      21.65       0.23   -60.61     16
   2       C    -0.07    -3.89     5.39      84.86       0.46    -3.44     16
   3       Si    0.97    47.55    29.56      68.60       2.03    49.58     16
   4       H    -0.29   -14.10     7.11      99.48       0.71   -13.39     16
   5       H    -0.30   -14.64     7.11      99.48       0.71   -13.93     16
   6       H    -0.26   -11.70     7.11      99.48       0.71   -10.99     16
   7       C    -0.24   -13.73    10.28      84.03       0.86   -12.87     16
   8       H     0.08     4.53     7.84      -2.91      -0.02     4.50     16
   9       N    -0.62   -32.10     3.39    -695.36      -2.35   -34.45     16
  10       C     0.13     4.99     6.60      86.51       0.57     5.57     16
  11       C    -0.09    -2.74     3.63      -9.44      -0.03    -2.78     16
  12       H     0.12     2.47     7.71      -2.39      -0.02     2.45     16
  13       C    -0.11    -2.70     6.25      31.74       0.20    -2.50     16
  14       C    -0.11    -2.81     6.92      32.00       0.22    -2.59     16
  15       C    -0.11    -4.02     5.53      31.74       0.18    -3.84     16
  16       C     0.16     5.70     0.91       4.96       0.00     5.70     16
  17       N    -0.70   -19.78     5.09    -416.81      -2.12   -21.90     16
  18       C     0.51    16.33     6.43      87.66       0.56    16.89     16
  19       O    -0.57   -25.11    13.79      -3.03      -0.04   -25.15     16
  20       C    -0.08    -1.89     2.51     -11.37      -0.03    -1.92     16
  21       H     0.13     2.16     8.12      -2.39      -0.02     2.14     16
  22       C    -0.09    -2.43     4.64      30.75       0.14    -2.29     16
  23       C    -0.03    -0.58     6.30      30.02       0.19    -0.39     16
  24       C    -0.02    -0.48     6.98      41.61       0.29    -0.19     16
  25       N    -0.53   -12.88    11.25    -187.90      -2.11   -14.99     16
  26       C     0.19     2.95    13.05     180.10       2.35     5.30     16
  27       N    -0.55    -7.85     6.15    -186.45      -1.15    -9.00     16
  28       H     0.44     3.12     8.96     -92.70      -0.83     2.29     16
  29       C    -0.03    -0.56     7.08      41.14       0.29    -0.27     16
  30       C    -0.03    -0.53     5.43      29.97       0.16    -0.37     16
  31       C     0.20    10.17     3.70      86.51       0.32    10.49     16
  32       H     0.24    15.44     8.31     -92.71      -0.77    14.67     16
  33       H     0.02     0.63     8.14      -2.39      -0.02     0.61     16
  34       H     0.06     2.03     7.63      -2.38      -0.02     2.02     16
  35       H     0.05     1.54     7.61      -2.39      -0.02     1.52     16
  36       H     0.10     1.80     8.14      -2.39      -0.02     1.78     16
  37       H     0.07     1.70     8.14      -2.39      -0.02     1.68     16
  38       H     0.09     1.81     8.14      -2.39      -0.02     1.79     16
  39       H     0.06     2.14     6.65      -2.38      -0.02     2.13     16
  40       H     0.03     1.50     7.61      -2.39      -0.02     1.48     16
  41       H     0.41     7.81     7.97     -92.71      -0.74     7.07     16
  42       H     0.05     1.42     8.11      -2.39      -0.02     1.40     16
  43       H     0.07     2.05     7.77      -2.38      -0.02     2.03     16
  44       H     0.07     1.73     8.14      -2.39      -0.02     1.71     16
  45       H     0.09     1.73     8.14      -2.39      -0.02     1.71     16
  46       H     0.25     1.51     8.06      -2.91      -0.02     1.49     16
  47       H     0.06     1.60     8.11      -2.39      -0.02     1.58     16
  48       H     0.11     1.47     7.90      -2.39      -0.02     1.45     16
  49       H     0.03     1.87     6.22      -2.39      -0.01     1.86     16
  50       H    -0.03    -1.46     8.06      -2.39      -0.02    -1.48     16

 Total:         -1.00   -87.08   384.52                  0.62   -86.46

 By element:

 Atomic #  1  Polarization:   20.16  SS G_CDS:  -0.60  Total:    19.56 kcal
 Atomic #  6  Polarization:    3.77  SS G_CDS:   6.74  Total:    10.51 kcal
 Atomic #  7  Polarization: -133.45  SS G_CDS:  -7.50  Total:  -140.95 kcal
 Atomic #  8  Polarization:  -25.11  SS G_CDS:  -0.04  Total:   -25.15 kcal
 Atomic # 14  Polarization:   47.55  SS G_CDS:   2.03  Total:    49.58 kcal

 Total:                      -87.08              0.62           -86.46 kcal

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

                                                                                 
 ZINC001099111662.mol2 50                                                        

 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                 97.547 kcal
 (2)  G-P(sol) polarization free energy of solvation               -87.083 kcal
 (3)  G-ENP(sol) elect.-nuc.-pol. free energy of system             10.463 kcal
 (4)  G-CDS(sol) cavity-dispersion-solvent structure free energy     0.624 kcal
 (5)  G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4)             -86.459 kcal
 (6)  G-S(sol) free energy of system = (1) + (5)                    11.088 kcal

     1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED

 Total computer time =       1.01 seconds