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

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

Восстановить пароль Регистрация
 
Рейтинг: Рейтинг темы: голосов - 11, средняя оценка - 4.91
Ролли
Сообщений: n/a
19.04.2010, 23:14     Решение системы уравнений методом Зейделя #1
Ребята, есть такая проблема, есть код который решает сис-му уравнений методом Зейделя, но при вводе отрицательных значений выдается какой-то странный результат, скажите. что не так?
Код


#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++
C++ Решение системы линейных уравнений матричным методом
Решение системы дифференциальных уравнений методом Эйлера C++
C++ Решение системы из двух нелинейных уравнений методом Ньютона
C++ Решение системы нелинейных уравнений методом Ньютона (С/С++)

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

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

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