jacobi.m

给出了线性方程组的CG法、Gauss_Seidel、SOR、Jacobi、fem求解

clear all
clc

format long

A=[10 3 1;2 -10 3;1 3 10];
b=[14;-5;14];
x=zeros(size(b));
y=zeros(size(b));
n=length(b);

k=0;
while k<6
    k=k+1;
    for i=1:n
        if i==1
            y(i)=(b(i)-A(i,i+1:n)*x(i+1:n))/A(i,i);
        elseif i>=2 & i<=n-1
            y(i)=(b(i)-A(i,1:i-1)*x(1:i-1)-A(i,i+1:n)*x(i+1:n))/A(i,i);
        else
            y(i)=(b(i)-A(i,1:i-1)*x(1:i-1))/A(i,i);
        end
    end
   x=y;
end
x