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

Разработать алгоритм численного решения СЛАУ методом верхней релаксации - C++

Войти
Регистрация
Восстановить пароль
Другие темы раздела
C++ Структуры. Сделал http://www.cyberforum.ru/cpp-beginners/thread292733.html
Сделал сам задачку, работает правильно, но нужно было через структуры, а я не понял как. #include <iostream> #include <fstream> #include <string> using namespace std; int main() { ifstream...
C++ Найти максимально возможную площадь сарая и где он может размещаться. Фермер хочет построить на своей земле как можно больший по площади сарай. Но на его участке есть деревья и хоз. постройки, которые он не хочет никуда переносить. Для простоты представим форму сеткой... http://www.cyberforum.ru/cpp-beginners/thread292719.html
C++ В массиве А из целых чисел.....
В массиве А из чисел найти номер первого нечетного элемента.А чётные элементы этого массива заменить их индексами.
задание по матрице C++
Здравствуйте! Извините, пожалуйста, за беспокойство( Я совсем новичок( кажется, просьбы помочь с лабами не самые лучшие, но вот не знаю, что делать, поэтому выкладываю. Если не сложно, помогите,...
C++ Шаблоны функций и классов http://www.cyberforum.ru/cpp-beginners/thread292684.html
Друзья,очень нужна ваша помощь в решении следующей задачи Создать шаблонную очередь mydeque с использованием стандартного STL-контейнера vector, таким образом, чтобы добавление в конец и в...
C++ Програмирование графики Друзья срочно требуется помощь в решении следующей задачи Нужно построить траекторию движения дыух заряженных частиц,соединеных пружинкой с жесткостью k,в однородном электрическом поле с... подробнее

Показать сообщение отдельно
Jupiter
Каратель
Эксперт С++
6559 / 3980 / 227
Регистрация: 26.03.2010
Сообщений: 9,273
Записей в блоге: 1
Завершенные тесты: 2
11.05.2011, 04:57
вот вам обычный Гаусс и идите с миром...
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
88
#include <iostream>
#include <iomanip>
 
void Gauss(double **A, double *B, unsigned dim)
{
    for (unsigned k = 0; k < dim; ++k)
    {
        if (A[k][k] == 0.0)
        {
            for (unsigned i = k + 1; i < dim; ++i)
            {
                if (A[i][k] != 0.0)
                {
                    for (unsigned j = 0; j < dim; ++j)
        {
            double tmp = A[i][j];
            A[i][j] = A[k][j];
            A[k][j] = tmp;
        }
        double tmp = B[i];
        B[i] = B[k];  
        B[k] = tmp;
                   break;
                }
            }
        }
 
        double div = A[k][k];
 
        for (unsigned j = 0; j < dim; ++j)
            A[k][j] /= div;
 
        B[k] /= div;
 
        for (unsigned i = k + 1; i < dim; ++i)
        {
            double multi = A[i][k];
 
            for (unsigned j = 0; j < dim; ++j)
                A[i][j] -= multi * A[k][j];
 
            B[i] -= multi * B[k];
        }
    }
 
    for (unsigned k = dim - 1; k > 0; --k)
    {
        for (unsigned i = k - 1; i + 1 > 0; --i)
        {
            double multi = A[i][k];
 
            for (unsigned j = 0; j < dim; ++j)
                A[i][j] -= multi * A[k][j];
 
            B[i] -= multi * B[k];
        }
    }
}
 
int main()
{
    double *A[5];
    double B[5] = {22.0, 14.0, 456.0, 0.0, 5.0};
 
    for (int i = 0; i < 5; ++i)
        A[i] = new double [5];
 
    A[0][0] = 2;    A[0][1] = 0;    A[0][2] = 0;    A[0][3] = -12;  A[0][4] = 3;
    A[1][0] = 0;    A[1][1] = 2;    A[1][2] = 0;    A[1][3] = 3;    A[1][4] = 2;
    A[2][0] = 0;    A[2][1] = 0;    A[2][2] = 38;   A[2][3] = 5;    A[2][4] = 7;
    A[3][0] = -12;  A[3][1] = 13;   A[3][2] = 5;    A[3][3] = 0;    A[3][4] = 0;
    A[4][0] = 3;[4][1] = 2; A   [4][2] = 7;     A[4][3] = 0;    A[4][4] = 0;
 
    Gauss(A, B, 5);
 
    std::cout << "X = {"
              << std::setprecision(5) << B[0] << '\n'
              << std::setprecision(5) << B[1] << '\n'
              << std::setprecision(5) << B[2] << '\n'
              << std::setprecision(5) << B[3] << '\n'
              << std::setprecision(5) << B[4] 
              << "}" << std::endl;
    for (int i = 0; i < 5; ++i)
        delete [] A[i];
    std::cin.get();
 
    return 0;
}
0
 
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru