############################################################ # Building 'newtsqrt' and sample output on SGIs ############################################################ einstein% pwd; ls /usr2/people/phy329/nonlin/ex2 Makefile newtsqrt.f einstein% make f77 -g -c newtsqrt.f f77 -g -L/usr/local/lib newtsqrt.o -lp329f -o newtsqrt einstein% newtsqrt usage: newtsqrt [] einstein% newtsqrt 10.0 Iter xn dxn log(dxn) rn log(rn) 1 3.5000000000000000E+00 1.500E+00 0.41 1.500E+01 2.71 2 3.1785714285714284E+00 3.214E-01 -1.13 2.250E+00 0.81 3 3.1623194221508828E+00 1.625E-02 -4.12 1.033E-01 -2.27 4 3.1622776604441363E+00 4.176E-05 -10.08 2.641E-04 -8.24 5 3.1622776601683795E+00 2.758E-10 -22.01 1.744E-09 -20.17 10.00000000000000 3.162277660168380 einstein% newtsqrt 10.0 1.0e-15 Iter xn dxn log(dxn) rn log(rn) 1 3.5000000000000000E+00 1.500E+00 0.41 1.500E+01 2.71 2 3.1785714285714284E+00 3.214E-01 -1.13 2.250E+00 0.81 3 3.1623194221508828E+00 1.625E-02 -4.12 1.033E-01 -2.27 4 3.1622776604441363E+00 4.176E-05 -10.08 2.641E-04 -8.24 5 3.1622776601683795E+00 2.758E-10 -22.01 1.744E-09 -20.17 6 3.1622776601683791E+00 2.809E-16 -35.81 1.776E-15 -33.96 10.00000000000000 3.162277660168379 einstein% newtsqrt 0.5 Iter xn dxn log(dxn) rn log(rn) 1 1.1250000000000000E+00 -8.750E-01 -0.13 -4.375E-01 -0.83 2 7.8472222222222221E-01 3.403E-01 -1.08 7.656E-01 -0.27 3 7.1094518190757128E-01 7.378E-02 -2.61 1.158E-01 -2.16 4 7.0711714297003669E-01 3.828E-03 -5.57 5.443E-03 -5.21 5 7.0710678126246607E-01 1.036E-05 -11.48 1.465E-05 -11.13 6 7.0710678118654757E-01 7.592E-11 -23.30 1.074E-10 -22.95 0.5000000000000000 0.707106781186547