我要画一个分段函数的图形,在运行的时候老是出问题,求哪位大仙帮忙,叩谢!clearall;E=30;v=0.2;kappa=3-4*v;mju=E/(2*(1+v));E11=80;E22=12;E33=20;v12=0.28;v13=0.24;v23=0.08;G12=6;v21=(E22/E11)*v12;v31=(E33/E11)*v13;v32=(E33/E2

　　我要画一个分段函数的图形,在运行的时候老是出问题,求哪位大仙帮忙,叩谢!

　　clearall;

　　E=30;

　　v=0.2;

　　kappa=3-4*v;

　　mju=E/(2*(1+v));

　　E11=80;

　　E22=12;

　　E33=20;

　　v12=0.28;

　　v13=0.24;

　　v23=0.08;

　　G12=6;

　　v21=(E22/E11)*v12;

　　v31=(E33/E11)*v13;

　　v32=(E33/E22)*v23;

　　w=(1-v13*v31)*(1-v23*v32)-(v21+v23*v31)*(v12+v32*v13);

　　c11=E11*(1-v23*v32)/w;

　　c12=E11*(v21+v23*v31)/w;

　　c22=E22*(1-v13*v31)/w;

　　c66=G12;

　　a=c22*c66;

　　b=c11*c22-c12*c12-2*c12*c66;

　　c=c11*c66;

　　s1=sqrt((b+sqrt(b^2-4*a*c))/(2*a));

　　s2=sqrt((b-sqrt(b^2-4*a*c))/(2*a));

　　alpha1=(c11-c66*s1^2)/(c12*s1+c66*s1);

　　alpha2=(c11-c66*s2^2)/(c12*s2+c66*s2);

　　lambda=0.2072;

　　a1=2*mju*(lambda+1)*(s1^lambda);

　　a2=2*mju*(lambda+1)*(s2^lambda);

　　b1=-(lambda+1)*c66*(1+alpha1/s1)*(s1^lambda);

　　b2=-(lambda+1)*c66*(1+alpha2/s2)*(s2^lambda);

　　t=tan(lambda*pi/2);

　　g1=((lambda+1)*a1*alpha1+(lambda-kappa)*b1*s1-(lambda-1)*a1*t^2-(lambda+kappa)*b1*t^2)/((1+kappa)*(1+t^2));

　　g2=((lambda+1)*a2*alpha2+(lambda-kappa)*b2*s2-(lambda-1)*a2*t^2-(lambda+kappa)*b2*t^2)/((1+kappa)*(1+t^2));

　　l1=(a1*alpha1+b1*s1-(a1+b1)*t^2)/((1+kappa)*(1+t^2));

　　l2=(a1*alpha2+b2*s2-(a2+b2)*t^2)/((1+kappa)*(1+t^2));

　　x1=(-((lambda+1)*alpha1*a1+(lambda-1)*a1+(lambda-kappa)*s1*b1+(lambda+kappa)*b1))/((1+kappa)*(t+1/t));

　　x2=(-((lambda+1)*alpha2*a2+(lambda-1)*a2+(lambda-kappa)*s2*b2+(lambda+kappa)*b2))/((1+kappa)*(t+1/t));

　　y1=(-(a1*alpha1+a1+b1*s1+b1))/((1+kappa)*(t+1/t));

　　y2=(-(a2*alpha2+a2+b2*s2+b2))/((1+kappa)*(t+1/t));

　　k1=(g1+(lambda+1)*l1)*sin(lambda*pi)+(x1+(lambda+1)*y1)*cos(lambda*pi);

　　k2=(g2+(lambda+1)*l2)*sin(lambda*pi)+(x2+(lambda+1)*y2)*cos(lambda*pi);

　　k3=-(x1+(lambda-1)*y1)*sin(lambda*pi)+(g1+(lambda-1)*l1)*cos(lambda*pi);

　　k4=-