share/contrib/gentran/test/graeffe.output

   1       subroutine graeff(a,b)
   2       real*8 a(11),b(11)
   3 c
   4 c     --  Graeffe Root-Squaring Method   --
   5 c     --  to Find Roots of a Polynomial  --
   6 c
   7       b(1)=a(11)**2
   8       b(2)=-a(10)**2+2.0*a(9)*a(11)
   9       b(3)=a(9)**2-2.0*a(8)*a(10)+2.0*a(7)*a(11)
  10       b(4)=-a(8)**2+2.0*a(7)*a(9)-2.0*a(6)*a(10)+2.0*a(5)*a(11)
  11       b(5)=a(7)**2-2.0*a(6)*a(8)+2.0*a(5)*a(9)-2.0*a(4)*a(10)+2.0*a(3)*a
  12      . (11)
  13       b(6)=-a(6)**2+2.0*a(5)*a(7)-2.0*a(4)*a(8)+2.0*a(3)*a(9)-2.0*a(2)*a
  14      . (10)+2.0*a(1)*a(11)
  15       b(7)=a(5)**2-2.0*a(4)*a(6)+2.0*a(3)*a(7)-2.0*a(2)*a(8)+2.0*a(1)*a(
  16      . 9)
  17       b(8)=-a(4)**2+2.0*a(3)*a(5)-2.0*a(2)*a(6)+2.0*a(1)*a(7)
  18       b(9)=a(3)**2-2.0*a(2)*a(4)+2.0*a(1)*a(5)
  19       b(10)=-a(2)**2+2.0*a(1)*a(3)
  20       b(11)=a(1)**2
  21       return
  22       end