jacobi.m

The False-Position method to solve a linear equation The Bisection method to solve linear equation

%3x1 -x2 +x3 =1
%x1 + 6x2 + 2x3 = 0
%4x1 + 3x2 + 7x3=4

clear all;
clc;
t1 = -5:0.5:5;
t2 = -5:0.5:5;
m1 = length(t1)
m2 = length(t2)
for i = 1:m1
	for j = 1:m2
		mat1(i,j) = 1-3.*t1(i) + t2(j);
mat2(i,j) = (1/2).*(-t1(i)-6.*t2(j));
mat3(i,j) = (1/7).*(4-4.*t1(i) - 3.*t2(j));
end
end
mesh(t1,t2,mat1);
hold on;
mesh(t1,t2,mat2);
mesh(t1,t2,mat3);

n = 20;
x1 = zeros(1,n);
x2 = zeros(1,n);
x3 = zeros(1,n);

x1(1) = 1;
x2(1) = 1;
x3(1) =1;
for i = 1:n
    x1(i+1) = 1/3 .*(1+x2(i)-x3(i));
	x2(i+1) = 1/6 .*(-2 .* x3(i) - x1(i));
	x3(i+1) = 1/7 .* (4 -4 .* x1(i) - 3.* x2(i));
end;
plot3(x1,x2,x3);