7
xsetech([0.,0.,0.5,0.5],[-1,1,-1,1]);
12
xsetech([0.5,0.,0.5,0.5],[-1,1,-1,1]);
13
h1=((1:ncolor)'/ncolor)*[1 0 0];
17
xsetech([0.,0.5,0.5,0.5],[-1,1,-1,1]);
18
h2=graycolormap(ncolor);
22
xsetech([0.5,0.5,0.5,0.5],[-1,1,-1,1]);
23
deff('[x,y,z]=knot(u,v)',['vv=ones(v)';'uu=ones(u);';
24
'x=(5.*cos(u)+cos(u).*cos(v))';
25
'y=(5.*sin(u)+sin(u).*cos(v))';
27
nx=60;Nx=(0:nx)/nx ;ny=20;Ny=(0:ny)/ny;
28
[xx,yy,zz]=eval3dp(knot,2*%pi*Nx,2*%pi*Ny);
29
XXX=[-xx xx];YYY=[-yy zz];ZZZ=[-zz yy];
33
hh=hotcolormap(ncolor);
35
kk1=modulo(kk1,60);kk2=kk1;
36
KKK=list(ZZZ,[kk1 kk2]);
37
plot3d(XXX,YYY,KKK,35,70," @ @ ",[2,1,4],[-6,6,-6,6,-6,6]);
40
// Change the window colormap