[求助]请帮我看看下面的程序哪里错了，我是初学者，还不懂！谢谢啦！
 

a=sqrt(123.6/239);
r=1.0184;
x=[0.93174:4.34471];
w=1/6/(1/a)^(2/3)*(6*x^(5/2)*(1/a)^(5/3)/(x^3/a^2)^(1/2)*atan((x^3/a^2)^(1/2))-...
3/2*x^(5/2)*(1/a)^(5/3)/(x^3/a^2)^(5/6)*3^(1/2)*log(1-3^(1/2)*(x^3/a^2)^(1/6)+...
(x^3/a^2)^(1/3))-3*x^(5/2)*(1/a)^(5/3)/(x^3/a^2)^(5/6)*atan(1/2*(x^3/a^2)^(1/6)/...
(1-1/2*3^(1/2)*(x^3/a^2)^(1/6)))-6*x^(5/2)*(1/a)^(5/3)/(x^3/a^2)^(5/6)*...
atan((x^3/a^2)^(1/6))+3/2*x^(5/2)*(1/a)^(5/3)/(x^3/a^2)^(5/6)*3^(1/2)*log(1+3^(1/2)*(x^3/a^2)^(1/6)+...
(x^3/a^2)^(1/3))-3*x^(5/2)*(1/a)^(5/3)/(x^3/a^2)^(5/6)*atan(1/2*(x^3/a^2)^(1/6)/(1+1/2*3^(1/2)*(x^3/a^2)^(1/6))));
h=(2/5*x^(5/2)-x)/(3*a.^3);
f=(2/3*x^(3/2)-x)/a;
k=x.*atan(1/a);
y=1./(1+a.^2)*(x/r-4*a^3/r^2.*(h-f+w-k)) ;
z=2*a^3*r*((x^1.5-1)/(3*a^3)-(x^0.5-1)/a+atan(x^0.5/a)-atan(1./a));
subplot(2,2,1);
plot(x,y)
subplot(2,2,2);
plot(x,z)
subplot(2,2,3);
plot(y,z)
主要是想利用x这个参数，来画y与z的关系图，但是不知道哪里错了。

它总是出现??? Error using ==> mpower
atrix must be square.是什么意思啊？
在这里万分感谢咯！

mpower的距阵必须是方阵

我也知道要是矩阵，但是不知道要怎么改，可以教我吗，谢谢你了！

输出的应该是x y z的关系,不是仅仅是y z了

但是我还是有点想不通，比如：
t=[0:pi/100:2*pi];
x=sin(t);
y=cos(t);
plot(x,y)
却能够画出关于x，y的图像，为什么那个就不能呢？

用plot3(x,y,z)输出