4. (A) Program for solving linear system of equations using Gauss Jordan method.

4. (A) Program for solving linear system of equations using Gauss Jordan method.

	A=[2,1,1,10;3,2,3,18;1,4,9,16]; for i=1:3
	j=i
	while (A(i,i)==0&j<=3)
	for k=1:4
	B(1,k)=A(j+1,k)
	A(j+1,k)=A(i,k)
	A(i,k)=B(1,k)
	end
	disp(A);
	j=j+1;
	end
	disp(A);
	for k=4:-1:i
	A(i,k)=A(i,k)/A(i,i)
	end
	disp(A)
	for k=1:3
	if (k~=i) then
	l=A(k,i)/A(i,i)
	for m=i:4
	A(k,m)=A(k,m)-l*A(i,m)
	end
	end
	end
	disp(A)
	end
	for i=1:3
	printf('\nx(%i)=%g\n',i,A(i,4))
	end