y0 =0.000578428927106762 y1 =0.998184239328006 sin(b)=-2.44929359829471E-16 cos(b) =1 npoints=51