帮看看程序有什么错，怎么不能生成三维图

[那年伊始]

function xiachen
c=13*360;
m=3.1;q=0.89;a=8;h=304;v1=0.92;
w0=m*q*1000*cosd(a);r=h/2.25;
t1=30/0.92;
t(4)=330/v1;
p1=1-exp(-c.*t(4));
p2=1-exp(-c.*(t(4)-t1));
p3=1-exp(-c.*(t(4)-2.*t1));
p4=1-exp(-c.*(t(4)-3.*t1));
p5=1-exp(-c.*(t(4)-4.*t1));
p6=1-exp(-c.*(t(4)-5.*t1));
p7=1-exp(-c.*(t(4)-6.*t1));
p8=1-exp(-c.*(t(4)-7.*t1));
p9=1-exp(-c.*(t(4)-8.*t1));
p10=1-exp(-c.*(t(4)-9.*t1));
p11=1-exp(-c.*(t(4)-10.*t1));
x=-300:30:700;
y=-200:30:500;
[X,Y]=meshgrid(x,y);
Z=-w0./4.*(erf((sqrt(pi)./r).*Y)-erf((sqrt(pi)./r).*(Y-124))).*(p1.*(erf((sqrt(pi)./r).*X)-erf((sqrt(pi)./r)*(X-v1.*t1)))+...
p2.*(erf((sqrt(pi)./r).*(X-v1.*t1))-erf((sqrt(pi)./r).*(X-2*v1.*t1)))+...
p3.*(erf((sqrt(pi)./r).*(X-2*v1.*t1))-erf((sqrt(pi)./r)*(X-3*v1.*t1)))+...
p4.*(erf((sqrt(pi)./r).*(X-3.*v1.*t1))-erf((sqrt(pi)./r)*(X-4.*v1.*t1)))+...
p5.*(erf((sqrt(pi)./r).*(X-4.*v1.*t1))-erf((sqrt(pi)./r)*(X-5.*v1.*t1)))+...
p6.*(erf((sqrt(pi)./r).*(X-5.*v1.*t1))-erf((sqrt(pi)./r)*(X-6.*v1.*t1)))+...
p7.*(erf((sqrt(pi)./r).*(X-6.*v1.*t1))-erf((sqrt(pi)./r)*(X-7.*v1.*t1)))+...
p8.*(erf((sqrt(pi)./r).*(X-7.*v1.*t1))-erf((sqrt(pi)./r)*(X-8.*v1.*t1)))+...
p9.*(erf((sqrt(pi)./r).*(X-8.*v1.*t1))-erf((sqrt(pi)./r)*(X-9.*v1.*t1)))+...
p10.*(erf((sqrt(pi)./r).*(X-9.*v1.*t1))-erf((sqrt(pi)./r)*(X-10.*v1.*t1)))+...
p11.*(erf((sqrt(pi)./r).*(X-10.*v1.*t1))-erf((sqrt(pi)./r)*(X-11.*v1.*t1))));;
surf(X,X,Z);
shading interp;
xlabel('走向方向(m)');
ylabel('倾向方向(m)');
zlabel('下沉(mm)');

[那年伊始]

各位大虾，怎么没人回答啊