function val = sinPhi(a,b)
if a==0
  val=1;
else
  val = -sign(a)*b/(sqrt(a^2+b^2));
end
  