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