defaultfont:="ptmr8r";
warningcheck:=0;
beginfig(1)
size=75;
path p[];
def pointtopair(expr x,y,z)=
(-x*cosd a + y*sind a,
-x*sind a * sind b -y*cosd a * sind b
+ z*cosd b)
enddef;
%
%basic path (the shape of the boundary)
%can be molded, can be constrained etc
p1:= (0,3size){right}..
{down}(1.1size,1.75size){down}..
(.35size,.75size)..(.175size,.375size)..
{left}origin;
%path with regular---nearly so---
%distributed points
n:=0;%number of points along the curve
p10:= point 0 of p1 hide(n:=n+1)..
for t:=1 upto 19: hide(n:=n+1)
point .05t of p1..endfor
point 1 of p1 hide(n:=n+1)..
for t:=1 upto 13: hide(n:=n+1)
point 1+t/14 of p1..endfor
point 2 of p1 hide(n:=n+1)..
for t:=1 upto 3: hide(n:=n+1)
point 2+t/4 of p1..endfor
point 3 of p1 hide(n:=n+1)..
for t:=1 upto 3: hide(n:=n+1)
point 3+t/4 of p1..endfor
origin;
%viewing angle parameters
b:=-10; a:=60;
p100:= for k=0 upto n-1:
pointtopair(0,xpart(point k of p10),
ypart(point k of p10))..
endfor pointtopair(0,0,0);
p200:= for k=0 upto n-1:
pointtopair(xpart(point k of p10), 0,
ypart(point k of p10))..
endfor pointtopair(0,0,0);
p300:= for k=0 upto n-1:
pointtopair(0,-xpart(point n-k of p10),
3size-ypart(point n-k of p10))..
endfor pointtopair(0,0,0);
p400:= for k=0 upto n-1:
pointtopair(-xpart(point n-k of p10),
0, 3size-ypart(point n-k of p10))..
endfor pointtopair(0,0,0);
%
fill (-1.5size,-size)--(-1.5size,5size)--
(1.5size,5size)--(1.5size,-size)--cycle;
%
drawoptions(withcolor white);
%
pickup pencircle scaled .5pt;
%Top ring and hang up (rope)
draw point 0 of p100..
point 0 of p100 + (0,.1size)..cycle;
draw point 0 of p100 + (0,.1size)..
point 0 of p100 + (0,3size);
%Draw boundary curves
draw p100; draw p200; draw p300; draw p400;
%
%Draw (partially hidden) regular surfaces
pickup pencircle scaled .1pt;
for k=0 step 1 until n:
draw point k of p200..point n-k of p300;
endfor
for k=0 upto n:
draw point k of p400..point n-k of p100;
endfor
fill p100..reverse p200..cycle withcolor black;
%Front
pickup pencircle scaled .1pt;
draw p100; draw p200;
draw point 0 of p100--origin;
%
%Draw regular surface which is in sight
for k=0 step 1 until n:
draw point k of p100..point n-k of p200;
endfor
endfig;
end;
