Метод Гаусса

[Maray]

Добрый день!

Помогите, пожалуйста, программа неправильно решает СЛАУ.

2X1-1X2+0X3=3
0X1+5X2+2X3=7
1X1-1X2+3X3=4

Правильный ответ:
X1=2
X2=1
X3=1

Ответ, который выдает программа:
X1=2
X2=1
X3=0

#include <iostream>
using namespace std;
void sysout(double **a, double *y, int n) {
	for (int i = 0; i < n; i++) {
		for (int j = 0; j < n; j++){
			cout << a[i][j] << "*x" << j;
			if (j < n - 1) {
				cout << " + ";
			}
		}
		cout << " = " << y[i] << endl;
	}
	return;
}
int * gauss(double **a, double *y, int n) {
	double max;
	int *x, k, index;
	const double eps = 0.00001;
	x = new int[n];
	k = 0;
	while (k < n) {

		max = abs(a[k][k]);
		index = k;
		for (int i = k + 1; i < n; i++) {
			if (abs(a[i][k]) > max) {
				max = a[i][k];
				index = i;
			}
		}

		if (max < eps) {

			cout << "Решение получить невозможно из-за нулевого столбца ";
			cout << index << "матрицы A" << endl;
			return 0;
		}
		for (int j = 0; j < n; j++) {
			double temp = a[k][j];
			a[k][j] = a[index][j];
			a[index][j] = temp;
			temp = y[k];
			y[k] = y[index];
			y[index] = temp;
		}

		for (int i = k; i < n; i++) {
			double temp = a[i][k];
			if (temp < eps) continue;
			for (int j = 0; j < n; j++) {
				a[i][j] = a[i][j] / temp;
			}
			y[i] = y[i] / temp;
			if (i == k)  continue;
			for (int j = 0; j < n; j++) {
				a[i][j] = a[i][j] - a[k][j];
			}
			y[i] = y[i] - y[k];
		}
		k++;
	}

	for (k = n - 1; k >= 0; k--) {
		x[k] = y[k];
		for (int i = 0; i < k; i++) {
			y[i] = y[i] - a[i][k] * x[k];
		}
	}
	return x;
}
int main() {
	double **a, *y;
	int *x, n, g;
	system("chcp 1251");
	system("cls");
	cout << "Введите количество уравнений: ";
	cin >> n;
	a = new double*[n];
	y = new double[n];
	for (int i = 0; i < n; i++) {
		a[i] = new double[n];
		for (int j = 0; j < n; j++) {
			cout << "a[" << i << "][" << j << "]= ";
			cin >> a[i][j];
		}
	}
	for (int i = 0; i < n; i++) {
		cout << "y[" << i << "]= ";
		cin >> y[i];
	}
	sysout(a, y, n);
	x = gauss(a, y, n);
	for (int i = 0; i < n; i++){
		cout << "x[" << i << "]=" << x[i] << endl;
	}
	cin.get(); cin.get();
	return 0;
}

[min@y™]

чот уж больно наворочено...

[Maray]

Подскажите, а в моей программе никак не исправить? Если нет, то нет ли у вас кода к блок-схеме?

[min@y™]

а в моей программе никак не исправить?

Думаешь, что найдётся извращенец, который будет копаться в твоём говнокоде?

нет ли у вас кода к блок-схеме?

Конечно есть. Вот он.

[Maray]

Сообщение от min@y™

Думаешь, что найдётся извращенец, который будет копаться в твоём говнокоде?

Можно было и нормально сказать. Я ещё только изучаю С++. Не надо так называть мою программу, я её все выходные писала. Я же просто спросила

[min@y™]

Можно было и нормально сказать. Я ещё только изучаю С++. Не надо так называть мою программу, я её все выходные писала. Я же просто спросила

Добро пожаловать ва реальный мир. Тут живые люди, а не роботы.

[eoln]

Во-первых, double *x;
Во-вторых, программа неверная. То что она выдаёт 2 правильных корня чистая случайность. На любых других тестах она удачно валится. Достаточно добавить в процедуру double * gauss(double **a, double *y, int n) после k++; строку sysout(a, y, n); Там видно, что в 3-ем уравнении второй член не зануляется