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

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

Восстановить пароль Регистрация
Другие темы раздела
C++ добыть цифру из числа http://www.cyberforum.ru/cpp-beginners/thread119297.html
Задано трехзначные число.Определить равен лы квадрат этого числа сумме кубов его цифр.
C++ Симплексный метод Помгите пожалуйста написать программу которая решает задачу симплексным методом по минимизации и максимизации!! http://www.cyberforum.ru/cpp-beginners/thread119289.html
C++ Сформировать одномерный массив целых чисел.
Мужики помогите, на экзамен задали , не знаю как делать, не получается(! вот задание: Сформировать одномерный массив целых чисел. Распечатать полученный массив. Поменять местами минимальный и максимальный элементы массива. Удалить из массива все элементы превышающие его среднее значение более, чем на 100%. Распечатать полученный массив.
Массивом C++
помогите плиз.. чет проблемы с массивом у меня) Задать два массива вещественных чисел с помощью датчика случайных чисел. Для каждого элемента первого массива найти элементы с равными значениями во втором массиве (указать их индексы), если таковых не обнаружится удалить элемент из первого массива. На экран вывести как исходный, так и результирующий массивы. На печать числа выводить с двумя...
C++ Строки замена символов... http://www.cyberforum.ru/cpp-beginners/thread119247.html
Дана строка символов. Преобразовать данную строку, заменив в ней все восклицательные знаки ! точками(.), а ; каждую точку многоточием(...).
C++ Написать функцию, которая удаляет из строки s1 все символы, встречающиеся в строке s2 Народ, очень требуется помощь : :wall: № 1. Написать и протестировать функцию DELETE(s1, s2), которая удаляет из строки s1 все символы, встречающиеся в строке s2. № 2. Написать и протестировать функцию ESCAPE(str1, str2), которая при копировании текста из str1 в str2 преобразует литеры "новая строка" и "табуляция" в видимые последовательности литер \n и \t. Сделать также функцию, выполняющую... подробнее

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


#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();
}
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
 
Текущее время: 20:49. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru