用matlab软件作锥面x^2+y^2=z^2和柱面(x-1)^2+y^2=1相交的图形，求大神帮忙写一下程序~~~

你看看程序，不懂得可以问我。

%程序

xgrid = -2:0.1:2;

ygrid = -2:0.1:2;

[x,y] = meshgrid(xgrid,ygrid);

z_p = sqrt(x.^2+y.^2);

z_n = -sqrt(x.^2+y.^2);

mesh(x,y,z_p);%锥面上半部分

hold on;

mesh(x,y,z_n);%锥面下半部分

hold on;

t=0:pi/20:2*pi;

x1=sin(t)+1;

y1=cos(t);

z1=linspace(-3,3,length(t));

X=meshgrid(x1);Y=meshgrid(y1);Z=[meshgrid(z1)]';

mesh(X,Y,Z)%圆柱面

xlabel('x')

ylabel('y')

zlabel('z')

axis equal