procedure sarcsin(var x,y: scalar);
var
  x0,x1,x2: scalar;
begin { sarcsin }
  sprod(x,x,x2);
  sdiff(sone,x2,x0);
  ssqrt(x0,x1);
  squot(x,x1,x0);
  sarctan(x0,y);
end { sarcsin };

