треугольный вид матрицы

[Pein95]

for (int k=0;k<n;k++)
       {
          for (int i=k+1;i<n;i++)
          {
              if (a[k][k] != 0)
             mu=a[i][k]/a[k][k];
             for (int j=0;j<n;j++)
                a[i][j]-=a[k][j]*mu;
           
          }
       }

помогите исправить етот код.
Пробовал найти чтото в инете но там как то всё очень сложно. Мне не нужно решить систему уравнений. Нужно только свести матрицу к треугольному виду.

[Pein95]

есть тут кто?)

[Stilet]

Я есть. Но я не понял задачи.

[Perchik71]

квадрат не может быть треугольником. поймите вашу задачу более правильно.

[EUGY]

"Треугольник будет выпит. Будь он параллелепипед, будь он круг, едрена вошь!"

На первом этапе осуществляется так называемый прямой ход, когда путём элементарных преобразований над строками систему приводят к ступенчатой или треугольной форме, либо устанавливают, что система несовместна. А именно, среди элементов первого столбца матрицы выбирают ненулевой, перемещают его на крайнее верхнее положение перестановкой строк и вычитают получившуюся после перестановки первую строку из остальных строк, домножив её на величину, равную отношению первого элемента каждой из этих строк к первому элементу первой строки, обнуляя тем самым столбец под ним. После того, как указанные преобразования были совершены, первую строку и первый столбец мысленно вычёркивают и продолжают пока не останется матрица нулевого размера. Если на какой-то из итераций среди элементов первого столбца не нашёлся ненулевой, то переходят к следующему столбцу и проделывают аналогичную операцию.

http://ru.wikipedia.org/wiki/%D0%9C%...81%D1%81%D0%B0

[Perchik71]

Сообщение от EUGY

http://ru.wikipedia.org/wiki/%D0%9C%...81%D1%81%D0%B0

ты перепил...... ты понимаешь, даже лучше трезвых.

[Pein95]

спс) всем кто откликнулся) проблему уже решил)

[Stilet]

Интересно как...

[Pein95]

void triangle_view(double a[][N],int n)
{
    double v;
    int row = 0;
    int k=0;
    while (k<n-1)
   {
       int i,j;
	   int gr = row;
        for (i=row+1; i<n;i++)
        {
            if (fabs(a[row][k]) < fabs(a[i][k]))
                gr = i;
        }

        if ( (row == gr) && (a[row][k] == 0)) {
			k++;
			continue;
		}


        if (row != gr)
        {
            for (j=0; j<n; j++)
            {
				
                v = a[row][j]; 
                a[row][j] = a[gr][j];
                a[gr][j] = v;
            }

        }

	
        for (i=row+1; i<n; i++)
        {
			
            v = 1.0*a[i][k]/a[row][k];
            a[i][k] = 0;
            if (v!=0) 
                for(j=k+1; j<n; j++)
                    a[i][j] = a[i][j] - v*a[k][j];
			
        }
    row++;
    k++;
    }
	
}

так)