: 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)!