function x = sor(A,b,omega,x0,e,N)
% 用途：用分量形式的SOR迭代解線性方程組Ax=b
% 格式：x = sor(A,b,omega,x0,e,N)  A為系數矩陣，b為右端向量，x返回解向量，x0為
%     初值向量（默認原點），e為精度（默認1e-4），設置迭代次數上限以防發散（默認
%     500），omega是松弛因子，一般取1~2之間的數（默認1.5）
n = length(b);
if nargin<6,N = 500;end
if nargin<5,e = 1e-4;end
if nargin<4,x0 = zeros(n,1);end
if nargin<3,omega = 1.5;end
x = x0;x0 = x+2*e;
k = 0;L = tril(A,-1);U = triu(A,1);
while norm(x0-x,inf)>e&k<N,
    k = k+1,x0 = x;
    for i = 1:n
        x1(i) = (b(i)-L(i,1:i-1)*x(1:i-1,1)-U(i,i+1:n)*x0(i+1:n,...
            1))/A(i,i);
        x(i) = (1-omega)*x0(i)+omega*x1(i);
    end
    disp(x')
end
if k == N,warning('以達迭代次數上限');end