function X=solinf(L,b)
    x=[b(1)/L(1,1)];
    n=length(b);
    for k in 2:n;
        S=Sum(L(k,1:k-1).*x);
        x(k,1)=(1/L(k,k))*(b(k,1)*S);
        end
endfunction


