1
function fac3d1(x,y,z,T,A,leg,flags,ebox)
5
write(%io(2),'WARNING : fac3d1 is obsolete see plot3d plot3d1 instead');
7
s_mat=['deff(''[x,y,z]=sph(alp,tet)'',[''x=r*cos(alp).*cos(tet)+orig(1)*ones(tet)'';';
8
'''y=r*cos(alp).*sin(tet)+orig(2)*ones(tet)'';';
9
'''z=r*sin(alp)+orig(3)*ones(tet)'']);';
11
'[x1,y1,z1]=eval3dp(sph,linspace(-%pi/2,%pi/2,40),linspace(0,%pi*2,20));';
13
'r=1/2;orig=[-1 0 0];';
14
'[x2,y2,z2]=eval3dp(sph,linspace(-%pi/2,%pi/2,40),linspace(0,%pi*2,20));';
16
'x=[x1 x2];y=[y1 y2];z=[z1 z2];';
18
write(%io(2),s_mat);execstr(s_mat);
21
case 3 then plot3d1(x,y,z)
22
case 4 then plot3d1(x,y,z,T)
23
case 5 then plot3d1(x,y,z,T,A)
24
case 6 then plot3d1(x,y,z,T,A,leg)
25
case 7 then plot3d1(x,y,z,T,A,leg,flags)
26
case 8 then plot3d1(x,y,z,T,A,leg,flags,ebox)