zz = 0.02
u = (xx - sw/2)*zz
v = (sh/2 - yy)*zz
'cmul uu, vv, u - (mx - sw/2)*0.01, v - (sh/2 - my)*0.01, u + 0.707, v - 0.707
'cdiv x, y, uu, vv, u - 1, v
sinz uu, vv, u, v
cdiv x, y, uu, vv, u - (mx - sw/2)*zz, v - (sh/2 - my)*zz
a
= (pi
+ _atan2(y
, x
))/(2*pi
)
'pset (xx, yy), hrgb&(a, m)
pset (xx
, yy
), hrgb&
(a
, log(5000*m
+ 1))
r
= 0.5 - 0.5*sin(2*pi
*h
- pi
/2) g
= (0.5 + 0.5*sin(2*pi
*h
*1.5 - pi
/2)) * -(h
< 0.66) b
= (0.5 + 0.5*sin(2*pi
*h
*1.5 + pi
/2)) * -(h
> 0.33)
n = 16
rr = 255*r - 0.15*mm - 35*p
gg = 255*g - 0.15*mm - 35*p
bb = 255*b - 0.15*mm - 35*p
x = xx
y = yy
x = xx
y = yy
'u + iv = (x + iy)*(a + ib)
sub cmul
(u
, v
, xx
, yy
, aa
, bb
) x = xx
y = yy
a = aa
b = bb
u = x*a - y*b
v = x*b + y*a
'u + iv = (x + iy)/(a + ib)
sub cdiv
(u
, v
, xx
, yy
, aa
, bb
) x = xx
y = yy
a = aa
b = bb
d = a*a + b*b
u = (x*a + y*b)/d
v = (y*a - x*b)/d