Форум программистов, компьютерный форум CyberForum.ru
Наши страницы

LU - разложение, для решения СЛАУ - C++

Войти
Регистрация
Восстановить пароль
Другие темы раздела
C++ Рекурсивная функция для вычисления суммы элементов арифметической прогрессии http://www.cyberforum.ru/cpp-beginners/thread745605.html
#include <stdio.h> #include <conio.h> int pr(int n, int a1, int d,int sum) { if(n==1) return a1; sum=a1+(n-1)*d+pr(n-1,a1,d); return sum; }
C++ Номер минимального простого числа в массиве одномерном C++ Упорядочить массив целых положительных чисел по возрастанию, определить номер минимального простого числа до и после упорядочивания массива. Создать функции для сортировки массива и определения... http://www.cyberforum.ru/cpp-beginners/thread745598.html
Переписать код с использование функций C++
Какие функции можно сделать из этого кода? #include "stdafx.h" #include <iostream> #include <conio.h> #include <windows.h> #include<iomanip> using namespace std;
C++ Написать программу, которая считывает три предложения и выводит их в обратном порядке
помогите пожалуйста со строками, нужно Написать программу, которая считывает три предложения и выводит их в обратном порядке
C++ Не получается знаковое умножение http://www.cyberforum.ru/cpp-beginners/thread745574.html
#ifndef INTEGERLONGS_HPP_INCLUDED #define INTEGERLONGS_HPP_INCLUDED //================================================================================================= namespace Longs { class...
C++ Объясните как работает код #include <stdio.h> #include <stdlib.h> #include <math.h> #define PRECISION 0.1e-6 double func1(size_t steps_count) { size_t i; double sum = .0; подробнее

Показать сообщение отдельно
ISunRiseI
0 / 0 / 0
Регистрация: 24.12.2012
Сообщений: 4

LU - разложение, для решения СЛАУ - C++

25.12.2012, 22:16. Просмотров 2196. Ответов 0
Метки (Все метки)

Здравствуйте. Кому не сложно можете проверить программу для решение СЛАУ, методом LU разложения. Если правильно, то можете скинуть .exe файл этой программы.
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
#include <iostream>
using namespace std;
 
int main ()
{
  int n;
  double sum = 0;
 
  cout << "Введите порядок матрицы\n n = ";
  cin >> n;
 
  double A [ n ][ n ];
  double L [ n ][ n ];
  double U [ n ][ n ];
 
//задаем матрицу A[][] ...
 
  for (int i  = 0; i < n; i++)
  {
    for (int j = 0; j < n; j++)
    {
      cout << "\na[" << i << "][" << j << "] = ";
      cin >> A [i][j];
 
      L [i][j] = 0;
      U [i][j] = 0;
 
      if (i == j)
        U [i][j] = 1;
    }
  }
 
//==============================================
 
//находим первый столбец L[][] и первую строку U[][]
 
  for (int i = 0; i < n; i++)
  {
    L [i][0] = A [i][0];
    U [0][i] = A [0][i] / L [0][0];
  }
 
//дальше вычисляем L[][], U[][] по формуле
 
  for (int i = 1; i < n; i++)
  {
      for (int j = 1; j < n; j++)
      {
          if (i >= j) //нижний треугольник
          {
              sum = 0;
              for (int k = 0; k < j; k++)
                  sum += L [i][k] * U [k][j];
 
              L [i][j] = A [i][j] - sum;
          }
          else // верхний
          {
              sum = 0;
              for (int k = 0; k < i; k++)
                  sum += L [i][k] * U [k][j];
 
              U [i][j] = (A [i][j] - sum) / L [i][i];
          }
      }
   }
 
//====================================================
   cout << "\n\n";
   for (int i = 0; i < n; i++)
   {
     for (int j = 0; j < n; j++)
       cout << "  " << L [i][j] << "   ";
     cout << "\n\n";
   }
 
   cout << "\n\n";
 
   for (int i  = 0; i < n; i++)
   {
     for (int j = 0; j < n; j++)
       cout << "  " << U [i][j] << "   ";
     cout << "\n\n";
   }
 
   return 0;
}
Добавлено через 5 часов 1 минуту
Кто-нибудь может помочь исправить программу? Надо чтобы задавался порядок матрицы,сама матрица и выводилось решение.
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
 
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru