Main Linear Direct 2.m
Main Linear Direct 2.m
m 1 of 3
C_mod(i,:) = C_mod(i,:)-(C_mod((i-counter),:)*(C_mod(i,(i-counter))/C_mod
((i-counter),(i-counter))));
counter = counter + 1;
end;
j = j+1;
counter = 1;
end;
C_mod
for i = 1:n
for j = 1:n
A_mod(i,j) = C_mod(i,j);
end;
end;
A_mod
Condition_Number = cond(A_mod)
S_b = 0;
x(n,1) = C_mod(n,k)/C_mod(n,n);
for i = (n-1):-1:1 %Backward Substitution
for j = n:-1:(i+1)
S_b = S_b + (C_mod(i,j)*x(j,1));
end;
x(i,1) = (C_mod(i,k)-S_b)/C_mod(i,i);
S_b = 0;
end;
x
A_mod = A;
pause();
%Gauss Jordan
fprintf ('\n\n');
disp('GAUSS JORDAN')
for i = 1:n
C_mod(i,:)= C_mod(i,:)/C_mod(i,i);
end;
C_mod
for i = 1:n
for j = 1:n
A_mod(i,j) = C_mod(i,j);
end;
end;
A_mod
Condition_Number = cond(A_mod)
k = n+1;
counter = 1;
j = n;
while (j>1)
for i = (j-1):-1:1
C_mod(i,:) = C_mod(i,:)-(C_mod((i+counter),:)*(C_mod(i,(i+counter))/C_mod
((i+counter),(i+counter))));
counter = counter + 1;
end;
j = j-1;
counter = 1;
end;
C_mod
for i = 1:n
21/1/19 12:10 PM C:\Users\admi...\main_linear_direct_2.m 3 of 3
for j = 1:n
A_mod(i,j) = C_mod(i,j);
end;
end;
A_mod
Condition_Number = cond(A_mod)
X(n,1) = C_mod(n,k)/C_mod(n,n);
for i = 1:n
X(i,1) = C_mod(i,k);
end;
X
A_mod = A;
end;