*

: 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*

Based on: Hypno4 - DarkstarAG