Найти ошибку в LU-разложении

 for (int i = 0; i < d; i++)
                   {
                       for (int j = 0; j <d; j++)
                       {
                           U[0,j]=A[0,j];
                       
                       if (U[0, 0] != 0)
                         {
                           L[j, 0] = A[j, 0] / U[0,0];
                         }
                       }
                   }
                   double sum = 0;
                   double summ = 0;
                   for (int i = 1; i < d; i++)
                   {
                       for (int j = i; j < d; j++)
                       {
                          for(int k1=0;k1<i-1;k1++)
                             {
                               sum+=L[i,k1]*U[k1,j];
                               summ += L[j, k1] * U[k1, i];
                             }
                          U[i, j] = A[i, j] - sum;
                          L[j, i] = (1 / U[i, i]) * (A[j, i] - summ);
                       }
                   summ = 0; sum = 0;
                   }
                
             
            }