Jacobi method to solve equation using MATLAB(mfile)
n=input('Enter number of equations,
n: ');
A =
zeros(n,n+1);
x1 =
zeros(n);
x2 =
zeros(n);
tol = input('Enter the tolerance, tol: ');
m = input('Enter maximum number of
iterations, m: ');
A=[4 2 3 8; 3
-5 2 -14; -2 3 8 27];
x1=[0 0 0];
k = 1;
while k <= m
err = 0;
for i = 1 : n
s = 0;
for j = 1 : n
s =
s-A(i,j)*x1(j);
end
s =
(s+A(i,n+1))/A(i,i);
if abs(s) > err
err =
abs(s);
end
x2(i) =
x1(i)+s;
end
if err <= tol
break;
else
k = k+1;
for i = 1 : n
x1(i) =
x2(i);
end
end
end
fprintf('Solution vector after %d
iterations is :\n', k-1);
for i = 1 : n
fprintf(' %11.8f \n', x2(i));
end
>> jacobimethod
Enter number of equations, n: 3
Enter the tolerance, tol: 0.001
Enter maximum number of iterations, m: 100
Solution vector after 41 iterations is :
-0.99919498
2.99974076
2.00028882
>>
May I have a question, for this code, how you can plot the point x in the figure plot?
ReplyDelete
ReplyDeleteThanks for sharing such a Excellent Blog! We are linking too this particularly great post on our site. Keep up the great writing.
matlab assignment help
شركة تنظيف منازل بالدمام
ReplyDeleteشركة تنظيف منازل بالقصيم
شركة تنظيف منازل بالاحساء
شركة تنظيف منازل بالقطيف
شركة تنظيف منازل بالجبيل