[求助]Jacobi迭代问题

yangdong28

程序如下：

1. A=[9 -4;3 -10];
   
2. B=[5 -7]';
   
3. B=eye(2)-diag(1/9,-1/10)*A;
   
4. g=diag(1/9,-1/10)*b;
   
5. x=[1 0]';
   
6. error=Bx+g;
   
7. while error>e-4
   
8. x=B*x+g-x
   
9. error=B*x+g-x
   
10. end

复制代码

结果输出x＝[10]',似乎没有进行迭代运算，请问问题出在哪里？
谢谢！

happy

参考下面帖子里的程序吧
http://forum.vibunion.com/thread-1561-1-1.html