Hi Michael,
I can't figure out the what are your variables for and working of your algorithm. But I have written a functional code for that. Before that you have to initialize the variables.
for k=1:n
if g(k,k) == 0 %if this g(k,k) becomes 0 it will display and then exit
disp('A is singular');
return;
end
b(k) = b(k)*a(k,k);
for m = k+1:n
if(k == n) % if k becomes n it will exit
return;
end
b(m) = b(m)-g(m,k)*b(k);
end
end
If you are doing forward substitution. Here is the function for forward which take lower triangular matrix L and column vector b and give result vector x:
function x=forwardsubstitution(L,b)
n = size(L,1);
for i=1:n
for j=1:i-1
b(i)=b(i)-L(i,j)*b(j);
end
b(i) = b(i)/L(i,i);
end
x=b;
end
Hope this helps!
