do i = 1, n sum = b(i) do j = 1, i-1 sum = sum - L(i,j)*x(j) enddo x(i) = sum / L(i,i) enddo