Метод хорд (не могу найти ошибку в программе)

@Blitzor DDD · Регистрация: 17.08.2015

Author24 — интернет-сервис помощи студентам

Здравствуйте,

вот написал программу, которая вычисляет корни ур-я методом хорд. Уравнение вида x^3-0.2*x^2-0.2*x-1.2=0
Интервал (1;1.5)

Код программы, алгоритм, не могу понять, почему постоянно выводит -1.#IND (вроде как на нуль нигде не делю)

Никак не могу понять, в чем косяк..

C++ (Qt)

#include <iostream>
#include <cmath>
using namespace std;
 
double f(double x)
{
    return  x*x*x-0.2*x*x-0.2*x-1.2; //Здесь вводим ту функцию, которую будем эксплуатировать
}
 
int main()
{
    double a,b,eps,X,x0;
       cout<<"Enter a= "<<endl;cin>>a;     //ввожу и считываю 1
       cout<<"Enter b= "<<endl;cin>>b;     //ввожу и считываю 1.5
       cout<<"Enter eps= "<<endl;cin>>eps; //нужная точность
       x0=a;
           do{
           //так ввели данные, что конец b - неподвижен
           X=x0-(f(x0*(b-x0))/(f(b)-f(x0)));
           x0=X;
       }while(abs(x0-X)<eps);
        cout<<endl<<"Solution is"<<x0;
}

@zss · 02.09.2015, 19:26

Сообщение от Blitzor DDD

X=x0-(f(x0*(b-x0))/(f(b)-f(x0)));

C++

1	X=x0-f(x0)*(b-x0)/(f(b)-f(x0));

@Blitzor DDD · 02.09.2015, 19:32 **[ТС]**

Да, лох я конечно.
Спасибо, но теперь вообще ничего не выводит, я поменял цикл, сделал чуть понятнее вроде бы

C++ (Qt)

 for(;;){
           //так ввели данные, что конец b - неподвижен
          // X=x0-(f(x0*(b-x0))/(f(b)-f(x0)));
           X=x0-f(x0)*(b-x0)/(f(b)-f(x0));
 
           x0=X;
           cout<<endl<<x0;
           if(abs(x0-X)<=eps) break;
          }

но теперь выводит результат неправильный. Выводит x=1.148
А результат должен быть х=1.2

Может быть, в алгоритме ошибка?

@zss · 02.09.2015, 19:44

Сообщение от Blitzor DDD

x0=X; // Тут x0 стало равно X
cout<<endl<<x0;
if(abs(x0-X)<=eps) // а тут проверяем, правда ли, что приравняли x0 к X

поменяйте 4 и 6 строки местами

@Blitzor DDD · 02.09.2015, 19:48 **[ТС]**

Благодарю Вас!
Впредь постараюсь внимательнее анализировать свой код

@zss Модератор 13706 / 10909 / 6473 Регистрация: 18.12.2011 Сообщений: 29,126
	02.09.2015, 19:44	4
	Сообщение было отмечено Blitzor DDD как решение Решение Сообщение от Blitzor DDD x0=X; // Тут x0 стало равно X cout<<endl<<x0; if(abs(x0-X)<=eps) // а тут проверяем, правда ли, что приравняли x0 к X поменяйте 4 и 6 строки местами 0

@Blitzor DDD 14 / 14 / 1 Регистрация: 17.08.2015 Сообщений: 460
	02.09.2015, 19:48 [ТС]	5
	Благодарю Вас! Впредь постараюсь внимательнее анализировать свой код 0

Опции темы

Метод хорд (не могу найти ошибку в программе)

Решение