-1

Ответ должен быть 0.1140. А мне выводит -5.

#include <iostream>

using namespace std;

float f(int x) {
    return 7 * x + 5 ^ x - 2;
}

float Halfy(float a, float b, float e) {
    while (b - a > 2 * e) {
        float c = (a + b) / 2;
        if (f(a) * f(c) < 0) {
            b = c;
        }
        else {
            return (a + b) / 2;
        }
    }
    return (a + b) / 2;
}

int main()
{
    cout << Halfy(-10, 10, 0.0001);
}
6
  • Пробовал вычислять F(-10) = 7 * (-10) + 5 ^ (-10) - 2. Выводил 75. Хотя должен вывести −71,999.
    – Felitsius
    10 июн 2022 в 17:22
  • Хорошо что вы экспериментируете. Плохо что вы не знаете язык, систему типов и что делает операция ^. 10 июн 2022 в 17:24
  • Спасибо, но теперь другая проблема) Он округляет где то до 0.
    – Felitsius
    10 июн 2022 в 17:34
  • "... систему типов ...". 10 июн 2022 в 17:35
  • else { return (a + b) / 2; } а не a=c надо?
    – pavel
    10 июн 2022 в 17:46

1 ответ 1

2

Есть подозрение, что под 5^x вы имеете в виду 5 в степени x, но компилятор-то считает это побитовым XOR!

Вобщем, после исправления всех ошибок получается вот такое:

#include <iostream>
using namespace std;
double f(double x)
{
    return 7 * x + pow(5,x) - 2;
}

double Halfy(double a, double b, double e)
{
    while (b - a > 2 * e)
    {
        double c = (a + b) / 2;
        if (f(a) * f(c) < 0)
        {
            b = c;
        }
        else
        {
            a = c;
        }
    }
    return (a + b) / 2;
}

int main()
{
    cout << Halfy(-10, 10, 0.0001) << endl;
}

Но при этом хорошо бы еще вначале проверять, что для исходных точек получаются значения разного знака, да и каждый раз вычислять лишнее значение функции тоже не обязательно. Да и саму функцию лучше передавать как параметр — тогда код можно использовать для разных функций:

double Halfy(double a, double b, double e, double (*f)(double))
{
    double fa;
    if (a > b) swap(a,b);
    if ((fa = f(a)) * f(b) > 0) throw runtime_error("Wrong data");
    
    while (b - a > 2 * e)
    {
        double c = (a + b) / 2;
        if (fa * f(c) < 0)
        {
            b = c;
        }
        else
        {
            a = c;
        }
    }
    return (a + b) / 2;
}

double f(double x)
{
    return 7 * x + pow(5,x) - 2;
}

int main()
{
    cout << Halfy(-10, 10, 0.0001, f) << endl;
}

Ваш ответ

Нажимая на кнопку «Отправить ответ», вы соглашаетесь с нашими условиями использования и подтверждаете, что прочитали и поняли наши политику конфиденциальности и нормы поведения.

Всё ещё ищете ответ? Посмотрите другие вопросы с метками или задайте свой вопрос.