function fibsq32()
#
#  test polynomial suggested by Goedecker
#  square of Fibocacci polynomial
#
    n = 32;
    p = fib(n)
    p = conv(p,p)
    z = [-.9665892387165332;
         -.9485092065867169-.1867222182129376*im;
         -.9485092065867169+.1867222182129376*im;
         -.8949223987510440-.3665688458356083*im;
         -.8949223987510440+.3665688458356083*im;
         -.8077651569138829-.5329094943253877*im;
         -.8077651569138829+.5329094943253877*im;
         -.6901871862102008-.6795950061744879*im;
         -.6901871862102008+.6795950061744879*im;
         -.5464383280299697-.8011754158225671*im;
         -.5464383280299697+.8011754158225671*im;
         -.3817161237254617-.8930914814237271*im;
         -.3817161237254617+.8930914814237271*im;
         -.2019801534661303-.9518321713734013*im;
         -.2019801534661303+.9518321713734013*im;
         -.1374081948052525e-1-.9750513807247308*im;
         -.1374081948052525e-1+.9750513807247308*im;
         .1761678077924798-.9616381291177512*im;
         .1761678077924798+.9616381291177512*im;
         .3608186423377573-.9117356125589683*im;
         .3608186423377573+.9117356125589683*im;
         .5334081767292508-.8267064135873833*im;
         .5334081767292508+.8267064135873833*im;
         .6874310990889808-.7090471025782057*im;
         .6874310990889808+.7090471025782057*im;
         .8167660514846921-.5622778345619681*im;
         .8167660514846921+.5622778345619681*im;
         .9156129648696007-.3909102172601782*im;
         .9156129648696007+.3909102172601782*im;
         .9783492503358520-.2007808902007197*im;
         .9783492503358520+.2007808902007197*im;
         1.999999999767169];
    z = [z 2*ones(n)]
    p, PolyZeros(z)
end
