%!

% x -> arcsin(x)

/asin { dup dup      % x x x
mul 1                   % x x^2 1
exch sub sqrt       % x sqrt(1-x^2)
atan
} def

% -----------------------------------------------------------

0.5 asin ==