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

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

Войти
Регистрация
Восстановить пароль
 
Рейтинг: Рейтинг темы: голосов - 11, средняя оценка - 4.91
Ролли
Сообщений: n/a
#1

Решение системы уравнений методом Зейделя - C++

19.04.2010, 23:14. Просмотров 1376. Ответов 0
Метки нет (Все метки)

Ребята, есть такая проблема, есть код который решает сис-му уравнений методом Зейделя, но при вводе отрицательных значений выдается какой-то странный результат, скажите. что не так?
Код


#include <stdio.h>
#include <math.h>
#define N    3
#define M    300
double A[N][N]={{4,1,-5},{3,4,2},{1,5,5}};
double B[N]={11,1,-10};
double X[N],S;
float Tolerance;
int i,j;
bool Seidel(void)
{
    double Y[N];
    int i,j,l;
    double copy;
    bool Tolerance_stop_flag;
    int k=0;
    bool Iteration_stop_flag=false;
    for(i=0;i<N;i++)
    {
        Y[i]=0;X[i]=0;
    }
    do
    {
        k++;
        for(i=0;i<N;i++)
        {
            if(!A[i][i])
            {
                l=i;
                do
                {
                    l++;
                    if((!A[l][i])&&(l==N-1))
                    {
                        printf("Ошибка приведения системы для %d-го уравнения\n",i);
                        return false;
                    }
                }while(!A[l][i]);
                copy=B[i];
                B[i]=B[l];
                B[l]=copy;
                for(j=0;j<N;j++)
                {
                    copy=A[i][j];
                    A[i][j]=A[l][j];
                    A[l][j]=copy;
                }
                printf("%d-ое и %d-ое ур-ия системы переставлены\n",i,l);
            }
            S=0;
            for(j=0;j<N;j++)if(j!=i)S+=A[i][j]*X[j];
            X[i]=(B[i]-S)/A[i][i];
        }
        i=0;
        Tolerance_stop_flag=false;
        do
        {
            if(fabs(X[i]-Y[i])>Tolerance)
                Tolerance_stop_flag=true;
            else i++;
        }while(!((i==N-1)||(Tolerance_stop_flag)));
        if(!Tolerance_stop_flag)
        {
            Iteration_stop_flag=true;
            printf("Число итераций: %d\n",k);
        }else 
            for(i=0;i<N;i++)Y[i]=X[i];
    }while(!((k==M)||(Iteration_stop_flag)));
    if(!Iteration_stop_flag)
        printf("Достигнуто заданное число итераций! %d\n",M);
    return true;
}
void main (void)
{
    printf("Метод Зейделя\n");
    printf("A                B\n");
    for(i=0;i<N;i++)
    {
        for(j=0;j<N;j++)printf("%d  ",(int)A[i][j]);
        printf("            %d\n",(int)B[i]);
    }
    do{
        printf("Допустимая точность решения? ");
        scanf("%f",&Tolerance);
    }while(!((Tolerance>0)&&(Tolerance<1)));
    if(Seidel())
    {
        printf("Вектор решения X        Вектор невязок B-AX\n");
        for(i=0;i<N;i++)
        {
            S=0;
            for(j=0;j<N;j++)S+=A[i][j]*X[j];
            printf("%f            %f\n",X[i],(B[i]-S));
        }
    }
    getchar();getchar();
}
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
19.04.2010, 23:14     Решение системы уравнений методом Зейделя
Посмотрите здесь:

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

Решение системы уравнений компактным методом исключений - C++
На С++ надо решить систему уравнений компактным методом исключений. Элементы матрицы- комплексные числа, введенные в символьной форме....

Решение системы линейных уравнений методом Гаусса - C++
необходимо решить данную задачу в visual studio c++, если можно с комментариями, в консольном режиме Реализуйте функцию решения...

Решение системы дифференциальных уравнений методом Эйлера - C++
В интернете нашел похожий код, но для 1 уравнения, может кто помочь сделать для 2? Тоесть для системы #include &quot;stdafx.h&quot; ...

Решение системы линейных уравнений методом главных элементов - C++
программа на С++ решения системы линейных уравнений методом главных элементов

Решение системы линейных алгебраических уравнений методом Гаусса - C++
Решить систему n линейных алгебраических уравнений методом Гаусса. Использовать указатели.

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

Решить систему линейных уравнений вида Ax=b методом Зейделя - C++
2)Решить систему линейных уравнений вида Ax=b методом Зейделя 6.9 0,0319 0,039 0,461 0,0191 6,0 0,0333 0,405 0,0134 0,0205 5,1...

Решение СЛАУ методом Зейделя - C++
Методом Зейделя решить с точностью до 0.001 заданную систему уравнений 3.3*x1+2.1*x2+2.8*x3=0.8 4.1*x1+3.7*x2+4.8*x3=5.7 ...

Помогите найти программу решения систем алгебраических уравнений методом Зейделя - C++
Ребята, помогите найти программу решения систем алгебраических уравнений методом Зейделя. Код нужен на С++ или на С#. Переделывать с...


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

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

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