# input opsmode afspc a, improved i i
# input type of run c, v, m: v
# input constants 721, 72, 84 72
# input elset filename: SGP4-VER.TLE

# Must set globals same as testmat.m's run

global opsmode = "i"; 
global tumin  mu  radiusearthkm  xke  j2  j3  j4  j3oj2

[tumin, mu, radiusearthkm, xke, j2, j3, j4, j3oj2] = getgravc(72);

TOL = 0.000001;

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 1.859667000000e-01,  1.844178495062e+04,  5.980929187319e-01,  4.722944544078e-02, 5);
assert(isequalRel(ainv,          7.386070850905e-01, TOL));
assert(isequalRel(ao,            1.353899820603e+00, TOL));
assert(isequalRel(con41,         1.048865087996e+00, TOL));
assert(isequalRel(con42,        -2.414775146659e+00, TOL));
assert(isequalRel(cosio,         8.264109324857e-01, TOL));
assert(isequalRel(cosio2,        6.829550293319e-01, TOL));
assert(isequalRel(einv,          5.377306797400e+00, TOL));
assert(isequalRel(eccsq,         3.458361350889e-02, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        9.654163864911e-01, TOL));
assert(isequalRel(posq,          1.708450473396e+00, TOL));
assert(isequalRel(rp,            1.102119538835e+00, TOL));
assert(isequalRel(rteosq,        9.825560475063e-01, TOL));
assert(isequalRel(sinio,         5.630674654676e-01, TOL));
assert(isequalRel(gsto,          3.469172342379e+00, TOL));
assert(isequalRel(no,            4.720630155918e-02, TOL));
# 5

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 1.450506000000e-01,  1.975491070959e+04,  2.000636014976e-01,  5.246109836613e-03, 4632);
assert(isequalRel(ainv,          1.707228925380e-01, TOL));
assert(isequalRel(ao,            5.857445273648e+00, TOL));
assert(isequalRel(con41,         1.881517177058e+00, TOL));
assert(isequalRel(con42,        -3.802528628430e+00, TOL));
assert(isequalRel(cosio,         9.800539401920e-01, TOL));
assert(isequalRel(cosio2,        9.605057256859e-01, TOL));
assert(isequalRel(einv,          6.894145904946e+00, TOL));
assert(isequalRel(eccsq,         2.103967656036e-02, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        9.789603234396e-01, TOL));
assert(isequalRel(posq,          3.288112440965e+01, TOL));
assert(isequalRel(rp,            5.007819322238e+00, TOL));
assert(isequalRel(rteosq,        9.894242383526e-01, TOL));
assert(isequalRel(sinio,         1.987316640953e-01, TOL));
assert(isequalRel(gsto,          1.716027084070e+00, TOL));
assert(isequalRel(no,            5.245868658927e-03, TOL));
# 4632

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 3.003500000000e-03,  2.063082412014e+04,  1.013301511794e+00,  6.791020666036e-02, 6251);
assert(isequalRel(ainv,          9.413191670523e-01, TOL));
assert(isequalRel(ao,            1.062338933490e+00, TOL));
assert(isequalRel(con41,        -1.602801925378e-01, TOL));
assert(isequalRel(con42,        -3.995330124370e-01, TOL));
assert(isequalRel(cosio,         5.290620024982e-01, TOL));
assert(isequalRel(cosio2,        2.799066024874e-01, TOL));
assert(isequalRel(einv,          3.329448976194e+02, TOL));
assert(isequalRel(eccsq,         9.021012250000e-06, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        9.999909789878e-01, TOL));
assert(isequalRel(posq,          1.128543648121e+00, TOL));
assert(isequalRel(rp,            1.059148198503e+00, TOL));
assert(isequalRel(rteosq,        9.999954894837e-01, TOL));
assert(isequalRel(sinio,         8.485831706513e-01, TOL));
assert(isequalRel(gsto,          3.673754967933e+00, TOL));
assert(isequalRel(no,            6.791803699335e-02, TOL));
# 6251

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 6.877146000000e-01,  2.063033215444e+04,  1.119778813470e+00,  8.748086888067e-03, 8195);
assert(isequalRel(ainv,          2.400881758466e-01, TOL));
assert(isequalRel(ao,            4.165136398215e+00, TOL));
assert(isequalRel(con41,        -4.300218866316e-01, TOL));
assert(isequalRel(con42,         5.003647771941e-02, TOL));
assert(isequalRel(cosio,         4.358815257110e-01, TOL));
assert(isequalRel(cosio2,        1.899927044561e-01, TOL));
assert(isequalRel(einv,          1.454091566472e+00, TOL));
assert(isequalRel(eccsq,         4.729513710532e-01, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        5.270486289468e-01, TOL));
assert(isequalRel(posq,          4.819032241803e+00, TOL));
assert(isequalRel(rp,            1.300711286171e+00, TOL));
assert(isequalRel(rteosq,        7.259811491677e-01, TOL));
assert(isequalRel(sinio,         9.000040530708e-01, TOL));
assert(isequalRel(gsto,          5.741801269032e-01, TOL));
assert(isequalRel(no,            8.748547019630e-03, TOL));
# 8195

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 7.069051000000e-01,  2.063056157475e+04,  1.127426846252e+00,  8.762146867822e-03, 9880);
assert(isequalRel(ainv,          2.403464391941e-01, TOL));
assert(isequalRel(ao,            4.160660766820e+00, TOL));
assert(isequalRel(con41,        -4.479141048878e-01, TOL));
assert(isequalRel(con42,         7.985684147962e-02, TOL));
assert(isequalRel(cosio,         4.289855844945e-01, TOL));
assert(isequalRel(cosio2,        1.840286317041e-01, TOL));
assert(isequalRel(einv,          1.414617039826e+00, TOL));
assert(isequalRel(eccsq,         4.997148204060e-01, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        5.002851795940e-01, TOL));
assert(isequalRel(posq,          4.332712683910e+00, TOL));
assert(isequalRel(rp,            1.219468451385e+00, TOL));
assert(isequalRel(rteosq,        7.073084048659e-01, TOL));
assert(isequalRel(sinio,         9.033113351973e-01, TOL));
assert(isequalRel(gsto,          2.019617116981e+00, TOL));
assert(isequalRel(no,            8.762667060363e-03, TOL));
# 9880

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 2.709710000000e-02,  2.023779417928e+04,  1.657329751109e-01,  5.069604863892e-03, 9998);
assert(isequalRel(ainv,          1.668720698242e-01, TOL));
assert(isequalRel(ao,            5.992614588250e+00, TOL));
assert(isequalRel(con41,         1.918349444333e+00, TOL));
assert(isequalRel(con42,        -3.863915740555e+00, TOL));
assert(isequalRel(cosio,         9.862976975087e-01, TOL));
assert(isequalRel(cosio2,        9.727831481111e-01, TOL));
assert(isequalRel(einv,          3.690431817427e+01, TOL));
assert(isequalRel(eccsq,         7.342528284100e-04, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        9.992657471716e-01, TOL));
assert(isequalRel(posq,          3.585871282662e+01, TOL));
assert(isequalRel(rp,            5.830232111491e+00, TOL));
assert(isequalRel(rteosq,        9.996328061701e-01, TOL));
assert(isequalRel(sinio,         1.649753069066e-01, TOL));
assert(isequalRel(gsto,          3.007604046768e+00, TOL));
assert(isequalRel(no,            5.069384722250e-03, TOL));
# 9998

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 7.318036000000e-01,  1.118729629788e+04,  8.166674822762e-01,  9.971844782556e-03, 11801);
assert(isequalRel(ainv,          2.619649072783e-01, TOL));
assert(isequalRel(ao,            3.817305189423e+00, TOL));
assert(isequalRel(con41,         4.062531798299e-01, TOL));
assert(isequalRel(con42,        -1.343755299716e+00, TOL));
assert(isequalRel(cosio,         6.846539709542e-01, TOL));
assert(isequalRel(cosio2,        4.687510599433e-01, TOL));
assert(isequalRel(einv,          1.366486855216e+00, TOL));
assert(isequalRel(eccsq,         5.355365089730e-01, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        4.644634910270e-01, TOL));
assert(isequalRel(posq,          3.143525080235e+00, TOL));
assert(isequalRel(rp,            1.023787509504e+00, TOL));
assert(isequalRel(rteosq,        6.815155838475e-01, TOL));
assert(isequalRel(sinio,         7.288682597402e-01, TOL));
assert(isequalRel(gsto,          1.265125075736e+00, TOL));
assert(isequalRel(no,            9.971131604593e-03, TOL));
# 11801

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 1.156200000000e-03,  2.063002844893e+04,  1.996377411601e-01,  4.314022552805e-03, 14128);
assert(isequalRel(ainv,          1.498506014108e-01, TOL));
assert(isequalRel(ao,            6.673313223874e+00, TOL));
assert(isequalRel(con41,         1.882014339061e+00, TOL));
assert(isequalRel(con42,        -3.803357231769e+00, TOL));
assert(isequalRel(cosio,         9.801384832532e-01, TOL));
assert(isequalRel(cosio2,        9.606714463538e-01, TOL));
assert(isequalRel(einv,          8.649022660439e+02, TOL));
assert(isequalRel(eccsq,         1.336798440000e-06, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        9.999986632016e-01, TOL));
assert(isequalRel(posq,          4.453299032044e+01, TOL));
assert(isequalRel(rp,            6.665597539125e+00, TOL));
assert(isequalRel(rteosq,        9.999993316006e-01, TOL));
assert(isequalRel(sinio,         1.983142799855e-01, TOL));
assert(isequalRel(gsto,          4.943902851905e+00, TOL));
assert(isequalRel(no,            4.313874514604e-03, TOL));
# 14128

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 5.596327000000e-01,  2.060567415771e+04,  1.083685404539e+00,  2.131535163460e-02, 16925);
assert(isequalRel(ainv,          4.347456405863e-01, TOL));
assert(isequalRel(ao,            2.300195577928e+00, TOL));
assert(isequalRel(con41,        -3.427179465249e-01, TOL));
assert(isequalRel(con42,        -9.547008912508e-02, TOL));
assert(isequalRel(cosio,         4.680747993911e-01, TOL));
assert(isequalRel(cosio2,        2.190940178250e-01, TOL));
assert(isequalRel(einv,          1.786886291670e+00, TOL));
assert(isequalRel(eccsq,         3.131887589093e-01, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        6.868112410907e-01, TOL));
assert(isequalRel(posq,          2.495768607552e+00, TOL));
assert(isequalRel(rp,            1.012930916124e+00, TOL));
assert(isequalRel(rteosq,        8.287407562626e-01, TOL));
assert(isequalRel(sinio,         8.836888491856e-01, TOL));
assert(isequalRel(gsto,          2.298863660699e+00, TOL));
assert(isequalRel(no,            2.131732105901e-02, TOL));
# 16925

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 7.864447000000e-01,  2.045279166667e+04,  2.155725972308e-01,  1.077308058718e-03, 20413);
assert(isequalRel(ainv,          5.942468133044e-02, TOL));
assert(isequalRel(ao,            1.682802461219e+01, TOL));
assert(isequalRel(con41,         1.862731633406e+00, TOL));
assert(isequalRel(con42,        -3.771219389011e+00, TOL));
assert(isequalRel(cosio,         9.768540719074e-01, TOL));
assert(isequalRel(cosio2,        9.542438778022e-01, TOL));
assert(isequalRel(einv,          1.271545221171e+00, TOL));
assert(isequalRel(eccsq,         6.184952661581e-01, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        3.815047338419e-01, TOL));
assert(isequalRel(posq,          4.121602829256e+01, TOL));
assert(isequalRel(rp,            3.593713844463e+00, TOL));
assert(isequalRel(rteosq,        6.176606947523e-01, TOL));
assert(isequalRel(sinio,         2.139068072733e-01, TOL));
assert(isequalRel(gsto,          4.071885638059e-01, TOL));
assert(isequalRel(no,            1.077283640185e-03, TOL));
# 20413

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 7.421690000000e-01,  2.063002341244e+04,  1.085156717098e+00,  8.782060201923e-03, 21897);
assert(isequalRel(ainv,          2.407096096629e-01, TOL));
assert(isequalRel(ao,            4.154383372564e+00, TOL));
assert(isequalRel(con41,        -3.463657886765e-01, TOL));
assert(isequalRel(con42,        -8.939035220576e-02, TOL));
assert(isequalRel(cosio,         4.667741107229e-01, TOL));
assert(isequalRel(cosio2,        2.178780704412e-01, TOL));
assert(isequalRel(einv,          1.347402006821e+00, TOL));
assert(isequalRel(eccsq,         5.508148245610e-01, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        4.491851754390e-01, TOL));
assert(isequalRel(posq,          3.482282274182e+00, TOL));
assert(isequalRel(rp,            1.071128819331e+00, TOL));
assert(isequalRel(rteosq,        6.702127837031e-01, TOL));
assert(isequalRel(sinio,         8.843765767810e-01, TOL));
assert(isequalRel(gsto,          4.912171010008e+00, TOL));
assert(isequalRel(no,            8.782535528761e-03, TOL));
# 21897

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 3.087230000000e-02,  2.054846235912e+04,  1.084697695505e+00,  6.962749021661e-02, 22312);
assert(isequalRel(ainv,          9.572126270093e-01, TOL));
assert(isequalRel(ao,            1.044699967158e+00, TOL));
assert(isequalRel(con41,        -3.452285162292e-01, TOL));
assert(isequalRel(con42,        -9.128580628459e-02, TOL));
assert(isequalRel(cosio,         4.671800094791e-01, TOL));
assert(isequalRel(cosio2,        2.182571612569e-01, TOL));
assert(isequalRel(einv,          3.239149658432e+01, TOL));
assert(isequalRel(eccsq,         9.530989072900e-04, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        9.990469010927e-01, TOL));
assert(isequalRel(posq,          1.089318592280e+00, TOL));
assert(isequalRel(rp,            1.012447676362e+00, TOL));
assert(isequalRel(rteosq,        9.995233369425e-01, TOL));
assert(isequalRel(sinio,         8.841622242231e-01, TOL));
assert(isequalRel(gsto,          6.267076519855e+00, TOL));
assert(isequalRel(no,            6.964539419426e-02, TOL));
# 22312
# *** error: t:= 494.202867 *** code =   1

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 7.541712000000e-01,  2.063055909107e+04,  1.108345161540e+00,  8.581775554476e-03, 22674);
assert(isequalRel(ainv,          2.370374917914e-01, TOL));
assert(isequalRel(ao,            4.218741906365e+00, TOL));
assert(isequalRel(con41,        -4.028688855919e-01, TOL));
assert(isequalRel(con42,         4.781475986499e-03, TOL));
assert(isequalRel(cosio,         4.461431438481e-01, TOL));
assert(isequalRel(cosio2,        1.990437048027e-01, TOL));
assert(isequalRel(einv,          1.325958880424e+00, TOL));
assert(isequalRel(eccsq,         5.687741989094e-01, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        4.312258010906e-01, TOL));
assert(isequalRel(posq,          3.309599096075e+00, TOL));
assert(isequalRel(rp,            1.037088260351e+00, TOL));
assert(isequalRel(rteosq,        6.566778518349e-01, TOL));
assert(isequalRel(sinio,         8.949616166056e-01, TOL));
assert(isequalRel(gsto,          2.003968969617e+00, TOL));
assert(isequalRel(no,            8.582332516240e-03, TOL));
# 22674

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 7.258491000000e-01,  2.062945752052e+04,  1.230387309486e-01,  9.857037897024e-03, 23177);
assert(isequalRel(ainv,          2.599056204118e-01, TOL));
assert(isequalRel(ao,            3.847550500891e+00, TOL));
assert(isequalRel(con41,         1.954813125046e+00, TOL));
assert(isequalRel(con42,        -3.924688541744e+00, TOL));
assert(isequalRel(cosio,         9.924402794872e-01, TOL));
assert(isequalRel(cosio2,        9.849377083487e-01, TOL));
assert(isequalRel(einv,          1.377696824312e+00, TOL));
assert(isequalRel(eccsq,         5.268569159708e-01, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        4.731430840292e-01, TOL));
assert(isequalRel(posq,          3.314008747501e+00, TOL));
assert(isequalRel(rp,            1.054809432615e+00, TOL));
assert(isequalRel(rteosq,        6.878539699887e-01, TOL));
assert(isequalRel(sinio,         1.227285282698e-01, TOL));
assert(isequalRel(gsto,          1.346832291139e+00, TOL));
assert(isequalRel(no,            9.853789469239e-03, TOL));
# 23177

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 9.728298000000e-01,  1.637649999999e+04,  5.017647066558e-01,  3.189367114872e-04, 23333);
assert(isequalRel(ainv,          2.639566378141e-02, TOL));
assert(isequalRel(ao,            3.788501051845e+01, TOL));
assert(isequalRel(con41,         1.305993571902e+00, TOL));
assert(isequalRel(con42,        -2.843322619837e+00, TOL));
assert(isequalRel(cosio,         8.767351504117e-01, TOL));
assert(isequalRel(cosio2,        7.686645239674e-01, TOL));
assert(isequalRel(einv,          1.027929037536e+00, TOL));
assert(isequalRel(eccsq,         9.463978197680e-01, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        5.360218023196e-02, TOL));
assert(isequalRel(posq,          4.123820314506e+00, TOL));
assert(isequalRel(rp,            1.029343312788e+00, TOL));
assert(isequalRel(rteosq,        2.315214465918e-01, TOL));
assert(isequalRel(sinio,         4.809734670776e-01, TOL));
assert(isequalRel(gsto,          3.849219568564e+00, TOL));
assert(isequalRel(no,            3.189177226396e-04, TOL));
# 23333

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 5.782022000000e-01,  2.062576535463e+04,  1.209984410530e-01,  1.953881109593e-02, 23599);
assert(isequalRel(ainv,          4.100810857141e-01, TOL));
assert(isequalRel(ao,            2.438542119685e+00, TOL));
assert(isequalRel(con41,         1.956292061635e+00, TOL));
assert(isequalRel(con42,        -3.927153436058e+00, TOL));
assert(isequalRel(cosio,         9.926886154337e-01, TOL));
assert(isequalRel(cosio2,        9.854306872116e-01, TOL));
assert(isequalRel(einv,          1.729498780876e+00, TOL));
assert(isequalRel(eccsq,         3.343177840848e-01, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        6.656822159152e-01, TOL));
assert(isequalRel(posq,          2.635083805981e+00, TOL));
assert(isequalRel(rp,            1.028571701290e+00, TOL));
assert(isequalRel(rteosq,        8.158935077050e-01, TOL));
assert(isequalRel(sinio,         1.207034083544e-01, TOL));
assert(isequalRel(gsto,          3.217495488339e+00, TOL));
assert(isequalRel(no,            1.952919673236e-02, TOL));
# 23599

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 2.664000000000e-03,  2.063104061740e+04,  6.725800805485e-02,  4.397272140132e-03, 24208);
assert(isequalRel(ainv,          1.517719882533e-01, TOL));
assert(isequalRel(ao,            6.588831124297e+00, TOL));
assert(isequalRel(con41,         1.986449532035e+00, TOL));
assert(isequalRel(con42,        -3.977415886724e+00, TOL));
assert(isequalRel(cosio,         9.977390326858e-01, TOL));
assert(isequalRel(cosio2,        9.954831773449e-01, TOL));
assert(isequalRel(einv,          3.753753753754e+02, TOL));
assert(isequalRel(eccsq,         7.096896000000e-06, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        9.999929031040e-01, TOL));
assert(isequalRel(posq,          4.341207939592e+01, TOL));
assert(isequalRel(rp,            6.571278478182e+00, TOL));
assert(isequalRel(rteosq,        9.999964515457e-01, TOL));
assert(isequalRel(sinio,         6.720731102433e-02, TOL));
assert(isequalRel(gsto,          5.037771726289e+00, TOL));
assert(isequalRel(no,            4.397108759492e-03, TOL));
# 24208

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 1.765000000000e-04,  1.976268057285e+04,  6.981317007977e-06,  4.375160345672e-03, 25954);
assert(isequalRel(ainv,          1.512627677760e-01, TOL));
assert(isequalRel(ao,            6.611012179023e+00, TOL));
assert(isequalRel(con41,         1.999999999854e+00, TOL));
assert(isequalRel(con42,        -3.999999999756e+00, TOL));
assert(isequalRel(cosio,         9.999999999756e-01, TOL));
assert(isequalRel(cosio2,        9.999999999513e-01, TOL));
assert(isequalRel(einv,          5.665722379603e+03, TOL));
assert(isequalRel(eccsq,         3.115225000000e-08, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        9.999999688477e-01, TOL));
assert(isequalRel(posq,          4.370547930815e+01, TOL));
assert(isequalRel(rp,            6.609845335374e+00, TOL));
assert(isequalRel(rteosq,        9.999999844239e-01, TOL));
assert(isequalRel(sinio,         6.981317007921e-06, TOL));
assert(isequalRel(gsto,          4.036986406591e-01, TOL));
assert(isequalRel(no,            4.374997775854e-03, TOL));
# 25954

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 3.319000000000e-04,  2.056074503247e+04,  2.862339973271e-04,  4.375271959478e-03, 26900);
assert(isequalRel(ainv,          1.512653401930e-01, TOL));
assert(isequalRel(ao,            6.610899752212e+00, TOL));
assert(isequalRel(con41,         1.999999754210e+00, TOL));
assert(isequalRel(con42,        -3.999999590351e+00, TOL));
assert(isequalRel(cosio,         9.999999590350e-01, TOL));
assert(isequalRel(cosio2,        9.999999180701e-01, TOL));
assert(isequalRel(einv,          3.012955709551e+03, TOL));
assert(isequalRel(eccsq,         1.101576100000e-07, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        9.999998898424e-01, TOL));
assert(isequalRel(posq,          4.370398590514e+01, TOL));
assert(isequalRel(rp,            6.608705594584e+00, TOL));
assert(isequalRel(rteosq,        9.999999449212e-01, TOL));
assert(isequalRel(sinio,         2.862339934185e-04, TOL));
assert(isequalRel(gsto,          1.971276525301e+00, TOL));
assert(isequalRel(no,            4.375109379983e-03, TOL));
# 26900

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 5.602877000000e-01,  2.062885818871e+04,  1.195051373450e+00,  8.973503578612e-03, 26975);
assert(isequalRel(ainv,          2.441946885467e-01, TOL));
assert(isequalRel(ao,            4.095093165013e+00, TOL));
assert(isequalRel(con41,        -5.960087184875e-01, TOL));
assert(isequalRel(con42,         3.266811974792e-01, TOL));
assert(isequalRel(cosio,         3.669656121548e-01, TOL));
assert(isequalRel(cosio2,        1.346637605042e-01, TOL));
assert(isequalRel(einv,          1.784797346078e+00, TOL));
assert(isequalRel(eccsq,         3.139223067713e-01, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        6.860776932287e-01, TOL));
assert(isequalRel(posq,          7.893582846451e+00, TOL));
assert(isequalRel(rp,            1.800662834302e+00, TOL));
assert(isequalRel(rteosq,        8.282980702795e-01, TOL));
assert(isequalRel(sinio,         9.302345077967e-01, TOL));
assert(isequalRel(gsto,          3.853994597295e+00, TOL));
assert(isequalRel(no,            8.973959238486e-03, TOL));
# 26975

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 8.840000000000e-05,  2.063178615833e+04,  1.717897912141e+00,  6.263454709052e-02, 28057);
assert(isequalRel(ainv,          8.922044041800e-01, TOL));
assert(isequalRel(ao,            1.120819394429e+00, TOL));
assert(isequalRel(con41,        -9.355502628635e-01, TOL));
assert(isequalRel(con42,         8.925837714392e-01, TOL));
assert(isequalRel(cosio,        -1.465716402042e-01, TOL));
assert(isequalRel(cosio2,        2.148324571216e-02, TOL));
assert(isequalRel(einv,          1.131221719457e+04, TOL));
assert(isequalRel(eccsq,         7.814560000000e-09, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        9.999999921854e-01, TOL));
assert(isequalRel(posq,          1.256236095295e+00, TOL));
assert(isequalRel(rp,            1.120720313995e+00, TOL));
assert(isequalRel(rteosq,        9.999999960927e-01, TOL));
assert(isequalRel(sinio,         9.892000577678e-01, TOL));
assert(isequalRel(gsto,          3.451783621544e+00, TOL));
assert(isequalRel(no,            6.267239862091e-02, TOL));
# 28057

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 4.850600000000e-03,  2.062957071136e+04,  9.552152089580e-01,  8.751201646284e-03, 28129);
assert(isequalRel(ainv,          2.401367388842e-01, TOL));
assert(isequalRel(ao,            4.164294079475e+00, TOL));
assert(isequalRel(con41,         2.868387191267e-04, TOL));
assert(isequalRel(con42,        -6.671447311985e-01, TOL));
assert(isequalRel(cosio,         5.774330664585e-01, TOL));
assert(isequalRel(cosio2,        3.334289462397e-01, TOL));
assert(isequalRel(einv,          2.061600626727e+02, TOL));
assert(isequalRel(eccsq,         2.352832036000e-05, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        9.999764716796e-01, TOL));
assert(isequalRel(posq,          1.734052916450e+01, TOL));
assert(isequalRel(rp,            4.144094754613e+00, TOL));
assert(isequalRel(rteosq,        9.999882357706e-01, TOL));
assert(isequalRel(sinio,         8.164380281199e-01, TOL));
assert(isequalRel(gsto,          2.059978512381e+00, TOL));
assert(isequalRel(no,            8.751201528747e-03, TOL));
# 28129

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 2.487000000000e-03,  2.062121788666e+04,  1.134423871224e+00,  7.190131350005e-02, 28350);
assert(isequalRel(ainv,          9.780078361195e-01, TOL));
assert(isequalRel(ao,            1.022486694961e+00, TOL));
assert(isequalRel(con41,        -4.640891584380e-01, TOL));
assert(isequalRel(con42,         1.068152640634e-01, TOL));
assert(isequalRel(cosio,         4.226546429265e-01, TOL));
assert(isequalRel(cosio2,        1.786369471873e-01, TOL));
assert(isequalRel(einv,          4.020908725372e+02, TOL));
assert(isequalRel(eccsq,         6.185169000000e-06, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        9.999938148310e-01, TOL));
assert(isequalRel(posq,          1.045466108483e+00, TOL));
assert(isequalRel(rp,            1.019943770551e+00, TOL));
assert(isequalRel(rteosq,        9.999969074107e-01, TOL));
assert(isequalRel(sinio,         9.062908213221e-01, TOL));
assert(isequalRel(gsto,          5.982608943601e+00, TOL));
assert(isequalRel(no,            7.192721985507e-02, TOL));
# 28350
# *** error: t:= 1560.000000 *** code =   1

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 6.249053000000e-01,  2.063181079184e+04,  4.977679026688e-01,  1.655781478256e-02, 28623);
assert(isequalRel(ainv,          3.672775780112e-01, TOL));
assert(isequalRel(ao,            2.722736316807e+00, TOL));
assert(isequalRel(con41,         1.316080099741e+00, TOL));
assert(isequalRel(con42,        -2.860133499568e+00, TOL));
assert(isequalRel(cosio,         8.786504992963e-01, TOL));
assert(isequalRel(cosio2,        7.720266999137e-01, TOL));
assert(isequalRel(einv,          1.600242468739e+00, TOL));
assert(isequalRel(eccsq,         3.905066339681e-01, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        6.094933660319e-01, TOL));
assert(isequalRel(posq,          2.753906139243e+00, TOL));
assert(isequalRel(rp,            1.021283961932e+00, TOL));
assert(isequalRel(rteosq,        7.807005610552e-01, TOL));
assert(isequalRel(sinio,         4.774654962260e-01, TOL));
assert(isequalRel(gsto,          3.606984295022e+00, TOL));
assert(isequalRel(no,            1.655279727416e-02, TOL));
# 28623

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 3.350000000000e-05,  2.063046683397e+04,  3.316125578789e-05,  4.375111781885e-03, 28626);
assert(isequalRel(ainv,          1.512616484965e-01, TOL));
assert(isequalRel(ao,            6.611061098038e+00, TOL));
assert(isequalRel(con41,         1.999999996701e+00, TOL));
assert(isequalRel(con42,        -3.999999994502e+00, TOL));
assert(isequalRel(cosio,         9.999999994502e-01, TOL));
assert(isequalRel(cosio2,        9.999999989003e-01, TOL));
assert(isequalRel(einv,          2.985074626866e+04, TOL));
assert(isequalRel(eccsq,         1.122250000000e-09, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        9.999999988778e-01, TOL));
assert(isequalRel(posq,          4.370612874389e+01, TOL));
assert(isequalRel(rp,            6.610839627491e+00, TOL));
assert(isequalRel(rteosq,        9.999999994389e-01, TOL));
assert(isequalRel(sinio,         3.316125578181e-05, TOL));
assert(isequalRel(gsto,          1.422713433736e+00, TOL));
assert(isequalRel(no,            4.374949216286e-03, TOL));
# 28626

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 3.039550000000e-02,  2.042202012661e+04,  1.683781961252e+00,  7.182099414601e-02, 28872);
assert(isequalRel(ainv,          9.775310460243e-01, TOL));
assert(isequalRel(ao,            1.022985412143e+00, TOL));
assert(isequalRel(con41,        -9.618654265547e-01, TOL));
assert(isequalRel(con42,         9.364423775911e-01, TOL));
assert(isequalRel(cosio,        -1.127453967210e-01, TOL));
assert(isequalRel(cosio2,        1.271152448177e-02, TOL));
assert(isequalRel(einv,          3.289960684970e+01, TOL));
assert(isequalRel(eccsq,         9.238864202500e-04, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        9.990761135798e-01, TOL));
assert(isequalRel(posq,          1.044566354001e+00, TOL));
assert(isequalRel(rp,            9.918912590484e-01, TOL));
assert(isequalRel(rteosq,        9.995379500448e-01, TOL));
assert(isequalRel(sinio,         9.936239105005e-01, TOL));
assert(isequalRel(gsto,          1.313288306073e+00, TOL));
assert(isequalRel(no,            7.187462824278e-02, TOL));
# 28872

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 1.584800000000e-03,  2.062426783845e+04,  1.438653958468e+00,  6.952270688787e-02, 29141);
assert(isequalRel(ainv,          9.565366465641e-01, TOL));
assert(isequalRel(ao,            1.045438252253e+00, TOL));
assert(isequalRel(con41,        -9.479193821401e-01, TOL));
assert(isequalRel(con42,         9.131989702335e-01, TOL));
assert(isequalRel(cosio,         1.317581342965e-01, TOL));
assert(isequalRel(cosio2,        1.736020595330e-02, TOL));
assert(isequalRel(einv,          6.309944472489e+02, TOL));
assert(isequalRel(eccsq,         2.511591040000e-06, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        9.999974884090e-01, TOL));
assert(isequalRel(posq,          1.092935649238e+00, TOL));
assert(isequalRel(rp,            1.043781441711e+00, TOL));
assert(isequalRel(rteosq,        9.999987442037e-01, TOL));
assert(isequalRel(sinio,         9.912818943402e-01, TOL));
assert(isequalRel(gsto,          6.574767745587e-02, TOL));
assert(isequalRel(no,            6.957163219428e-02, TOL));
# 29141
# *** error: t:= 440.000000 *** code =   6

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 2.025790000000e-02,  2.063128732010e+04,  8.998830356820e-01,  6.867101959232e-02, 29238);
assert(isequalRel(ainv,          9.481900688830e-01, TOL));
assert(isequalRel(ao,            1.054640870873e+00, TOL));
assert(isequalRel(con41,         1.595385835540e-01, TOL));
assert(isequalRel(con42,        -9.325643059233e-01, TOL));
assert(isequalRel(cosio,         6.217015853162e-01, TOL));
assert(isequalRel(cosio2,        3.865128611847e-01, TOL));
assert(isequalRel(einv,          4.936345820643e+01, TOL));
assert(isequalRel(eccsq,         4.103825124100e-04, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        9.995896174876e-01, TOL));
assert(isequalRel(posq,          1.111354643684e+00, TOL));
assert(isequalRel(rp,            1.033276061575e+00, TOL));
assert(isequalRel(rteosq,        9.997947876877e-01, TOL));
assert(isequalRel(sinio,         7.832541980834e-01, TOL));
assert(isequalRel(gsto,          3.089091736929e-01, TOL));
assert(isequalRel(no,            6.866301602203e-02, TOL));
# 29238

[ainv, ao, con41, con42, cosio, cosio2, einv,...
 eccsq,  method, omeosq, posq, rp, rteosq, sinio,...
 gsto, no]...
= initl( 8.673100000000e-03,  1.123298708465e+04,  1.271358913676e+00,  7.006731262088e-02, 88888);
assert(isequalRel(ainv,          9.614298171958e-01, TOL));
assert(isequalRel(ao,            1.040117523000e+00, TOL));
assert(isequalRel(con41,        -7.389556198424e-01, TOL));
assert(isequalRel(con42,         5.649260330707e-01, TOL));
assert(isequalRel(cosio,         2.949827001467e-01, TOL));
assert(isequalRel(cosio2,        8.701479338586e-02, TOL));
assert(isequalRel(einv,          1.152990280292e+02, TOL));
assert(isequalRel(eccsq,         7.522266361000e-05, TOL));
assert(isequal(method,           "n"));
assert(isequalRel(omeosq,        9.999247773364e-01, TOL));
assert(isequalRel(posq,          1.081681709329e+00, TOL));
assert(isequalRel(rp,            1.031096479711e+00, TOL));
assert(isequalRel(rteosq,        9.999623879609e-01, TOL));
assert(isequalRel(sinio,         9.555025937244e-01, TOL));
assert(isequalRel(gsto,          1.082901416696e-01, TOL));
assert(isequalRel(no,            7.010615566653e-02, TOL));
# 88888

