function lgd100()
#
#  test polynomial suggested by Goedecker
#   Legendre polynomial of degree 100
#
    p = lgd(100);
    z = [-.9997137267734412; -.9984919506395958; -.9962951347331251; -.9931249370374435;
         -.9889843952429917; -.9838775407060570; -.9778093584869183; -.9707857757637063;
         -.9628136542558155; -.9539007829254917; -.9440558701362560; -.9332885350430795;
         -.9216092981453340; -.9090295709825297; -.8955616449707270; -.8812186793850184;
         -.8660146884971646; -.8499645278795913; -.8330838798884008; -.8153892383391763;
         -.7968978923903145; -.7776279096494955; -.7575981185197072; -.7368280898020207;
         -.7153381175730564; -.6931491993558020; -.6702830156031410; -.6467619085141293;
         -.6226088602037078; -.5978474702471787; -.5725019326213812; -.5465970120650942;
         -.5201580198817631; -.4932107892081909; -.4657816497733580; -.4378974021720315;
         -.4095852916783015; -.3808729816246300; -.3517885263724217; -.3223603439005292;
         -.2926171880384720; -.2625881203715035; -.2323024818449740; -.2017898640957360;
         -.1710800805386033; -.1402031372361140; -.1091892035800611; -.7806858281343664e-1;
         -.4687168242159163e-1; -.1562898442154308e-1; .1562898442154308e-1; .4687168242159163e-1;
         .7806858281343664e-1; .1091892035800611; .1402031372361140; .1710800805386033;
         .2017898640957360; .2323024818449740; .2625881203715035; .2926171880384720;
         .3223603439005292; .3517885263724217; .3808729816246300; .4095852916783015;
         .4378974021720315; .4657816497733580; .4932107892081909; .5201580198817631;
         .5465970120650942; .5725019326213812; .5978474702471787; .6226088602037078;
         .6467619085141293; .6702830156031410; .6931491993558020; .7153381175730564;
         .7368280898020207; .7575981185197072; .7776279096494955; .7968978923903145;
         .8153892383391763; .8330838798884008; .8499645278795913; .8660146884971646;
         .8812186793850184; .8955616449707270; .9090295709825297; .9216092981453340;
         .9332885350430795; .9440558701362560; .9539007829254917; .9628136542558155;
         .9707857757637063; .9778093584869183; .9838775407060570; .9889843952429917;
         .9931249370374435; .9962951347331251; .9984919506395958; .9997137267734412]
    z = [z ones(100)]
    # Legendre polynomial, roots are sensitive, see plot.
    
    p, PolyZeros(z)
end
