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

Зейдель - C++

Восстановить пароль Регистрация
 
qswd
0 / 0 / 0
Регистрация: 16.02.2013
Сообщений: 171
08.04.2013, 19:46     Зейдель #1
прошу помощи...вроди все по формуле... но не работает...
http://ru.wikipedia.org/wiki/%D0%9C%...B5%D0%BB%D1%8F
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
#include <iostream>
#include <math.h>
using namespace std;
 
const int size = 3;
 
int main()
{
    double matrix[size][size+1];
    double answer[size];
    matrix[0][0] = 1;   matrix[0][1] = 1;   matrix[0][2] = -2;  matrix[0][3] = 6;
    matrix[1][0] = 2;   matrix[1][1] = 3;   matrix[1][2] = -7;  matrix[1][3] = 16;
    matrix[2][0] = 5;   matrix[2][1] = 2;   matrix[2][2] = 1;   matrix[2][3] = 16;
 
    double oldanswer[size] = {0};
 
    while ( true )
    {
        for ( int i = 0; i < size; ++i )
        {
            answer[i] = matrix[i][size];
            for ( int j = 0; j < size; ++j )
            {
                if ( j < i )
                    answer[i] -= matrix[i][j] * answer[j];
                if ( j > i )
                    answer[i] -= matrix[i][j] * oldanswer[j];
            }
            answer[i] /= matrix[i][i];
        }
 
        int count = 0;
        for ( int k = 0; k < size; ++k)
            if( fabs( answer[k] - oldanswer[k] ) < 0.01 )
                ++count;
 
        if ( count == size )
            break;
 
        for ( int k = 0; k < size; ++k )
           oldanswer[k] = answer[k];
    }
 
    for ( int i = 0 ; i < size; ++i )
        cout << answer[i] << "\t";
    cout << endl;
 
    return 0;
}
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
08.04.2013, 19:46     Зейдель
Посмотрите здесь:

Scilab Методы решения СЛАУ (Гаусс, Зейдель)

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Ответ Создать тему
Опции темы

Текущее время: 13:17. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru