解方程
第7章matlab解方程与函数极值71线性方程组求解72非线
D=diag(diag(A))。 %求 A的对角矩阵 L=tril(A,1)。 %求 A的下三角阵 U=triu(A,1)。 %求 A的上三角阵 B=D\(L+U)。 f=D\b。 y=B*x0+f。 n=1。 %迭代次数 while norm(yx0)=eps x0=y。 y=B*x0+f。 n=n+1。 end 例 75 用 Jacobi迭代法求解下列线性方程组。 设迭代初值为 0