*

: pi 3.1416926 ; : pi/2 pi 2 / ; : pi/4 pi 4 / ;
: n .5 - 2 * ; : ^2 dup * ; : inv -1 pow ;
: sinc dup >r sin r> / ;
: (x<1)! 1 < * ;
: (x<0.95)! 0.95 / 1 < * ;
: (x/(x+1)) dup 1 + / ;
: (1/(1-x)) 1 - negate inv ;
: (x^2/(x^2+1)) ^2 (x/(x+1)) ;
: x x n ; : y y n ; : e 0.001 ;
: r ( x y - r ) ^2 >r ^2 r> + sqrt ;
: a ( y x - a ) atan2 pi/2 / ;
: m 1 mod abs 0.2 < ;
: rr x y r (1/(1-x)) sqrt sqrt 10 * m ;
: aa y x a 8 * m ;
: x' rr ;
: y' aa ;
x' y' or x y r 0.98 > or x y r (x<1)!