3
kp=xget("pixmap");xset("pixmap",1);
4
xset("wwpc"); // clean pixmap
6
//first plot, to fix boundaries
7
plot3d1(t,t,sin(t)'*cos(t),35,45," ",[1,2,4]);
8
xset("wshow"); // show pixmap
10
if driver()=='Pos' then st=4;else st=2;end;
11
for i=35:st:80, // loop on theta angle
13
plot3d1(t,t,sin(t)'*cos(t),i,45," ",[1,0,4])
16
for i=45:st:80, //loop on alpha angle
18
plot3d1(t,t,sin(t)'*cos(t),80,i," ",[1,0,4])