Как правильно вывести результат

[Ленащщ]

Решение дифференциального уравнения методом Эйлера-Коши
Не могу понять, как правильно вывести результат

double func(double x, double y) 
{
    if (((pow(y, 3) + x + 1))!= 0)
        return (3 * pow(y, 2)) / (pow(y, 3) + x + 1);
}
 
int main()
{
    double x0, y0, xn;
    setlocale(0, "");
    cout << "х0=";
    cin >> x0;
    cout << endl;
    cout << "y0=";
    cin >> y0;
    cout << endl;
    cout << "xn=";
    cin >> xn;
    cout << endl;
 
    const int n = 10;   
    double h = (xn - x0) / n;
    double X[(int)n];
    double Y1[(int)n];
    double Y[(int)n];
 
    X[0] = x0; Y[0] = 1;
    for (int i = 1; i <= n; i++) 
    {
        X[i] = x0 + i * h;
        Y1[i] = Y[i - 1] + h * func(X[i - 1], Y[i - 1]);
        Y[i] = Y[i - 1] + h * (func(X[i - 1], Y[i - 1]) + func(X[i], Y1[i])) / 2.0;
    }
    return 0;
}

[digitalis]

Сама функция какая-то странная:

 if (((pow(y, 3) + x + 1))!= 0)
        return (3 * pow(y, 2)) / (pow(y, 3) + x + 1);}

а если = 0 - что на выходе ? где else ? Если = 0, то - хавайся у бульбу ?
Не говоря уже о том, что проверять double на = или != - всё равно, что переходить улицу в час пик с завязанными глазами.
Но учебные задания бывают и ещё более идиотическими.

[BLACK_RAIN]

Сообщение от digitalis

Не говоря уже о том, что проверять double на = или != - всё равно, что переходить улицу в час пик с завязанными глазами

а как тогда?

[digitalis]

Сообщение от BLACK_RAIN

а как тогда?

Если вопрос об переходе улицы - то на зелёный свет с открытыми глазами.
Если об этой странной функции, то она д.быть определена для любого значения аргумента.
А если о проверке на равенство real, то это азбучная истина: 2*2 м.быть не равно 2+2 для real, поэтому проверяют на
Abs( разность двух значений ) < eps.

[BLACK_RAIN]

Сообщение от digitalis

Abs( разность двух значений ) < eps.

подробнее, пожалуйста

[digitalis]

Ну какое тут "подробно" . Из-за особенностей машинной арифметики значение, полученное одним путём, может быть не равным тому же значению, полученному другим путём, хотя вычисляя на бумажке - знаем, что одно и то же. Не приходилось видеть, к примеру, в результатах 2.99999, когда мы точно знаем, что д.быть 3.0 ? И если тупо сравнивать на =, то равные значения могут оказаться хоть на 1 бит отличающимися, и тогда прокрутим (или недокрутим) один лишний цикл, например. Поэтому для real типично: if Abs(var1 - var2) < eps
Я спросил у Яндекса: "сравнение float на равенство"... Первая же ссылка на Хабр.

[BLACK_RAIN]

Сообщение от digitalis

if Abs(var1 - var2) < eps

error: 'eps' was not declared in this scope
к сожалению

[BDA]

Нужно самому объявить eps со значением желаемой точности сравнения.

[BLACK_RAIN]

Сообщение от BDA

Нужно самому объявить eps со значением желаемой точности сравнения.

а как понять, какой точности достаточно?

[digitalis]

Это каждый решает сам - сколько точности ему достаточно для данной задачи.
"
Каждый выбирает для себя -
Женщину, религию, дорогу...
" © Ю. Д. Левитанский
ИЧСХ : ТС давно слиняла, а дяди тут рвут рубахи в поисках докопаться до истины