for ( i = 0; i < n; i++ )
  {
    sum = b[i];
    for ( j = 0; j < i; j++ )
      sum -= L[i][j]*x[j];
    x[i] = sum / L[i][i];
  }