1      subroutine graeff(a,b)
2      real*8 a(11),b(11)
3c
4c     --  Graeffe Root-Squaring Method   --
5c     --  to Find Roots of a Polynomial  --
6c
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
23