do_gauss.m

辛普森方法

%do_gauss
A = [3 2;1 2];
b = [1 -1]';
% x = gauss(A,b)
% x = jacobi(A,b)
% x = gauseid(A,b)
x = SOR(A,b,1.0)
x1 = A\b %for comparison with the result of backslash operation