SOR超松弛迭代-求三角对阵

此处是我的计算方法的作业
K为迭代次数

%---逐次超松弛迭代法-----
%---successive over-reaxation iteration method
clear;clc;
A=zeros (1000,1000);
b=zeros(1000,1);
for i=1:1000
    for j=1:1000
        if(i==j) 
            A(i,j)=8;
        else if(i==j-1)
                A(i,j)=1;
            else if(i==j+1)
                    A(i,j)=6;
                end
            end
        end
    end
end
b(1)=9;
b(1000)=14;
for i=2:999
    b(i)=15;
end
N=length(b); %解向量的维数
fprintf('库函数计算结果：');
x=inv(A)*b %库函数计算结果
x=zeros(N,1);%迭代初始值
%-----(A=D-E-F)------
D=diag(diag(A));
E=-tril(A,-1);%下三角
F=-triu(A,1);%上三角
w=1.1; %松弛因子，一般0<w<2
B=inv(D-w*E)*[(1-w)*D+w*F];g=w*inv(D-w*E)*b;
eps=0.00001;%相邻解的距离小于该数时，结束迭代
%--------开始迭代-------
for k=1:1000 %最大迭代次数为1000
fprintf('第%d次迭代：',k);
y=B*x+g;
if abs(x-y)<eps
break;
end
x=y
end
x