function y = qq(s)

  global Q

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