Surface with holes 1

t=linspace(-%pi,%pi,40);z=sin(t)'*cos(t);
z1=find(abs(z) > 0.5);
z(z1)=%inf*z1;
plot3d1(t,t,z);