U x is x%sqrt(dot x x)

S r s is if 0>d then inf else if 0>t then inf else if 0<u then u else t
 where (d is sqr(s@1)+sqr b-dot v v)
 where (b is dot v(r@1))
 where (v is (*s)-*r)
 where (t is b+e)
 where (e is sqrt d)
 where u is b-e

I r h o is if ~type(*o) then f(S r(*o)) else g(S r o)
 where (f x is if ~x<*h then h else h(I r)/o@1)
 where g x is if ~x<*h then h else [x;U((*r)-(*o)-x*r@1)]

T r o d z l is if inf=*h then 0. else if ~0>g then 0. else if inf=i then -g else 0.
 where (g is dot(h@1)l)
 where (h is I r z o)
 where (i is *I j z o)
 where j is [(*r)+((r@1)**h)+d*h@1;-l]

N n o i is 0 f/+vs 4(!16)
 where f x y is x+T[0 0 -4.;U((i+(y%4)-n%2),n)]o(sqrt(2^-42))[inf;0 0 0.](U -1 -3 2.)

R k n is _0.5+15.9375*(N(n*1.)(C k 1.(0 -1 0.)))',f(vs n(!n*n)) 
 where f x is+[|*x;*|x]

C k r c is if 1=k then [c;r] else [[c;r*3];[[c;r]],(C(k-1)(r%2))'[+c+(-3 3(vs 2(2 3 6 7)))*r%sqrt 12]] 

T r o d z l
 where (r is [0 0 -4;-0.4364358 0.2182179 0.8728716])
 where (o is C 3(1.)(0 -1 0.))
 where (d is sqrt(2^-42)
 where (z is [inf;0 0 0.])
 where (l is -0.2672612 -0.8017837 0.5345225)


\

T 9 where T a is h+a where (g is a) where h is I 0 where (I b is g+b where g is 0)
T 9 where T a is h+a where (g is a) where h is I 0 where (I b is q+b where q is 0)
S(1)2 where S r s is d where (d is s+b+v) where (b is v+r) where (v is s-r)
S(1)2 where S r s is d s r where (d s r is s+b s r+v s r) where (b s r is v s r+r) where (v s r is s-r)
f 3 where f x is g where (g is h+1) where h is x+2