function fib30()
#
#  test polynomial suggested by Goedecker
#
    n = 30;
    p = fib(30)
    z = [-.9644254204388659;
         -.9439446810207413-.1984294062681023*im; -.9439446810207413+.1984294062681023*im;
         -.8833406981383851-.3885717512211873*im; -.8833406981383851+.3885717512211873*im;
         -.7850939960122262-.5624744046156477*im; -.7850939960122262+.5624744046156477*im;
         -.6532262984053521-.7128393846266873*im; -.6532262984053521+.7128393846266873*im;
         -.4931368697016981-.8333156264808382*im; -.4931368697016981+.8333156264808382*im;
         -.3113834310599106-.9187504784857605*im; -.3113834310599106+.9187504784857605*im;
         -.1154175582352257-.9653888511274815*im; -.1154175582352257+.9653888511274815*im;
         .8671260421202705e-1-.9710099019085114*im; .8671260421202705e-1+.9710099019085114*im;
         .2866757308842761-.9349927196966469*im; .2866757308842761+.9349927196966469*im;
         .4761639390711853-.8583045796885356*im; .4761639390711853+.8583045796885356*im;
         .6471568601471165-.7434108410885516*im; .6471568601471165+.7434108410885516*im;
         .7920805342279742-.5941279591518842*im; .7920805342279742+.5941279591518842*im;
         .9037637390940364-.4155315939426088*im; .9037637390940364+.4155315939426088*im;
         .9752028356220179-.2142999451221434*im; .9752028356220179+.2142999451221434*im;
         1.999999999068677];

    z = [z ones(n)]
    p, PolyZeros(z)
end
