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

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
Yuregg
Сообщений: n/a
#1

Метод Гаусса-Зейделя - C++

06.01.2014, 17:07. Просмотров 707. Ответов 0
Метки нет (Все метки)

Ребят, нужна помощь. Программа вроде бы написана правильно, корни выдает подходящие. По идее, должно выполняться около 10-20 итераций, для нахождения корней. Но, вместо этого, выполняется от 40 до 600 итераций. Если выполняется более 100 итераций, корни выдает в виде -1.IND. Что в программе написано не так?
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
#include "stdafx.h"
#include "math.h"
#include <iostream>
#include <iomanip>
#include <time.h>
using namespace std;
 
bool converge(double *xk, double *xkp, double eps, int n)
{
    for (int i = 0; i < n; i++) 
    {
        if (fabs(xk[i] - xkp[i]) >= eps) 
            return false;
    }
    return true;
}
 
int _tmain(int argc, _TCHAR* argv[])
{
    srand((unsigned)time(NULL));
    double *p, *x, **a, *b, eps;
    int n,i,j,k=0;
    cin >> n;
    cin >> eps,
    p = new double[n];
    x = new double[n];
    a = new double*[n];
    b = new double[n];
    for (i = 0; i < n; i++)
    {
        a[i] = new double[n];
        b[i] = ((double)rand()/RAND_MAX)*100;
        x[i] = 0.0;
        for (j = 0; j < n; j++)
        a[i][j] = ((double)rand()/RAND_MAX)*100;
    
    }
    cout << setw (13) << "A" << setw (16) << "B";
    cout<< endl;
    for (i = 0; i < n; i++)
    {
        cout<< endl;
        for (j = 0; j < n; j++)
            cout << a[i][j] << " ";
        cout << " " << b[i];
    }
    do
    {
        k++;
        for (int i = 0; i < n; i++)
         p[i] = x[i];
 
        for (int i = 0; i < n; i++)
        {
            double var = 0;
            for (int j = 0; j < i; j++)
                var += (a[i][j] * x[j]);
            for (int j = i + 1; j < n; j++)
                var += (a[i][j] * p[j]);
            x[i] = (b[i] - var) / a[i][i];
        }
    }
    while (!converge(x, p, eps, n));
    cout<< endl;
    for (i = 0; i < n; i++)
    {
        cout << x[i] << endl;
    }
    cout << k << endl;
    system ("pause");
 
    return 0;
}
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
06.01.2014, 17:07
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Метод Гаусса-Зейделя (C++):

СЛАУ. Метод обратной матрицы, метод Гаусса, метод Крамера, метод Зейделя - C++
Помогите ребят. Не могу построить алгоритмы для этих методов Язык C++

Метод Гаусса-Зейделя - C++
Здравствуйте. Задание дали нормальное такое) Вообще не понимаю что хотят. Объясните пожалуйста языком для домохозяек или же есть у кого...

Метод Зейделя-Гаусса С++ - C++
Для написания программы, решающей систему линейных уравнения методом итерации или Зейделя, использовалась среда разработки Visual Studio...

Метод Гаусса-Зейделя: неверный результат - C++
Здравствуйте. Я вот написал(точнее фрагмент взял с Википедии) программу решения СЛАУ методом Гаусса-Зейделя, но результат не тот, который...

Параллельный алгоритм, решенный методом Гаусса-Зейделя - C++
Всем привет.;)Нашла в инете текст программы, кот. мне нужна.Написана она на С++(параллельный алгоритм решенный методом...

Метод итераций и метод Зейделя - C++
Здравсвуйте программисты! Спасибо всем за помощь в предыдущих темах, осталась последняя лаба, которую нужно решить по предмету &quot;Численные...

Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
06.01.2014, 17:07
Привет! Вот еще темы с ответами:

Метод Зейделя в С++ - C++
всем доброго времени суток... у меня такая проблема, есть программа по решению систем линейных уравнений методом Зейделя написали мне...

Метод Зейделя - C++
Ребята, каждый из вас уже делал лабы по матрицам. Если у кого-то есть готовая лаба по методу Зейделя, прошу, напишите ее тут и желательно с...

Метод Зейделя - C++
Я написал прогу по методу Зейделя, так вот при матрице 3х3 выдает после 7 итерации правильный ответ) а при матрице 4х4 уже на второй...

метод Зейделя - C++
Ребят помогите пожалуйста, а то зачот простаивает)Нужна программка реализующая метод Зейделя(или Зейзеля препод так написал что...


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

Или воспользуйтесь поиском по форуму:
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru