matlab不连续分段函数,分段函数画出的曲面不连续如何办？

代码如下：

a0=0.1025;a1=-0.8403;a2=-0.0910;b1=-0.4507;b2=-0.1018;

[t,xi]=meshgrid(-pi/3:0.1:pi/3,-8:0.1:0.1025*3/sqrt(3));

sigmam=(1/3)*sqrt(3)*xi;

rhot=(1/(2*a2))*((-a1)-sqrt((a1)^2-4*a2*(a0-sigmam)));

rhoc=(1/(2*b2))*((-b1)-sqrt((b1)^2-4*b2*(a0-sigmam)));

u=2*rhoc.*(rhoc.^2-rhot.^2).*cos(t)+rhoc.*(2*rhot-rhoc).*sqrt(4*(rhoc.^2-rhot.^2).*(cos(t)).^2+5*rhot.^2-4*rhot.*rhoc);

v=4*(rhoc.^2-rhot.^2).*(cos(t)).^2+(rhoc-2*rhot).^2;

rho=u./v;

x=((1/sqrt(3))*xi+sqrt(2/3)*rho.*cos(t));

y=((1/sqrt(3))*xi+sqrt(2/3)*rho.*cos(t-2*pi/3));

z=((1/sqrt(3))*xi+sqrt(2/3)*rho.*cos(t+2*pi/3));

mesh(x,y,z)

hold on

a0=0.1025;a1=-0.8403;a2=-0.0910;b1=-0.4507;b2=-0.1018;

[t,xi]=meshgrid(pi/3:0.1:pi,-8:0.1:0.1025*3/sqrt(3));