deff('[x,y,z]=sph(alp,tet)',['x=r*cos(alp).*cos(tet)+orig(1)*ones(tet)'; 'y=r*cos(alp).*sin(tet)+orig(2)*ones(tet)'; 'z=r*sin(alp)+orig(3)*ones(tet)']); r=1;orig=[0 0 0]; x=linspace(-%pi/2,%pi/2,40);y=linspace(0,%pi*2,20); x(5:8)=%inf*ones(5:8); x(30:35)=%inf*ones(30:35); [x1,y1,z1]=eval3dp(sph,x,y); plot3d1(x1,y1,z1)