23 October 2023 6:26:53.246 PM mus_test(): FORTRAN77 version Test mus(). musl_test() FORTRAN77 version Test a linear problem. Condition number = 1.33153 Amplification factor = 2.21335 I T APPROX. SOL. EXACT SOL. ABS. ERROR 1 0.0000 0.100000001D+01 0.100000000D+01 -0.120757435D-07 0.100000001D+01 0.100000000D+01 -0.149754595D-07 0.100000001D+01 0.100000000D+01 -0.130719142D-07 2 0.6000 0.182211882D+01 0.182211880D+01 -0.230910786D-07 0.182211882D+01 0.182211880D+01 -0.186150277D-07 0.182211880D+01 0.182211880D+01 0.276482237D-08 3 1.2000 0.332011694D+01 0.332011692D+01 -0.162950098D-07 0.332011695D+01 0.332011692D+01 -0.299702663D-07 0.332011690D+01 0.332011692D+01 0.253191095D-07 4 1.8000 0.604964745D+01 0.604964746D+01 0.189447835D-07 0.604964752D+01 0.604964746D+01 -0.521154053D-07 0.604964743D+01 0.604964746D+01 0.319208597D-07 5 2.4000 0.110231763D+02 0.110231764D+02 0.450974760D-07 0.110231764D+02 0.110231764D+02 -0.360646197D-07 0.110231764D+02 0.110231764D+02 0.539663603D-08 6 3.0000 0.200855369D+02 0.200855369D+02 0.716162063D-08 0.200855369D+02 0.200855369D+02 -0.169556280D-07 0.200855369D+02 0.200855369D+02 -0.136451916D-07 7 3.6000 0.365982345D+02 0.365982344D+02 -0.159335229D-07 0.365982345D+02 0.365982344D+02 -0.192572429D-07 0.365982344D+02 0.365982344D+02 -0.500946129D-08 8 4.2000 0.666863311D+02 0.666863310D+02 -0.193062135D-07 0.666863311D+02 0.666863310D+02 -0.313411306D-07 0.666863310D+02 0.666863310D+02 0.170772125D-07 9 4.8000 0.121510418D+03 0.121510418D+03 0.102887725D-07 0.121510418D+03 0.121510418D+03 -0.503275857D-07 0.121510417D+03 0.121510418D+03 0.372507571D-07 10 5.4000 0.221406416D+03 0.221406416D+03 0.489649210D-07 0.221406416D+03 0.221406416D+03 -0.360826959D-07 0.221406416D+03 0.221406416D+03 0.207050732D-07 11 6.0000 0.403428793D+03 0.403428793D+03 0.120757591D-07 0.403428793D+03 0.403428793D+03 0.149754555D-07 0.403428793D+03 0.403428793D+03 0.130718831D-07 musl_test(): Normal end of execution. musn_test(): FORTRAN77 version Test a nonlinear problem. ******************** START TOLERANCE : 1.00000D-02 SHOOTING POINTS : 0.00000D+00 0.10000D+00 0.20000D+00 0.30000D+00 0.40000D+00 0.50000D+00 0.60000D+00 0.70000D+00 0.80000D+00 0.90000D+00 0.10000D+01 NORMALIZED RESIDUE : 8.107378887D+00 COND.NUMBER = 1.77127D+01 AMPLI.FACTOR = 9.12773D+00 PREDICTED DAMPING FACTOR : 1.0000000 NORMALIZED RESIDUE : 3.894499706D+00 COND.NUMBER = 1.98034D+01 AMPLI.FACTOR = 1.07774D+01 ************************************************************ **************** 1 ITERATION COMPLETED **************** ******************** FULL NEWTON STEP ******************** ************************************************************ PREDICTED DAMPING FACTOR : 1.0000000 NORMALIZED RESIDUE : 7.332734293D-01 COND.NUMBER = 2.53341D+01 AMPLI.FACTOR = 1.17548D+01 ************************************************************ **************** 2 ITERATION COMPLETED **************** ******************** FULL NEWTON STEP ******************** ************************************************************ PREDICTED DAMPING FACTOR : 1.0000000 NORMALIZED RESIDUE : 6.999174185D-02 COND.NUMBER = 2.88354D+01 AMPLI.FACTOR = 1.19275D+01 ************************************************************ **************** 3 ITERATION COMPLETED **************** ******************** FULL NEWTON STEP ******************** ************************************************************ PREDICTED DAMPING FACTOR : 1.0000000 NORMALIZED RESIDUE : 6.005897918D-04 COND.NUMBER = 2.91954D+01 AMPLI.FACTOR = 1.19391D+01 ************************************************************ **************** 4 ITERATION COMPLETED **************** ******************** FULL NEWTON STEP ******************** ************************************************************ ******************** NEW TOLERANCE : 1.00000D-04 NORMALIZED RESIDUE : 6.005897918D-04 COND.NUMBER = 2.91954D+01 AMPLI.FACTOR = 1.19391D+01 PREDICTED DAMPING FACTOR : 1.0000000 NORMALIZED RESIDUE : 4.540102662D-08 COND.NUMBER = 2.91986D+01 AMPLI.FACTOR = 1.19392D+01 ************************************************************ **************** 5 ITERATION COMPLETED **************** ******************** FULL NEWTON STEP ******************** ************************************************************ ******************** NEW TOLERANCE : 1.00000D-06 NORMALIZED RESIDUE : 4.540102662D-08 COND.NUMBER = 2.91986D+01 AMPLI.FACTOR = 1.19392D+01 RETURN FROM MUSN : NO ITERATIONS = 5 MUSN: IERROR = 0 A = 0.0000 B = 1.0000 REQUIRED TOLERANCE = 1.00000D-06 START TOLERANCE = 1.00000D-02 CONDITION NUMBER = 2.91986D+01 AMPLIFICATION FACTOR = 1.19392D+01 K-PARTITIONING = 2 I T Y1 Y2 Y3 Y4 Y5 1 0.0000 1.00000D+00 1.00000D+00 1.00000D+00 -1.00000D+01 9.67963D-01 2 0.1000 1.00701D+00 9.93036D-01 1.27014D+00 -9.99304D+00 1.24622D+00 3 0.2000 1.02560D+00 9.75042D-01 1.47051D+00 -9.97504D+00 1.45280D+00 4 0.3000 1.05313D+00 9.49550D-01 1.61931D+00 -9.94955D+00 1.60610D+00 5 0.4000 1.08796D+00 9.19155D-01 1.73140D+00 -9.91915D+00 1.72137D+00 6 0.5000 1.12900D+00 8.85737D-01 1.81775D+00 -9.88574D+00 1.80994D+00 7 0.6000 1.17554D+00 8.50676D-01 1.88576D+00 -9.85068D+00 1.87957D+00 8 0.7000 1.22696D+00 8.15025D-01 1.93990D+00 -9.81503D+00 1.93498D+00 9 0.8000 1.28262D+00 7.79653D-01 1.98190D+00 -9.77965D+00 1.97819D+00 10 0.9000 1.34161D+00 7.45374D-01 2.01050D+00 -9.74537D+00 2.00827D+00 11 1.0000 1.40232D+00 7.13102D-01 2.02032D+00 -9.71310D+00 2.02032D+00 musn_test(): Normal end of execution. mus_test(): Normal end of execution. 23 October 2023 6:26:53.247 PM