編程貓拋物線怎麼做

x=linspace(-10,10,100)';
y=x.^2;
z=y;
n=100;
t=(0:n)/n*2*pi;
X=x*ones(1,n+1);
Y=y*cos(t);
Z=z*sin(t);
surf(X,Y,Z)