首先提一点建议:function[]=J()因为没有输入输出,所以就没必要弄成一个函数文件,弄成一个脚本文件就可以了
我调试了下,错误出在xx=(b-A*xold+D.*xold)./D;
我把处理这条程序前所有变量的值给你列一下,你就知道为什么错了
K>> whos
Name Size Bytes Class Attributes
A 99x100 79200 double
D 99x1 792 double
b 99x1 792 double
i 1x1 8 double
itr 1x1 8 double
m 1x1 8 double
xold 99x1 792 double
xx 99x1 792 double
A 99*100 xold 99*1,他两无法相乘