1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
|
.PS
# worm.m4
gen_init
define(`svgcolor',`sprintf("rgb(%g,%g,%g)",`$1',`$2',`$3')')
`define case { exec sprintf("$%g",floor($1+0.5)+1); }'
`define hsvtorgb {
h_ = $1; s_ = $2; v_ = $3
f_ = h_/60-floor(h_/60)
p_ = v_*(1-s_)
q_ = v_*(1-f_*s_)
t_ = v_*(1-(1-f_)*s_)
case(pmod(floor(h_/60),6)+1,
$4 = v_; $5 = t_; $6 = p_,
$4 = q_; $5 = v_; $6 = p_,
$4 = p_; $5 = v_; $6 = t_,
$4 = p_; $5 = q_; $6 = v_,
$4 = t_; $5 = p_; $6 = v_,
$4 = v_; $5 = p_; $6 = q_)
$4 = int(($4)*255+0.5)
$5 = int(($5)*255+0.5)
$6 = int(($6)*255+0.5)
}'
Worm: [
linethick_(0.1)
dimen = 1.25 # scale parameter; try dimen=2 for a larger picture
w = 0.15 * dimen # worm radius
nhues = 265
S: dimen,0
for inx = 1 to nhues +2 do {
theta = inx*360/nhues
R: dimen*cosd(3*theta), dimen*sind(2*theta)
D: R-(S.x,S.y)
t = sqrt(max(0,w^2/(D.x^2+D.y^2)-1/4))
Ra: 0.5<S,R> + (t*D.y,-t*D.x)
Rb: 0.5<S,R> - (t*D.y,-t*D.x)
if inx <= nhues then {
hsvtorgb(theta+120,1,1,r,g,b)
rgbfill( r,g,b,
arc cw from Ra to Rb rad w with .c at S
arc ccw from Rb to Ra rad w with .c at R )
} else {
arc cw from Ra to Rb rad w with .c at S
arc ccw from Rb to Ra rad w with .c at R }
S: R
}
]
"svg_fsize(R,120)" at 0.9<Worm.se,Worm.ne> ljust
"svg_fsize(G,120)" at Worm.e ljust
"svg_fsize(B,120)" at 0.1<Worm.se,Worm.ne> ljust
"svg_fsize(Y,120)" at 0.1<Worm.sw,Worm.nw> rjust
"svg_fsize(M,120)" wid 0.2 at Worm.w rjust
"svg_fsize(C,120)" at 0.9<Worm.sw,Worm.nw> rjust
box invis ht Worm.ht wid Worm.wid+0.4 at Worm
[ thicklines_
C: circle rad 1.5 at 0,0
smallrad = C.rad/2
for t=0 to 330 by 30 do {
line from Rect_(C.rad-0.05,90-t) to Rect_(C.rad+0.05,90-t)
sprintf("%g",t) at Rect_(C.rad+0.2,90-t)
}
"svg_fsize(R,120)" at Rect_(C.rad+0.45,90-0)
"svg_fsize(Y,120)" wid 0.2 at Rect_(C.rad+0.45,90-60)
"svg_fsize(G,120)" at Rect_(C.rad+0.45,90-120)
"svg_fsize(C,120)" at Rect_(C.rad+0.45,90-180)
"svg_fsize(B,120)" at Rect_(C.rad+0.45,90-240)
"svg_fsize(M,120)" at Rect_(C.rad+0.45,90-300)
linethick_(0)
for t = 1 to 90 do {
hsvtorgb(t/90*360,1,1,r,g,b)
line shaded svgcolor(r,g,b) outlined svgcolor(r,g,b) \
from C to Rect_(C.rad-0.05,90-(t-0.5)/90*360) \
then to Rect_(C.rad-0.05,90-(t+0.5)/90*360) \
then to C
}
R: circle diam C.rad at Rect_(C.rad*3/8,90) color "rgb(255,0,0)"
G: circle diam C.rad at Rect_(C.rad*3/8,-30) color "rgb(0,255,0)"
B: circle diam C.rad at Rect_(C.rad*3/8,-150) color "rgb(0,0,255)"
RG: cintersect(R,R.rad,G,G.rad)
GR: cintersect(R,R.rad,G,G.rad,R)
GB: cintersect(G,G.rad,B,B.rad)
BG: cintersect(G,G.rad,B,B.rad,R)
BR: cintersect(B,B.rad,R,R.rad)
RB: cintersect(B,B.rad,R,R.rad,R)
# The following could be done other ways depending on the postprocessor
`define lozenge {
arc from $1$2 to $2$1 with .c at $2
arc from $2$1 to $1$2 with .c at $1 }'
# yellow
linethick = 0.1; rgbfill(255,255,0,lozenge(R,G))
linethick = 0.4; rgbdraw(255,255,0,lozenge(R,G))
# cyan
linethick = 0.1; rgbfill(0,255,255,lozenge(G,B))
linethick = 0.4; rgbdraw(0,255,255,lozenge(G,B))
# magenta
linethick = 0.1; rgbfill(255,0,255,lozenge(B,R))
linethick = 0.4; rgbdraw(255,0,255,lozenge(B,R))
`define white {
arc cw from RB to GR with .c at R
arc cw from GR to BG with .c at G
arc cw from BG to RB with .c at B }'
linethick = 0.1; rgbfill(255,255,255,white)
linethick = 0.4; rgbdraw(255,255,255,white)
setrgb(255,255,255)
line from RB to GR then to BG then to RB
resetrgb
] with .w at Worm.e+(0.5,0)
.PE
|