LU-разложение - может кто увидит ошибку (программа на С++)

@KROL64 · Регистрация: 04.10.2015

Author24 — интернет-сервис помощи студентам

Не могу понять этап нахождения самих корней с помощью LU разложения может кто увидит ошибку

Кликните здесь для просмотра всего текста

C++

void find_y(int n, double **L, double *B, double *y)
{
    double sum;
    y[0] = B[0] / L[0][0];
    for (int i = 1; i < n; i++) 
    {
        sum = 0;
        for (int j = 0; j < i; j++)
            sum += L[i][j] * y[j];
        y[i] = (B[i] - sum) / L[i][i];
    }
}
 
 
void find_x(int n, double **U, double *y, double *x)
{
    double sum;
    x[n - 1] = y[n - 1];
    for (int i = n - 2; i >= 0; i--) {
        sum = 0;
        for (int j = n-1 ; j >= i+1; j--)
            sum += x[j] * U[i][j];
        x[i] = y[i] - sum;
    }
}

B- Свободные члены L и U соответственно матрицы

@Том Ардер · 12.12.2016, 19:27

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

@ROS88 · 04.02.2017, 12:23

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

for (int j = n-1 ; j >= i+1; j--)

Попробуй поменять этот цыкл на слледующий:

C++

1	for (int j = i+1 ; j >= n; j++)

@Том Ардер $Эксперт по математике/физике$ 4217 / 3412 / 396 Регистрация: 15.06.2009 Сообщений: 5,818
	12.12.2016, 19:27	2
	KROL64, в этой ветке форума обсуждают математику, нужны уравнения и формулы. Метод известный, всё есть в учебниках. А реализация в конкретном коде - для других разделов. Примеров кода для данного метода - несть числа, на любом языке. 0

	04.02.2017, 12:23

Опции темы