flexin.m

matlab程序

function x=Flexin(A,b,x,N,w)
temp1=0;
temp2=0;
temp3=0;
temp4=0; 
for k=1:N
    for j=2:4
      temp1=temp1+A(1,j)*x(j);
    end
    x(1)=(1-w)*x(1)+w*(b(1)-temp1)/A(1,1);
    temp1=0;
    for j=3:4
       temp2=temp2+A(2,j)*x(j);
    end
      x(2)=(1-w)*x(1)+w*(b(2)-A(2,1)*x(1)-temp2)/A(2,2);
      temp2=0;
    for j=1:2
       temp3=temp3+A(3,j)*x(j);
    end
       x(3)=(1-w)*x(1)+w*(b(3)-temp3-A(3,4)*x(4))/A(3,3);
       temp3=0;
    for j=1:3
       temp4=temp4+A(4,j)*x(j);
    end
       x(4)=(1-w)+w*(b(4)-temp4)/A(4,4);
       temp4=0;
end