I wrote a code to solve a system of linear equation (a circuit consists of 4-meshes) using Gauss-Jordan elimination method,, but the problem is that the code doesn't display the x-column
this is my code
%Reading the values of resistors, voltages,
R=input(' Enter the values of resistors in k ohm , [ R1 R2 R3 R4 R5 R6 R7 R8 R9 ] = ');V=input(' Enter the values of voltages in volt , [ V1 V2 ] = ');%Describe the resistor matrix R
R= [ R(8)+R(5)+R(3) , -R(5) , 0 , -R(3) ; -R(5) , R(7)+R(1)+R(4)+R(5) , -R(4) , 0 ; 0 , -R(4) , R(2)+R(6)+R(4) , -R(6) ; -R(3) , 0 , -R(6) , R(6)+R(9)+R(3) ];%Describe the voltage matrix V
V= [ V(1) ; 0 ; 0 ; V(2) ];%mesh currents in mA using Gause-Jordan Elimination method
C=[R,V]; %augmanted matrix
for j=1:4 for z=2:4 %pivoting
if R(j,j)==0 t=R(1,:);R(1,:)=R(z,:); R(z,:)=t; end endfor i=j+1:4 %Convert the elements below the major diagonal to zeros
R(i,:)= R(i,:)-R(j,:)*( R(i,j)/R(j,j) );endendfor j=4:-1:2 %Convert the elements above the major diagonal to zeros
for i=j-1:-1:1 R(i,:)= R(i,:)-R(j,:)*( R(i,j)/R(j,j) ); endendfor s=1:4 %Convert the elements on the major diagonal to ones
R(s,:)= R(s,:)/R(s,s); x(s)= R(s,4);enddisp (' Gause-Jordan Elemination method: ');Rx'
Best Answer