07 June 2024   9:09:33.169 PM
 
toms493_test():
  Fortran90 version
  Test toms493().
 
test01():
  p(x) = x^4 + 5.6562x^3 + 5.8854x^2
             + 7.3646x   + 6.1354
 
  Approximate roots:
    -1.0,
    -4.674054017161709, 
     0.0089 + 1.1457 i,
     0.0089 - 1.1457 i.

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

    0.8927008581E-02     1.145674888         0.00000       -0.145717E-15
    0.8927008581E-02    -1.145674888         0.00000        0.145717E-15
    -1.000000000         0.000000000       -0.888178E-15     0.00000    
    -4.674054017         0.000000000       -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.
 
07 June 2024   9:09:33.169 PM