1set term epslatex color solid
2set output 'eg7.tex'
3unset border
4set dummy u,v
5unset key
6unset encoding
7set parametric
8set view 60, 30, 1.2, 1.3
9set isosamples 50, 20
10set hidden3d offset 1 trianglepattern 3 undefined 1 altdiagonal bentover
11set noxtics
12set noytics
13set noztics
14set title "Interlocking Tori - PM3D surface with depth sorting"
15set urange [ -3.14159 : 3.14159 ] noreverse nowriteback
16set vrange [ -3.14159 : 3.14159 ] noreverse nowriteback
17set pm3d depthorder
18x=.10; y=.15; dy=.05
19set label "left torus:"                           at screen x,y; y=y-dy
20set label '$x=\cos u+\frac{1}{2}\cos u \cos v$'   at screen x,y; y=y-dy
21set label '$y=\sin u+\frac{1}{2}\sin u \cos v$'   at screen x,y; y=y-dy
22set label "$z=\\frac{1}{2}\\sin v$"               at screen x,y
23x=.65; y=.08
24set label '\parbox{2.5in}{right torus:\\$x=1+\cos u+\fr\
25ac{1}{2}\cos u \cos v$\\$y=\frac{1}{2}\sin v$\\\
26$z=\sin u + \frac{1}{2}\sin u \cos v$}' at screen x,y left
27set pm3d interpolate 1,1 flush begin noftriangles nohidden3d corners2color mean
28splot cos(u)+.5*cos(u)*cos(v),sin(u)+.5*sin(u)*cos(v),.5*sin(v) with pm3d, \
29    1+cos(u)+.5*cos(u)*cos(v),.5*sin(v),sin(u)+.5*sin(u)*cos(v) with pm3d
30
31