==================================================== @(#) plfadr.f 1.6 from 96/07/19 (mexx 1.2) ==================================================== Start example Pendulum (default options) Enter number of bodies (default 8): Given: 8 NS, NF, T, H, K, KOPT: 0 0 0.0000000E+00 1.0000000E-03 0 3 NS, NF, T, H, K, KOPT: 1 5 1.0000000E-03 1.4444444E-02 2 3 NS, NF, T, H, K, KOPT: 2 9 1.5444444E-02 2.6585170E-02 2 3 NS, NF, T, H, K, KOPT: 3 20 4.2029614E-02 5.7276012E-02 4 4 NS, NF, T, H, K, KOPT: 4 36 9.9305627E-02 1.0185275E-01 5 5 Restart step with H = 0.0509263768 KOPT = 5 NS, NF, T, H, K, KOPT: 5 67 1.5023200E-01 5.0926377E-02 5 5 NS, NF, T, H, K, KOPT: 6 83 2.0115838E-01 6.9652183E-02 5 6 Restart step with H = 0.0320663714 KOPT = 6 NS, NF, T, H, K, KOPT: 7 128 2.3322475E-01 3.2066371E-02 6 5 NS, NF, T, H, K, KOPT: 8 144 2.6529112E-01 5.1687711E-02 5 6 NS, NF, T, H, K, KOPT: 9 176 3.1697883E-01 4.9092419E-02 7 7 Restart step with H = 0.0245462094 KOPT = 7 NS, NF, T, H, K, KOPT: 10 214 3.4152504E-01 2.4546209E-02 7 6 NS, NF, T, H, K, KOPT: 11 230 3.6607125E-01 5.1214897E-02 5 6 Restart step with H = 0.0100487682 KOPT = 5 NS, NF, T, H, K, KOPT: 12 261 3.7612002E-01 1.0048768E-02 5 4 NS, NF, T, H, K, KOPT: 13 272 3.8616879E-01 2.2368527E-02 4 5 NS, NF, T, H, K, KOPT: 14 288 4.0853732E-01 2.8529224E-02 5 6 Restart step with H = 0.0169907577 KOPT = 6 NS, NF, T, H, K, KOPT: 15 342 4.2552807E-01 1.6990758E-02 6 6 NS, NF, T, H, K, KOPT: 16 358 4.4251883E-01 2.0158450E-02 5 6 NS, NF, T, H, K, KOPT: 17 381 4.6267728E-01 2.9353226E-02 6 7 Restart step with H = 0.0146766131 KOPT = 7 NS, NF, T, H, K, KOPT: 18 419 4.7735390E-01 1.4676613E-02 7 6 NS, NF, T, H, K, KOPT: 19 435 4.9203051E-01 2.7823195E-02 5 6 NS, NF, T, H, K, KOPT: 20 458 5.1985370E-01 4.3880056E-02 6 7 Restart step with H = 0.0219400281 KOPT = 7 NS, NF, T, H, K, KOPT: 21 512 5.4179373E-01 2.1940028E-02 7 6 NS, NF, T, H, K, KOPT: 22 528 5.6373376E-01 3.6099143E-02 5 6 Restart step with H = 0.021122368 KOPT = 6 NS, NF, T, H, K, KOPT: 23 582 5.8485613E-01 2.1122368E-02 6 6 NS, NF, T, H, K, KOPT: 24 598 6.0597850E-01 4.0023118E-02 5 6 Restart step with H = 0.020011559 KOPT = 6 NS, NF, T, H, K, KOPT: 25 627 6.2599005E-01 2.0011559E-02 6 6 NS, NF, T, H, K, KOPT: 26 650 6.4600161E-01 1.8779578E-02 6 6 NS, NF, T, H, K, KOPT: 27 666 6.6478119E-01 4.4527116E-02 5 6 Restart step with H = 0.00970735342 KOPT = 5 NS, NF, T, H, K, KOPT: 28 697 6.7448854E-01 9.7073534E-03 5 5 NS, NF, T, H, K, KOPT: 29 708 6.8419590E-01 1.2077496E-02 4 5 NS, NF, T, H, K, KOPT: 30 719 6.9627339E-01 1.5183764E-02 4 5 NS, NF, T, H, K, KOPT: 31 730 7.1145716E-01 1.8756060E-02 4 5 NS, NF, T, H, K, KOPT: 32 746 7.3021322E-01 3.7426177E-02 5 6 Restart step with H = 0.0187130886 KOPT = 6 NS, NF, T, H, K, KOPT: 33 784 7.4892631E-01 1.8713089E-02 6 6 NS, NF, T, H, K, KOPT: 34 800 7.6763939E-01 4.4369467E-02 5 6 NS, NF, T, H, K, KOPT: 35 816 8.1200886E-01 9.7619454E-02 5 6 NS, NF, T, H, K, KOPT: 36 848 9.0962832E-01 8.2697314E-02 7 6 NS, NF, T, H, K, KOPT: 37 864 9.9232563E-01 1.1471527E-01 5 6 NS, NF, T, H, K, KOPT: 38 871 1.0000000E+00 3.3602523E-02 3 4 o.k. Position values at the end point: 1.0E+00 P(I),I=1,2,.. -4.8603595155E-02 -3.9293008768E-02 -8.9093120552E-01 -1.5757224018E-01 -9.4782027957E-02 -1.3086686161E+00 -2.7587184759E-01 -8.7529434072E-02 -1.9553847637E+00 -3.8611968484E-01 -2.9881841998E-02 -2.1497922695E+00 -4.9857855921E-01 2.1310010297E-02 -1.8461484708E+00 -6.1376399352E-01 8.6909507078E-03 -1.0772059154E+00 -7.2845646578E-01 -3.9570730107E-02 -1.2677845658E+00 -8.3978782081E-01 -9.3371015204E-02 -9.7350062917E-01 C=================================== TOL 1.0E-03 Statistics: NSTEP, NACCPT, NREJCT 49 38 11 NFPROB, NFCN, NGCN 2080 871 44 NMULT, NDEC, NSOL 1126 1160 1165 C=================================== End of example Pendulum (default options) Bye