21.14 scherk.input
Scherk's minimal surface, defined by:
Scherk's minimal surface
.
See: A Comprehensive Introduction to Differential Geometry, Vol. 3,
by Michael Spivak, Publish Or Perish, Berkeley, 1979, pp. 249-252.
(xOffset, yOffset):DFLOAT Offsets for a single piece
of Scherk's minimal surface
drawScherk(m,n) == Draw Scherk's minimal surface
free xOffset, yOffset on an by patch
space := createThreeSpace()
for i in 0..m-1 repeat
xOffset := i*%pi
for j in 0 .. n-1 repeat
rem(i+j, 2) = 0 => 'iter Draw only odd patches
yOffset := j*%pi
drawOneScherk(space) Draw a patch
makeViewport3D(space, "Scherk's Minimal Surface")
scherk1(u,v) == The first patch that makes
x := cos(u)/exp(v) up a single piece of
point [xOffset + acos(x), yOffset + u, v, abs(v)] Scherk's minimal surface
scherk2(u,v) == The second patch
x := cos(u)/exp(v)
point [xOffset - acos(x), yOffset + u, v, abs(v)]
scherk3(u,v) == The third patch
x := exp(v) * cos(u)
point [xOffset + u, yOffset + acos(x), v, abs(v)]
scherk4(u,v) == The fourth patch
x := exp(v) * cos(u)
point [xOffset + u, yOffset - acos(x), v, abs(v)]
drawOneScherk(s) == Draw the surface by
breaking it into four
patches and then drawing
the patches
makeObject(scherk1,-%pi/2..%pi/2,0..%pi/2,space==s,
var1Steps == 28, var2Steps == 28)
makeObject(scherk2,-%pi/2..%pi/2,0..%pi/2,space==s,
var1Steps == 28, var2Steps == 28)
makeObject(scherk3,-%pi/2..%pi/2,-%pi/2..0,space==s,
var1Steps == 28, var2Steps == 28)
makeObject(scherk4,-%pi/2..%pi/2,-%pi/2..0,space==s,
var1Steps == 28, var2Steps == 28)
void()