a = input('Enter the Augmented Matrix > ')
x(1)=0; y(1)=0; z(1)=0;
for i=2:100
x(i)=( a(10) - a(4)*y(i-1) - a(7)*z(i-1) )/a(1);
y(i)=( a(11) - a(2)*x(i-1) - a(8)*z(i-1) )/a(5);
z(i)=( a(12) - a(3)*x(i-1) - a(6)*y(i-1) )/a(9);
e1=abs(x(i)-x(i-1));
e2=abs(y(i)-y(i-1));
e3=abs(z(i)-z(i-1));
if e1<=0.001 && e2<=0.001 && e3<=0.001
Iteration=i
x = x(i)
y = y(i)
z = z(i)
break
end
end
x(1)=0; y(1)=0; z(1)=0;
for i=2:100
x(i)=( a(10) - a(4)*y(i-1) - a(7)*z(i-1) )/a(1);
y(i)=( a(11) - a(2)*x(i-1) - a(8)*z(i-1) )/a(5);
z(i)=( a(12) - a(3)*x(i-1) - a(6)*y(i-1) )/a(9);
e1=abs(x(i)-x(i-1));
e2=abs(y(i)-y(i-1));
e3=abs(z(i)-z(i-1));
if e1<=0.001 && e2<=0.001 && e3<=0.001
Iteration=i
x = x(i)
y = y(i)
z = z(i)
break
end
end
No comments:
Post a Comment