function fibsq24()
#
#  test polynomial suggested by Goedecker
#  square of Fibocacci polynomial
#
    n = 24;
    p = fib(n)
    p = conv(p, p)
    z = [-.9558467190743565;
         -.9244083656460181-.2442575811931259*im;
         -.9244083656460181+.2442575811931259*im;
         -.8320677914606713-.4727784060949767*im;
         -.8320677914606713+.4727784060949767*im;
         -.6846250033032301-.6707962641470166*im;
         -.6846250033032301+.6707962641470166*im;
         -.4913433576989534-.8254218906657225*im;
         -.4913433576989534+.8254218906657225*im;
         -.2643730603203030-.9264240880183555*im;
         -.2643730603203030+.9264240880183555*im;
         -.1800084671534713e-1-.9668298465915332*im;
         -.1800084671534713e-1+.9668298465915332*im;
         .2322158170552186-.9432943625623267*im;
         .2322158170552186+.9432943625623267*im;
         .4703449187206489-.8561989307649443*im;
         .4703449187206489+.8561989307649443*im;
         .6808619003199168-.7094507040001145*im;
         .6808619003199168+.7094507040001145*im;
         .8490194865023282-.5100588510854889*im;
         .8490194865023282+.5100588510854889*im;
         .9602996918859324-.2681678908797330*im;
         .9602996918859324+.2681678908797330*im;
         1.999999940395313];
    z = [z 2*ones(n)]
    p, PolyZeros(z)
end
