20240314  202750.080
 
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.892701E-02     1.14567         0.00000       -0.145717E-15
    0.892701E-02    -1.14567         0.00000        0.145717E-15
    -1.00000         0.00000       -0.888178E-15     0.00000    
    -4.67405         0.00000       -0.408562E-13     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.809017        0.587785        0.444089E-15   -0.499600E-15
   -0.809017       -0.587785        0.444089E-15    0.499600E-15
    0.309017        0.951057       -0.888178E-15   -0.199840E-14
    0.309017       -0.951057       -0.888178E-15    0.199840E-14
     1.00000         0.00000       -0.277556E-14     0.00000    
 
toms493_test():
  Normal end of execution.
 
20240314  202750.080