用最小二乘法解回归方程出错，请大侠帮忙看看

[smartpig]

回归模式是这样的：
1856
简化成用最小二乘法解的式子是这样的：
1857
下面这个程序不知道是什么错，请各位大侠帮忙看看

function nafit(x,y,z,m,n)
A=zeros(m,m);
for r=1:m
c=0;
for i=0:n
for j=0:i
c=c+1;
A(r,c)=x.^(i-j).*y.^j;
end
end
end
a=A\z;
p=fliplr(a');

谢谢啦！！！

[smartpig]

输入下面的参数值：
x=[1000
1200
1400
1600
1800
2000
2200
2400
2600
2800
];
>> y=[354.6
362.8
356.6
351
345.2
343.6
346.6
341.6
335.8
332
];
>> z=[11.12
12.86
14.78
16.28
17.81
19.85
21.71
23.54
24.98
28.22
];
>> z=[288.2
283.4
284.2
284.6
284.6
290
288.2
292.8
306
310.2
]
m=10;n=3;
显示出现
Subscripted assignment dimension mismatch.
Error in ==> nafit at 8
A(r,c)=x.^(i-j).*y.^j;
的错误，请高手帮忙看看了，谢谢了！！