text7.m

来自「fit programme for a serials of data」· M 代码 · 共 12 行

%text work
A=[4 -1 0 -1 0 0;-1 4 -1 0 -1 0;0 -1 4 -1 0 -1;-1 0 -1 4 -1 0;0 1 0 1 4 -1;0 0 -1 0 -1 4]
b=[0 5 -2 5 -2 6];
X0=[1;2;3;6;5;4];
disp(sprintf('jacob iteration'));
[X,n]=jacobi(A,b,X0,0.0001,10000)
disp(sprintf('Gao-Sai iteration'));
[X,n]=gaosai(A,b,X0,0.0001,10000)
disp(sprintf('chaosong iteration'));
[X,n]=chaosong(A,b,X0,1.334,0.0001,10000)
[X,n]=chaosong(A,b,X0,1.95,0.0001,10000)
[X,n]=chaosong(A,b,X0,0.95,0.0001,10000)