function x=infsolve(L,b);
n=size(L,1);
x=zeros(n,1);
for k=1:n
    y=b(k);
    for j=1:k-1
        y=y-L(k,j)*x(j);
    end
    x(k)=y/L(k,k);
end
