20240318  093346.580
 
toms493_test():
  Fortran77 version
  Test toms493().
 
test01():
  p(x) = x^4 + 5.6562x^3 + 5.8854x^2
             + 7.3646x   + 6.1354
 
  Approximate roots:
    -1.001,
    -4.7741, 
     0.0089 + 1.1457 i,
     0.0089 - 1.1457 i.

      real(z)         imag(z)       real(f(z))      imag(f(z))

    0.892719E-02     1.14567       -0.317724E-05    0.485611E-06
    0.892719E-02    -1.14567       -0.317724E-05   -0.485611E-06
    -1.00000         0.00000       -0.360489E-05     0.00000    
    -4.67405         0.00000       -0.516219E-05     0.00000    
 
test02():
  p(x) = x^4 - 8 x^3 + 24 x^2 - 32 x + 16
 
  Exact roots:
    2 (multiplicity 4)

      real(z)         imag(z)       real(f(z))      imag(f(z))

     2.00000         0.00000         0.00000         0.00000    
     2.00000         0.00000         0.00000         0.00000    
     2.00000         0.00000         0.00000         0.00000    
     2.00000         0.00000         0.00000         0.00000    
 
test03():
  p(x) = x^5 - 1
 
  Exact roots:
    0.3090 + 0.9510 i
   -0.8090 + 0.5877 i
   -0.8090 - 0.5877 i
    0.3090 - 0.9510 i
    1

      real(z)         imag(z)       real(f(z))      imag(f(z))

   -0.809018        0.587786        0.463270E-05   -0.145367E-05
   -0.809018       -0.587786        0.463270E-05    0.145367E-05
    0.309018        0.951057        0.143467E-05   -0.249098E-05
    0.309018       -0.951057        0.143467E-05    0.249098E-05
     1.00000         0.00000        0.162164E-06     0.00000    
 
toms493_test():
  Normal end of execution.
 
20240318  093346.580