метода Гауса

[Виктория Нисова]

Пожалуйста помогите мне написать программу, которая бы решала СЛАУ методом Гауса.
Нужно:
1)написать через функцию сам метод Гауса, затем ее вызвать;
2)отдельной функцие написать проверку так, что бы в ней выводилось на экран есть вообще решениЯ или их вовсе нет, если оно одно то вывести А*Х=В, где А-исходная матрица, Х-вектор неизвестных, В-то что должно получиться, сравнить свободные члены исходной матрицы с полученными.

Программа у меня написано, но не полностью... Посмотрите, пожалуйста.

PHP код:

#include "stdafx.h"
#include "conio.h"
#include <string.h>
#include "clocale"

#include <stdio.h> // Описания функций ввода-вывода
#include <math.h> // Описания математических функций
#include <stdlib.h> 

const int n=3;
void ReadData(float A[n][n+1])// функция считывания матрицы с файла
{
FILE *f;
f=fopen("input (2).txt","r");

{
printf("\t\tПрограмма находи решения системы уравнений по методу Гауса\n");
printf("Считанная матрица\n");
for (int i=0;i<n;i++)
{

    for(int j=0;j<n+1;j++)
    {
        fscanf(f,"%f ",&A[i][j]);
        printf("%f ",A[i][j]);
    }
printf("\n");
}
printf("\n\n");
fclose(f);
}
}

void check (float A[n][n+1],float R[n][n],float *Sh ,float *X)
{
    int i,j;
    float* C;
    C=new float[n];
    printf ("Проверка:\n");
    for(j=0; j<n; j++)
        {
            C[j]=0;
            for(i=0;i<n;i++) 
            C[j]=R[j][i]*X[i]+C[j];
        } 
    printf("\n(");
    for (i=0;i<n;i++)
    {   
        for(j=0; j<n; j++)
        printf ("%5.2f ",R[i][j]);
        printf(") * ( %5.6f ) = %5.2f",X[i],C[i]);
        printf("\n(");
    }
    
    for(j=0;j<n;j++)
    {printf("%4.2f ",C[j]);}
    delete C;
}


/* ------------------------------------------------- */

void Gauss(float X[n],float A[n][n+1])
{
    float B[n],temp,r; // рабочий массив свободных членов
    int i,j,k; // рабочие переменные

    // копирование в рабочую матрицу A
    for( i = 0; i < n; i++ ) 
    {
        B[i] = A[i][n]; // копирование свободных членов
        X[i] = 0;
    }

    // прямой ход
    for(i = 0;i < n;i++)
    {
        if(A[i][i] != 0)
        {
            // вычитаем строку из последующих
            for( k = (i + 1); k < n; k++ )
            {
                r = A[k][i] / A[i][i];
                for( j = i; j < n; j++ ) A[k][j] -= A[i][j] * r;
                B[k] -= B[i] * r;
            }
        }
        else 
        {
            // ищем не нулевой элемент для строки
            for(j = (i + 1);j < n;j++)
            {
                if(A[j][i] != 0) break;
            }
            // если не найден то решений нет
            
        if(j == n)
            {
                printf("единственного решения нет\n");
                exit(1);
            }
        
            
        else // если найден то меняем строки местами
            { 
                // матрица
                for(k = i;k < n;k++)
                {
                    temp = A[i][k];
                    A[i][k] = A[j][k];
                    A[j][k] = temp;
                }
                // правая часть
                temp = B[i];
                B[i] = B[j];
                B[j] = temp;
                // понижаем i на единицу чтобы цикл повторил проверки и попал в вычитание строки из последующих
                i--;
            }
        }
    }

    // обратный ход
    for( i = n; i >= 0; i-- )
    {
        float s = 0;
        for(j = (i + 1);j < n;j++) s += A[i][j] * X[j];
        X[i] = ( B[i] - s ) / A[i][i];
    }
    // вывод треугольной матрицы
    for( i = 0; i < n; i++ )
    { 
        for( j = 0; j < n; j++ ) printf("%8.2lf ", A[i][j] );
        printf( "%8.2lf", B[i] );
        printf( "\n" );
    }
}

void main( void )
{
float A[n][n+1], Sh[n], R[n][n];
setlocale(LC_ALL,"rus");
ReadData(A);

float X[n];// корни системы линейных уравнений
short i, j; // рабочие переменные
printf( "Матрица в разобранном виде\n" );

for( i = 0; i < n; i++ ) // цикл строк
{
printf( "\n" ); // новая строка на мониторе
for( j = 0; j < n+1; j++ )
 printf( "%8.2f", A[i][j] );
//printf( " %8.2lf", B[i] );

}
//printf( "\n\n" );


for( i = 0; i < n; i++ )
for( j = 0; j < n; j++ )
R[i][j]=A[i][j];



/*for( i = 0; i < n; i++ ) 
{
printf( "\n" ); //копия исходной матрицы.
for( j = 0; j < n; j++ )
 printf( "%8.2lf", R[i][j] );
}*/

for( i = 0; i < n; i++ )
    {
        Sh[i] = A[i][n]; 
    }

printf( "\nТреугольная форма:\n" );
Gauss(X,A);// решение методом Гаусса

printf( "\n" );
printf( "Решения: \n" );
printf( "\n X= " );

for( i = 0; i < n; i++ )
printf( "%10.6f", X[i] );
printf( "\n" );
    check(A,R,Sh,X);
    printf("\n\nСледовательно, решениями уравнения будут:\n");
for(i=0;i<n;i++)
printf( "%f \n",X[i]);
scanf("%d",&i);
} 


[RUSt88]

что не написано? на чем остановились?
если что-то работает не так, то указывайте, что выдает неправильно, в каком месте, если ошибки рантайм - то в какой строке произошла ошибка, что сказал отладчик
попробуйте пошагово пройтись

а правильно написано-неправильно - эт к препу

[Виктория Нисова]

Мне в ней нужно дописать (куда-то вставить) проверку на наличие корней:
если две строчки совпадают - то корней очень много,
а если две строчки различаются только свободными членами - то нет корней совсем...