Vectorize the following MATLAb code - run the vectorized code:

r=0:1/14:1;th=0:2*pi/22:2*pi;

for j=1:15

for i =1:23

x(i,j)=r(j)*sin(th(i));

y(i,j)=r(j)*cos(th(i));

z(i,j)=sqrt (abs(1-x(i,j)^2-y(i,j)^2));

end

end

surf (x,y,z)