Метод Гаусса в С++.

[marko)]

Привет всем.
у кого-нибудь есть программа на С++. Решение матрицы с помощью метода Гаусса с выбором главного элемента?
написанная без циклов DO-WHILE. с использованием printf, scanf, вместо cout << " " << xn[i] << " \n".

поищите, пожалуйста. очень выручите меня)

[marko)]

Вот программа, которую написал я.
она просто отнимает строки.
а вот с выбором гл. элемента. у меня что-то проблемы.))

#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#define n 3

void main()
{
float mass[n][n],sc[n],x[n];
float h;
int i,j,k;


for(i=1;i<=n;i++) // заполнение матрицы
{
for(j=1;j<=n;j++)
{
printf("vvedite ellement A[%d,%d] :",i,j);
scanf("%f",&mass[i][j]);
}
}

for(i=1;i<=n;i++)
{
printf("vvedite %d-ellement stolbca svobodnih chlenov",i);
scanf("%f",&sc[i]);
x[i]=1;
}

printf("\n");
for(i=1;i<=n;i++)
{
printf("|");
for(j=1;j<=n;j++)
{
printf("%2.1f ",mass[i][j]);
}
printf("|");
printf("%2.1f",sc[i]);
printf("\n");
}


for(i=1;i<n;i++) // исключение переменных
{
for(j=i+1;j<=n;j++)
{
mass[j][i]=(-1.0)*mass[j][i]/mass[i][i];
for(k=i+1;k<=n;k++)
{
mass[j][k]=mass[j][k]+mass[j][i]*mass[i][k];
}
sc[j]=sc[j]+mass[j][i]*sc[i];
}
}
x[n]=sc[n]/mass[n][n];

for(i=n-1;i>=1;i--) // нахождение корней
{
h=sc[i];
for(j=i+1;j<=n;j++) h=h-x[j]*mass[i][j];
x[i]=h/mass[i][i];
}

printf("korni uravneniy: \n");
for(i=1;i<=n;i++) printf("x[%d] = %1.3f\n",i,x[i]);// вывод корней

}