VPF::Решение уравнения методом простых итераций

Eisfee

Дата 25.3.2011, 15:33 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 2
Регистрация: 25.3.2011

Репутация: нет
Всего: нет

Здравствуйте. Прошу помочь найти и устранить ошибку в программе, которая должна решать уравнение вида f(x)=0 (у=-х*х+2*sin(x)), методом простой итерации. Интервал на котором находится корень [1;2], погрешность вычисления 0.001, корень найденный с помощью графика Excel равен приблизительно 1.40, там же был определен интервал.

Код


#include <conio.h>
#include <math.h>
#include <iostream.h>
#define e 0.001
double f(double x) 
{
 
        return   -x*x+2*sin(x);
}
 
double fi(double x, double L) 
{
         return  x+L*(-x*x+2*sin(x));
}
  
int main() 
{
        int n=0;
        double x,y,c,b,L=2;
        cout<<"x="; cin>>x;
        do 
       {
                y=fi(x,L);
                b=fabs(x-y);
                x=y;
                n+=1;
        }
        while (b>=e);
                cout<<"c="<<x<<"\n";
                cout<<"n="<<n<<"\n";
                getch();
        return 0;
}

Проблема скорее всего заключается в том, что я плохо поняла этот метод и реализовала его с ошибками.

bsa

Дата 25.3.2011, 15:45 (ссылка)

(нет голосов)

Загрузка ...

Эксперт

Профиль
Группа: Модератор
Сообщений: 9185
Регистрация: 6.4.2006
Где: Москва, Россия

Репутация: 63
Всего: 196

в строке 24 ты вычисляешь модуль разницы x и y. А на самом деле, тебе нужен модуль разницы текущего и предыдущего значения y.

--------------------

Правильно заданный вопрос - половина ответа
Отзывы о работе модератора.

Eisfee

Дата 25.3.2011, 15:52 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 2
Регистрация: 25.3.2011

Репутация: нет
Всего: нет

bsa, а как тогда записать эту разницу? Подскажи, пожалуйста.

0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| C/C++: Общие вопросы \| Следующая тема »