untitled3.asv

牛顿插值法

function ja=jacobi(A,b,x0)
D=diag(diag(A));
U=-triu(A,1);
L=-tril(A,-1);
C=inv(D);
B=L+U;