: x x .5
- 2
* ;
: y y .5
- 2
* ;
: ^2
dup * ;
: z^2
over over z* ;
: a^2+b^2
( a b -> a^2+b^2 ) ^2
swap ^2
+ ;
: a/r,b/r
( a b r -> a/r b/r )
dup >r / swap r> / swap ;
: a-b,a+b
( a b -> a-b a+b )
over over + >r - r> ;
: 1/z
( a b -> [a-b]/r [a+b]/r, r=a^2+b^2 )
over over a^2+b^2
>r a-b,a+b
r> a/r,b/r
;
: zdup
over over ;
: z/ 1/z
z* ;
: riman
( z -> 1/[1+z^2] ) z^2 1 0
z+ 1/z
;
x y riman
t 7
* sin t 3
* sin z*