function y = pp(s)
global Q

     y = tan(s) -sqrt(Q-s.^2)./s ;