Где ошибка в коде? (Реализация метода хорд)

[WhiteLedy]

(Нажмите, чтобы показать/скрыть)

Код: [Выделить]

# include <stdio.h>
# include <math.h>

double func(double x)
{
    return ((pow(2,x))+1)-pow(x,5);
}
double chord_method(a,b,entered)
{
    while(fabs(b - a) > entered)
    {
        a = b - (b - a) * func(b)/(func(b) - func(a));
        b = a - (a - b) * func(a)/(func(a) - func(b));
    }
    return b;
}
int main(){
    double entered=0.0,a=0.0,b=0.0;
    a=5.0;
    b=-5.0;
    printf("enter accuracy\n");
    scanf("%le",&entered);
    printf("%le",chord_method(a,b,entered));
    return 0;
}


это попытка реализации метода хорд для нахождения точки пересечения графиков два в степени икс плюс один и графика функции икс в пятой.
просто когда ввожу точность вычисления 0.0001 выдаёт какую-то муть

если брать пример из википедии тоже ничего вразумительного не выдаёт

(Нажмите, чтобы показать/скрыть)

Код: [Выделить]

# include <stdio.h>
# include <math.h>

double func(double x)
{
    return pow(x,3)+((-18)*x)-(83);//точность - 0.001
}
double chord_method(a,b,entered)
{
    while(fabs(b - a) > entered)
    {
        a = b - ((b - a) * func(b)/(func(b) - func(a)));
        b = a - ((a - b) * func(a)/(func(a) - func(b)));
    }
    return b;
}
int main(){
    double entered=0.0,a=0.0,b=0.0;
    a=5.0;
    b=-5.0;
    printf("enter accuracy\n");
    scanf("%le",&entered);
    printf("%le",chord_method(a,b,entered));
    return 0;
}


Самое интересное что почти то же самое но на си++ отлично работает(см. коммент ниже)

[Señor_Gaga]

(Нажмите, чтобы показать/скрыть)

Код: [Выделить]

# include <stdio.h>
# include <math.h>

double func(double x)
{
    return ((pow(2,x))+1)-pow(x,5);
}
double chord_method(a,b,entered)
{
    while(fabs(b - a) > entered)
    {
        a = b - (b - a) * func(b)/(func(b) - func(a));
        b = a - (a - b) * func(a)/(func(a) - func(b));
    }
    return b;
}
int main(){
    double entered=0.0,a=0.0,b=0.0;
    a=5.0;
    b=-5.0;
    printf("enter accuracy\n");
    scanf("%le",&entered);
    printf("%le",chord_method(a,b,entered));
    return 0;
}


это попытка реализации метода хорд для нахождения точки пересечения графиков два в степени икс плюс один и графика функции икс в пятой.
просто когда ввожу точность вычисления 0.0001 выдаёт какую-то муть

Интересно было бы для начала посмотреть вывод.

[WhiteLedy]

0.000000e+00
выводит при вводе 0.001
при работе программы №2( с примером функции из википедии)
при работе первой выдаёт другое значение -2.147484е+09
а согласно той же википедии при работе программы номер два должно выводиться число около пяти, ну, примерно, пять

[Señor_Gaga]

Прикинул как это выглядит графически:


У вас мне кажется где-то ошибка в формуле. Где, пока не вижу.
Пользователь решил продолжить мысль 27 Марта 2014, 20:56:53:Возможно ошибка в строке 6
return (pow(2,x))+1)-pow(x,5);

Так правильно?
return (pow(x,2))+1)-pow(x,5);
Пользователь решил продолжить мысль 27 Марта 2014, 21:34:23:Думаю топикастер зачет не сдаст.

В первых двух строках директивы #include записаны с ошибкой.

Вы хоть компилировать-то эту прогу пробовали?

[WhiteLedy]

изначально задача была найти пересечение функций (2^x)+1 и x^5 у меня нет проблемы с формулой
Пользователь решил продолжить мысль 27 Марта 2014, 21:53:42:

Так правильно?
return (pow(x,2))+1)-pow(x,5);

для нахождения точки пересечения графиков два в степени икс плюс

нет,как я уже говорила два в степени икс а не икс в степени два
Пользователь решил продолжить мысль 27 Марта 2014, 21:55:25:

В первых двух строках директивы #include записаны с ошибкой.

code::blocks это ест
Пользователь решил продолжить мысль 27 Марта 2014, 22:23:26:Вот что странно

(Нажмите, чтобы показать/скрыть)

Код: [Выделить]

#include <stdio.h>
#include <math.h>

double func(double x){
    return ((pow(2,x))+1)-pow(x,5);
}
double chord_method(double a,double b,double entered){
    while(fabs(b - a) > entered){
        a = b - (b - a) * func(b)/(func(b) - func(a));
        b = a - (a - b) * func(a)/(func(a) - func(b));
    }
    return b;
}
int main(){
    double entered,a,b;
    a=5;
    b=-5;
    printf("enter accuracy\n");
    scanf("%le",&entered);
    printf("%le",chord_method(a, b, entered));
    return 0;
}


моя реализация

(Нажмите, чтобы показать/скрыть)

Код: [Выделить]

#include <iostream>
#include <cmath>

double func(double x)
{
    return ((pow(2,x))+1)-pow(x,5);
}

double find(double infinum, double supremum, double epsilon)
{
    while (fabs(supremum - infinum) > epsilon)
    {
        infinum = supremum - (supremum - infinum) * func(supremum) / (func(
                supremum) - func(infinum));
        supremum = infinum - (infinum - supremum) * func(infinum) / (func(
                infinum) - func(supremum));
    }

    return supremum;
}

int main()
{
    double a = -5, b = 5;
    std::cout << find(a, b, 0.001) << std::endl;
    return 0;
}



реализация на плюсах
и последняя работает правильно а моя выдаёт вместо числа 0.6  число 6 

[WhiteLedy]

Проблема решена:

(Нажмите, чтобы показать/скрыть)

Код: [Выделить]

#include <stdio.h>
#include <math.h>

double func(double x){
    return ((pow(2,x))+1)-pow(x,5);
}
double chord_method(double a,double b,double entered){
    while(fabs(b - a) > entered){
        a = b - (b - a) * func(b)/(func(b) - func(a));
        b = a - (a - b) * func(a)/(func(a) - func(b));
    }
    return b;
}
int main(){
    double entered,a,b;
    a=5;
    b=-5;
    printf("enter accuracy\n");
    scanf("%lf",&entered);
    printf("%lf",chord_method(a, b, entered));
    return 0;
}


я заменила на %lf то, что было %le и всё заработало