请问 高斯－赛德尔迭代 实例如何编程？？？

sunny819

刚接触matlab ，遇到“高斯－赛德尔迭代编程”，头大了！请高手指点！<BR>

happy

<P>function x=gaussdiedai(A,b,n)</P>
<P>%用高斯－赛德尔迭代解线性方程组</P>
<P>Len=length(A);</P>
<P>%求A的维数</P>
<P>U=-triu(A,1);B=A+U;</P>
<P>if (norm(B)==0)</P>
<P>x='Error!';  </P>
<P>%若B不可逆，给出错误信息，程序中止</P>
<P>else</P>
<P>     B2=inv(B)*U;f2=inv(B)*b;</P>
<P>    rou=max(abs(eig(B2)))</P>
<P>    %否则计算B2谱半径</P>
<P>    if (rou&gt;=1)</P>
<P>        x='Not convergent!'; </P>
<P>%谱半径&gt;=1，不收敛，中止程序</P>
<P>    else</P>
<P>        %否则进行n次高斯－赛德尔迭代，输出结果</P>
<P>        x=ones(Len,1);</P>
<P>        for i=1:n</P>
<P>            x=B2*x+f2;</P>
<P>        end</P>
<P>    end</P>
<P>end</P>

supernew

还有有没有其他的？