reb01

这是程序是老师提供的源程序

               COMPUTATION IN CARTISIAN COORDINATES
               ****************************************


 URELAX= .40 VRELAX= .40 TRELAX= .80 PRELAX= .70 RA=.100E+05   PR= .70

 NTIMES:  2  2  2  2
 LBLK(I):   T   T   T   T
   L1=42  M1=42  I1= 8  I2=14  J1= 8  J2=28  a=67.00  b= 7.00  d= 6.50
   af=  .229E-04  SRAPR=   .8373E+02  TM= 1.00  DT= .2000E+02
   ITER      SMAX        SSUM     V(15,18)    T(15,18)    P(15,18)    NUm
      0  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
      1  0.0000E+00  0.0000E+00  0.0000E+00  3.1502E-01  0.0000E+00 -4.2694E+01
      2  1.4727E-04  2.7915E-10  5.2395E-03  5.1937E-01 -1.9871E-03 -2.3793E+01
      3  8.2879E-05  5.0517E-10  1.3211E-02  6.3792E-01 -1.2787E-03 -1.6615E+01
      4  5.8208E-05 -5.4395E-10  2.2186E-02  7.0994E-01  1.6415E-03 -1.3215E+01
      5  4.2728E-05  1.2118E-09  3.1344E-02  7.5698E-01 -7.7168E-03 -1.1302E+01
      6  4.5708E-05 -3.0556E-10  4.0338E-02  7.8978E-01 -1.9533E-02 -1.0092E+01
      7  4.7783E-05  3.1793E-10  4.9034E-02  8.1388E-01 -3.2927E-02 -9.2657E+00
      8  4.6326E-05  1.0475E-09  5.7385E-02  8.3230E-01 -1.0438E-01 -8.6719E+00
      9  4.1235E-05  1.7224E-10  6.5390E-02  8.4681E-01 -1.0831E-01 -8.2284E+00
     10  3.6786E-05  7.8356E-10  7.3064E-02  8.5850E-01 -1.2466E-01 -7.8870E+00
     11  3.3540E-05 -1.0560E-09  8.0435E-02  8.6808E-01 -1.2896E-01 -7.6175E+00
     12  3.0626E-05  1.3688E-09  8.7537E-02  8.7601E-01 -1.3169E-01 -7.4001E+00
     13  2.8903E-05  2.6299E-09  9.4407E-02  8.8261E-01 -1.3629E-01 -7.2217E+00
     14  2.9634E-05  2.3184E-10  1.0108E-01  8.8812E-01 -1.4401E-01 -7.0729E+00
     15  3.0270E-05  8.3451E-10  1.0759E-01  8.9270E-01 -1.4827E-01 -6.9474E+00
     16  3.0698E-05  1.5347E-09  1.1396E-01  8.9649E-01 -1.5411E-01 -6.8405E+00
     17  3.1533E-05 -6.3581E-09  1.2021E-01  8.9957E-01 -1.5330E-01 -6.7488E+00
     18  3.4189E-05  3.1383E-09  1.2636E-01  9.0203E-01 -1.4843E-01 -6.6698E+00
     19  3.6761E-05 -7.4465E-10  1.3242E-01  9.0392E-01 -1.4790E-01 -6.6015E+00
     20  3.9646E-05 -1.7324E-09  1.3840E-01  9.0531E-01 -1.3881E-01 -6.5425E+00
     21  4.1597E-05  1.0002E-09  1.4430E-01  9.0625E-01 -1.3558E-01 -6.4917E+00
     22  4.2487E-05  4.9751E-09  1.5012E-01  9.0677E-01 -1.4276E-01 -6.4481E+00
     23  4.3174E-05  8.5351E-09  1.5587E-01  9.0691E-01 -1.2294E-01 -6.4112E+00
     24  4.5403E-05 -1.2918E-09  1.6155E-01  9.0672E-01 -1.2580E-01 -6.3802E+00
     25  4.7098E-05 -6.4193E-09  1.6716E-01  9.0623E-01 -1.8607E-01 -6.3547E+00
     26  4.9572E-05  2.9949E-10  1.7271E-01  9.0547E-01 -2.3886E-01 -6.3342E+00
     27  5.1835E-05 -1.3660E-08  1.7820E-01  9.0449E-01 -2.3673E-01 -6.3184E+00
     28  5.3890E-05  1.5596E-08  1.8363E-01  9.0331E-01 -2.7690E-01 -6.3068E+00
     29  5.5679E-05 -4.4892E-09  1.8900E-01  9.0196E-01 -2.8787E-01 -6.2993E+00
     30  5.7376E-05  7.0498E-09  1.9430E-01  9.0047E-01 -2.8047E-01 -6.2954E+00
     31  5.9003E-05 -8.6760E-09  1.9954E-01  8.9887E-01 -2.6787E-01 -6.2948E+00
     32  6.0139E-05  5.7087E-09  2.0471E-01  8.9718E-01 -2.9476E-01 -6.2972E+00
     33  6.2762E-05  1.5707E-09  2.0980E-01  8.9542E-01 -3.0447E-01 -6.3023E+00
     34  6.4828E-05  1.4626E-09  2.1480E-01  8.9362E-01 -3.2195E-01 -6.3098E+00
     35  6.8582E-05 -8.3530E-09  2.1970E-01  8.9178E-01 -2.9544E-01 -6.3194E+00
     36  7.6034E-05 -2.7520E-09  2.2449E-01  8.8993E-01 -2.8323E-01 -6.3307E+00
     37  8.2287E-05  1.7001E-08  2.2916E-01  8.8808E-01 -3.5518E-01 -6.3436E+00
     38  8.6943E-05 -9.4134E-09  2.3369E-01  8.8624E-01 -3.3197E-01 -6.3577E+00
     39  8.9708E-05 -1.9493E-08  2.3807E-01  8.8443E-01 -3.8280E-01 -6.3728E+00
     40  9.0396E-05  1.7721E-08  2.4230E-01  8.8265E-01 -3.8097E-01 -6.3886E+00
     41  8.9842E-05  5.0176E-09  2.4635E-01  8.8093E-01 -3.6281E-01 -6.4049E+00
     42  9.1214E-05 -5.4450E-09  2.5022E-01  8.7925E-01 -3.7391E-01 -6.4214E+00
     43  9.2751E-05  4.2723E-10  2.5389E-01  8.7765E-01 -3.8030E-01 -6.4381E+00
     44  9.9100E-05 -1.0635E-08  2.5736E-01  8.7611E-01 -3.6823E-01 -6.4546E+00
     45  1.0431E-04 -2.0631E-08  2.6062E-01  8.7465E-01 -3.5188E-01 -6.4709E+00
     46  1.0818E-04 -9.8783E-09  2.6366E-01  8.7327E-01 -3.4290E-01 -6.4868E+00
     47  1.1060E-04 -6.2621E-09  2.6648E-01  8.7198E-01 -3.2535E-01 -6.5023E+00
     48  1.1152E-04  2.1719E-09  2.6908E-01  8.7078E-01 -3.5311E-01 -6.5171E+00
     49  1.1098E-04 -1.1904E-08  2.7144E-01  8.6968E-01 -3.4335E-01 -6.5312E+00
     50  1.0908E-04  9.4658E-09  2.7358E-01  8.6867E-01 -3.7044E-01 -6.5445E+00
     51  1.0597E-04 -1.2069E-08  2.7550E-01  8.6776E-01 -3.8257E-01 -6.5569E+00
     52  1.0185E-04  3.5186E-09  2.7719E-01  8.6695E-01 -3.9090E-01 -6.5685E+00
     53  9.6922E-05  2.2271E-08  2.7867E-01  8.6624E-01 -4.0186E-01 -6.5791E+00
     54  9.1416E-05 -3.1414E-09  2.7993E-01  8.6563E-01 -4.3246E-01 -6.5887E+00
     55  8.5538E-05 -7.8991E-09  2.8099E-01  8.6511E-01 -4.3514E-01 -6.5974E+00
     56  7.9485E-05  1.7986E-08  2.8184E-01  8.6469E-01 -4.1796E-01 -6.6050E+00
     57  7.3417E-05  7.2069E-09  2.8251E-01  8.6437E-01 -4.3114E-01 -6.6115E+00
     58  6.7477E-05  9.0938E-10  2.8299E-01  8.6414E-01 -4.3000E-01 -6.6171E+00
     59  6.1770E-05 -3.3812E-09  2.8330E-01  8.6399E-01 -4.3766E-01 -6.6216E+00
     60  5.6375E-05 -1.3620E-08  2.8343E-01  8.6394E-01 -4.0900E-01 -6.6251E+00
     61  5.1348E-05  2.0638E-08  2.8342E-01  8.6396E-01 -6.5968E-01 -6.6276E+00
     62  4.6712E-05  6.2297E-09  2.8325E-01  8.6406E-01 -6.5847E-01 -6.6292E+00
     63  4.2498E-05  4.8229E-09  2.8295E-01  8.6423E-01 -6.8297E-01 -6.6300E+00
     64  3.8689E-05 -2.2738E-08  2.8252E-01  8.6446E-01 -5.0092E-01 -6.6298E+00
     65  3.5280E-05  7.9378E-09  2.8198E-01  8.6476E-01 -5.2162E-01 -6.6288E+00
     66  3.2242E-05 -1.8860E-09  2.8133E-01  8.6512E-01 -5.2128E-01 -6.6271E+00
     67  2.9560E-05  8.7107E-10  2.8058E-01  8.6553E-01 -5.2102E-01 -6.6246E+00
     68  2.7202E-05 -3.2666E-08  2.7974E-01  8.6599E-01 -4.9600E-01 -6.6215E+00
     69  2.5140E-05 -3.6681E-09  2.7883E-01  8.6649E-01 -4.8914E-01 -6.6178E+00
     70  2.3340E-05 -2.4539E-09  2.7784E-01  8.6702E-01 -4.8414E-01 -6.6135E+00
     71  2.1779E-05  3.0707E-10  2.7680E-01  8.6759E-01 -4.7503E-01 -6.6087E+00
     72  2.0425E-05  1.2030E-08  2.7571E-01  8.6819E-01 -4.4161E-01 -6.6034E+00
     73  1.9255E-05  1.1279E-08  2.7458E-01  8.6881E-01 -4.5486E-01 -6.5977E+00
     74  1.8246E-05 -2.4200E-08  2.7341E-01  8.6945E-01 -4.3290E-01 -6.5917E+00
     75  1.7376E-05  1.3782E-08  2.7222E-01  8.7011E-01 -4.5152E-01 -6.5853E+00
     76  1.6627E-05 -3.4306E-09  2.7100E-01  8.7078E-01 -4.4796E-01 -6.5786E+00
     77  1.5979E-05 -1.6100E-08  2.6978E-01  8.7146E-01 -6.5070E-01 -6.5717E+00
     78  1.5418E-05 -6.9747E-11  2.6854E-01  8.7214E-01 -6.5215E-01 -6.5646E+00
     79  1.4930E-05 -1.2761E-08  2.6730E-01  8.7282E-01 -6.3692E-01 -6.5573E+00
     80  1.4506E-05  9.0515E-09  2.6607E-01  8.7351E-01 -3.6523E-01 -6.5498E+00
     81  1.4130E-05  2.8343E-08  2.6484E-01  8.7419E-01 -3.8317E-01 -6.5423E+00
     82  1.3803E-05 -2.4257E-08  2.6363E-01  8.7487E-01 -3.7038E-01 -6.5347E+00
     83  1.3505E-05  4.8748E-09  2.6242E-01  8.7555E-01 -3.2554E-01 -6.5270E+00
     84  1.3237E-05  1.5443E-08  2.6124E-01  8.7621E-01 -3.4247E-01 -6.5193E+00
     85  1.2987E-05 -4.6893E-09  2.6007E-01  8.7687E-01 -3.4113E-01 -6.5115E+00
     86  1.2755E-05 -3.6506E-09  2.5893E-01  8.7751E-01 -3.4031E-01 -6.5038E+00
     87  1.2535E-05 -5.3160E-10  2.5782E-01  8.7815E-01 -3.4493E-01 -6.4961E+00
     88  1.2323E-05  6.8572E-09  2.5672E-01  8.7877E-01 -1.6923E+00 -6.4885E+00
     89  1.2120E-05  2.2747E-08  2.5566E-01  8.7937E-01 -1.6548E+00 -6.4809E+00
     90  1.1922E-05  1.2347E-08  2.5463E-01  8.7997E-01 -1.6534E+00 -6.4734E+00
     91  1.1719E-05 -6.0048E-09  2.5362E-01  8.8054E-01 -1.6569E+00 -6.4659E+00
     92  1.1519E-05 -1.8247E-09  2.5265E-01  8.8111E-01 -1.6554E+00 -6.4586E+00
     93  1.1319E-05  1.3505E-08  2.5171E-01  8.8166E-01 -1.6566E+00 -6.4513E+00
     94  1.1199E-05  2.1827E-08  2.5079E-01  8.8219E-01 -1.7113E+00 -6.4442E+00
     95  1.1124E-05 -1.9849E-08  2.4991E-01  8.8271E-01 -1.7289E+00 -6.4371E+00
     96  1.1050E-05  5.2338E-09  2.4907E-01  8.8321E-01 -1.7241E+00 -6.4302E+00
     97  1.0951E-05 -6.7221E-09  2.4825E-01  8.8370E-01 -1.7076E+00 -6.4234E+00
     98  1.0858E-05 -7.9955E-10  2.4746E-01  8.8418E-01 -1.7085E+00 -6.4167E+00
     99  1.0729E-05  1.6601E-08  2.4671E-01  8.8464E-01 -1.7254E+00 -6.4101E+00
    100  1.0632E-05 -1.8176E-08  2.4598E-01  8.8509E-01 -1.7115E+00 -6.4036E+00
    101  1.0512E-05 -1.7131E-08  2.4528E-01  8.8552E-01 -1.7037E+00 -6.3972E+00
    102  1.0365E-05  1.1939E-08  2.4461E-01  8.8594E-01 -1.6766E+00 -6.3910E+00
    103  1.0241E-05 -1.4549E-08  2.4397E-01  8.8634E-01 -1.6741E+00 -6.3848E+00
    104  1.0101E-05  4.0638E-09  2.4336E-01  8.8674E-01 -1.6633E+00 -6.3788E+00
    105  9.9664E-06 -6.8502E-09  2.4277E-01  8.8712E-01 -1.6709E+00 -6.3729E+00
    106  9.8341E-06 -5.8205E-10  2.4221E-01  8.8749E-01 -1.6707E+00 -6.3671E+00
    107  9.6994E-06  2.6659E-09  2.4167E-01  8.8785E-01 -1.6724E+00 -6.3615E+00
    108  9.5489E-06  2.4707E-08  2.4116E-01  8.8819E-01 -1.7301E+00 -6.3559E+00
    109  9.4146E-06 -9.4376E-09  2.4066E-01  8.8853E-01 -1.7216E+00 -6.3505E+00
    110  9.2936E-06  5.1215E-09  2.4019E-01  8.8885E-01 -1.7329E+00 -6.3451E+00
    111  9.1552E-06 -3.7129E-09  2.3975E-01  8.8917E-01 -1.7286E+00 -6.3399E+00
    112  9.0157E-06  1.2691E-08  2.3932E-01  8.8948E-01 -1.7236E+00 -6.3347E+00
    113  8.8734E-06  9.8048E-10  2.3891E-01  8.8977E-01 -1.7256E+00 -6.3297E+00
    114  8.7285E-06  7.3322E-09  2.3852E-01  8.9006E-01 -1.7751E+00 -6.3248E+00
    115  8.5827E-06 -1.0677E-08  2.3815E-01  8.9034E-01 -1.7094E+00 -6.3200E+00
    116  8.4380E-06 -3.1667E-08  2.3779E-01  8.9061E-01 -1.6705E+00 -6.3153E+00
    117  8.3051E-06  4.1829E-09  2.3746E-01  8.9087E-01 -1.6633E+00 -6.3106E+00
    118  8.1556E-06  8.1225E-09  2.3714E-01  8.9113E-01 -1.6334E+00 -6.3061E+00
    119  8.0245E-06  2.3059E-08  2.3683E-01  8.9137E-01 -1.6165E+00 -6.3017E+00
    120  7.8733E-06 -2.2792E-08  2.3654E-01  8.9161E-01 -1.7539E+00 -6.2973E+00
    121  7.7289E-06 -2.2026E-08  2.3626E-01  8.9185E-01 -1.7683E+00 -6.2931E+00
    122  7.5972E-06 -8.6284E-09  2.3599E-01  8.9207E-01 -1.7741E+00 -6.2889E+00
    123  7.4554E-06 -1.6066E-08  2.3574E-01  8.9229E-01 -1.6498E+00 -6.2848E+00
    124  7.3120E-06  6.4620E-09  2.3550E-01  8.9250E-01 -1.6563E+00 -6.2808E+00
    125  7.1786E-06  1.0425E-09  2.3528E-01  8.9271E-01 -1.6574E+00 -6.2769E+00
    126  7.0305E-06  1.0285E-10  2.3506E-01  8.9291E-01 -1.6593E+00 -6.2731E+00
    127  6.9066E-06  1.4116E-08  2.3485E-01  8.9310E-01 -1.6714E+00 -6.2694E+00
    128  6.7703E-06 -2.4885E-09  2.3466E-01  8.9329E-01 -1.6665E+00 -6.2657E+00
    129  6.6425E-06  3.2614E-09  2.3447E-01  8.9348E-01 -1.6690E+00 -6.2622E+00
    130  6.5075E-06 -3.3781E-09  2.3429E-01  8.9366E-01 -1.6733E+00 -6.2587E+00
    131  6.3780E-06  3.6186E-09  2.3413E-01  8.9383E-01 -1.6798E+00 -6.2553E+00
    132  6.2631E-06 -1.1880E-10  2.3397E-01  8.9400E-01 -1.6804E+00 -6.2519E+00
    133  6.1402E-06 -1.6996E-08  2.3382E-01  8.9416E-01 -1.6119E+00 -6.2487E+00
    134  5.9960E-06  1.3992E-08  2.3367E-01  8.9432E-01 -1.6406E+00 -6.2455E+00
    135  5.8781E-06  2.4411E-08  2.3354E-01  8.9447E-01 -1.6624E+00 -6.2424E+00
    136  5.7596E-06 -2.1919E-09  2.3341E-01  8.9462E-01 -1.6597E+00 -6.2393E+00
    137  5.6241E-06 -1.8455E-08  2.3329E-01  8.9477E-01 -1.6379E+00 -6.2363E+00
    138  5.5026E-06  7.6432E-09  2.3317E-01  8.9491E-01 -1.6394E+00 -6.2334E+00
    139  5.3979E-06 -1.4102E-08  2.3307E-01  8.9505E-01 -1.6082E+00 -6.2306E+00
    140  5.2615E-06  1.5043E-08  2.3296E-01  8.9518E-01 -1.6200E+00 -6.2278E+00
    141  5.1429E-06 -7.9690E-09  2.3287E-01  8.9531E-01 -1.6248E+00 -6.2251E+00
    142  5.0455E-06  1.1987E-08  2.3277E-01  8.9544E-01 -1.5893E+00 -6.2225E+00
    143  4.9189E-06  7.2677E-09  2.3269E-01  8.9556E-01 -1.5943E+00 -6.2199E+00
    144  4.8133E-06 -1.1056E-08  2.3261E-01  8.9568E-01 -1.4321E+00 -6.2174E+00
    145  4.7114E-06 -6.3381E-09  2.3253E-01  8.9580E-01 -1.4192E+00 -6.2149E+00
    146  4.6656E-06 -1.5263E-08  2.3246E-01  8.9591E-01 -1.3358E+00 -6.2125E+00
    147  4.6062E-06 -3.0061E-09  2.3239E-01  8.9602E-01 -1.3156E+00 -6.2102E+00
    148  4.5542E-06 -2.3605E-09  2.3232E-01  8.9613E-01 -1.3088E+00 -6.2079E+00
    149  4.4932E-06  1.0096E-09  2.3226E-01  8.9623E-01 -1.3097E+00 -6.2057E+00
    150  4.4328E-06 -1.9835E-08  2.3221E-01  8.9633E-01 -1.3324E+00 -6.2035E+00