function lgd50()
#
#  test polynomial suggested by Goedecker
#   Legendre polynomial of degree 50
#
    p = lgd(50);
    z = [-.9988664044200711; -.9940319694320907; -.9853540840480059; -.9728643851066921;
         -.9566109552428079; -.9366566189448779; -.9130785566557919; -.8859679795236130;
         -.8554297694299461; -.8215820708593359; -.7845558329003993; -.7444943022260685;
         -.7015524687068223; -.6558964656854394; -.6077029271849502; -.5571583045146501;
         -.5044581449074642; -.4498063349740388; -.3934143118975651; -.3355002454194374;
         -.2762881937795320; -.2160072368760418; -.1548905899981459; -.9317470156008614e-1;
         -.3109833832718888e-1; .3109833832718888e-1; .9317470156008614e-1; .1548905899981459;
         .2160072368760418; .2762881937795320; .3355002454194374; .3934143118975651;
         .4498063349740388; .5044581449074642; .5571583045146501; .6077029271849502;
         .6558964656854394; .7015524687068223; .7444943022260685; .7845558329003993;
         .8215820708593359; .8554297694299461; .8859679795236130; .9130785566557919;
         .9366566189448779; .9566109552428079; .9728643851066921; .9853540840480059;
         .9940319694320907; .9988664044200711]
    z = [z ones(50)]
    p, PolyZeros(z)
end
