Выдаёт -nan(ind)

@Xoluz_1 · Регистрация: 29.09.2019

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

Header.h файл

C++
1
2
3
4
5
6
7
#pragma once
double f(double x);
double h(double a, double b, int n);
double hi(int i, int n, double a, double b);
double xi(int i, int n, double b, double a);
double sigma(int n, double a, double b, int i);
double integral(double a, double E, double b, int i, int n);

main.cpp файл

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
#include "Header.h"
#include <iostream>
using namespace std;
int main()
{
    int i = 0;
    double n;
    double a, b, E;
    n = 2;
    cout << "input interval [a,b] =" << endl;
    cin >> a;
    cin >> b;
    cout << "input E=" << endl;
    cin >> E;
    cout << "integral is" << integral(a, E, b, n, i) << endl;
}

func.cpp файл

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
#include "Header.h"
#include <stdio.h>
#include <cmath>
double f(double x)
{
    return x;
}
double h(double a, double b, int n)
{
    return (b - a) / n;
}
double hi(int i, int n, double b, double a)
{
    return i * (b - a) / n;
}
double xi(double a, double b, int i, int n)
{
    return a + hi(i, a, b, n) - h(a, b, n) / 2;
}
double sigma(int n, double a, double b, int i)
{
    double sum = 0;
    for (i = 1; i < n; i++)
    {
        sum += f(xi(a, b, i, n));
    }
    return h(a, b, n) * sum;
}
double integral(double a, double b, double E, int n, int i)
{
    double s1, s2;
    s1 = sigma(n, a, b, i);
    s2 = s1;
    do
    {
        s1 = s2;
        n *= 2;
        s2 = sigma(a, n, b, i);
    } while (fabs(s1 - s2) > E);
    return s2;
}

@AndryS1 · 13.01.2020, 19:37

Пользуйтесь отладчиком, добавьте везде, где есть деление, проверку на 0

А зачем вот это?

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

double f(double x)
{
return x;
}

@FFPowerMan · 13.01.2020, 19:52

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

C++
1
cout << "integral is" << integral(a, E, b, n, i) << endl;

- вот одна из ошибок. Должно быть integral(a, b, E, n, i).

Добавлено через 1 минуту
Интеграл вообще так не считается: много функций. Все пишется в одном цикле.

@AndryS1 · 13.01.2020, 20:05

C++
1
2
3
4
5
6
7
8
double hi(int i, int n, double b, double a)
{
    return i * (b - a) / n;
}
double xi(double a, double b, int i, int n)
{
    return a + hi(i, a, b, n) - h(a, b, n) / 2;
}

тут тоже напутаны значения

Добавлено через 31 секунду

C++
1
2
3
4
5
6
7
8
double hi(int i, int n, double b, double a)
{
    return i * (b - a) / n;
}
double xi(double a, double b, int i, int n)
{
    return a + hi(i, n, b, a) - h(a, b, n) / 2;
}

Добавлено через 1 минуту

C++
1
        s2 = sigma(a, n, b, i);

тут тоже напутано

@FFPowerMan · 13.01.2020, 20:06

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

C++
1
double hi(int i, int n, double a, double b);

- вот 2-ая ошибка. Должно быть так: double hi(int i, double a, double b, int n);

@AndryS1 · 13.01.2020, 20:06

C++
1
    double n;

хотя везде он int

@FFPowerMan · 13.01.2020, 20:10

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
#include <iostream>
#include <stdio.h>
#include <cmath>
using namespace std;
 
double f(double x);
double h(double a, double b, int n);
double hi(int i, int n, double a, double b);
double xi(int i, int n, double b, double a);
double sigma(int n, double a, double b, int i);
double integral(double a, double E, double b, int i, int n);
 
int main()
{
    int i = 0;
    double n;
    double a, b, E;
    n = 2;
    cout << "Input interval [a,b] =" << endl;
    cin >> a;
    cin >> b;
    cout << "Input E=" << endl;
    cin >> E;
    cout << "Integral = " << integral(a, b, E, n, i) << endl;
}
 
double f(double x)
{
    return x;
}
double h(double a, double b, int n)
{
    return (b - a) / n;
}
double hi(int i, double a, double b, int n)
{
    return i * (b - a) / n;
}
double xi(double a, double b, int i, int n)
{
    return a + hi(i, a, b, n) - h(a, b, n) / 2;
}
double sigma(int n, double a, double b, int i)
{
    double sum = 0;
    for (i = 1; i < n; i++)
    {
        sum += f(xi(a, b, i, n));
    }
    return h(a, b, n) * sum;
}
double integral(double a, double b, double E, int n, int i)
{
    double s1, s2;
    s1 = sigma(n, a, b, i);
    s2 = s1;
    do
    {
        s1 = s2;
        n *= 2;
        s2 = sigma(n, a, b, i);
    }while(fabs(s1 - s2) > E);
    return s2;
}

@AndryS1 · 13.01.2020, 20:26

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

double n;

Мне кажется целесообразнее будет использовать int

Новые блоги и статьи Все статьи Все блоги /
Unity 4D GameUnited 13.06.2025 Четырехмерное пространство. . . Звучит как что-то из научной фантастики, правда? Однако для меня, как разработчика со стажем в игровой индустрии, четвертое измерение давно перестало быть абстракцией из. . .	SSE (Server-Sent Events) в ASP.NET Core и .NET 10 UnmanagedCoder 13.06.2025 Кажется, Microsoft снова подкинула нам интересную фичу в новой версии фреймворка. Работая с превью . NET 10, я наткнулся на нативную поддержку Server-Sent Events (SSE) в ASP. NET Core Minimal APIs. Эта. . .	С днём независимости России! Hrethgir 13.06.2025 Решил побеседовать, с утра праздничного дня, с LM о завоеваниях. То что она написала о народе, представителем которого я являюсь сам сначала возмутило меня, но дальше только смешило. Это чисто. . .	Лето вокруг. kumehtar 13.06.2025 Лето вокруг. Наполненное бурями и ураганами событий. На фоне магии Жизни, священной и вечной, неумелой рукой человека рисуется панорама душевного непокоя. Странные серые краски проникают и. . .	Популярные LM модели ориентированы на увеличение затрат ресурсов пользователями сгенерированного кода (грязь -заслуги чистоплюев). Hrethgir 12.06.2025 Вообще обратил внимание, что они генерируют код (впрочем так-же ориентированы разработчики чипов даже), чтобы пользователь их использующий уходил в тот или иной убыток. Это достаточно опытные модели,. . .
Топ10 библиотек C для квантовых вычислений bytestream 12.06.2025 Квантовые вычисления - это та область, где теория встречается с практикой на границе наших знаний о физике. Пока большая часть шума вокруг квантовых компьютеров крутится вокруг языков высокого уровня. . .	Dispose и Finalize в C# stackOverflow 12.06.2025 Работая с C# больше десяти лет, я снова и снова наблюдаю одну и ту же историю: разработчики наивно полагаются на сборщик мусора, как на волшебную палочку, которая решит все проблемы с памятью. Да,. . .	Повышаем производительность игры на Unity 6 с GPU Resident Drawer GameUnited 11.06.2025 Недавно копался в новых фичах Unity 6 и наткнулся на GPU Resident Drawer - штуку, которая заставила меня присвистнуть от удивления. По сути, это внутренний механизм рендеринга, который автоматически. . .	Множества в Python py-thonny 11.06.2025 В Python существует множество структур данных, но иногда я сталкиваюсь с задачами, где ни списки, ни словари не дают оптимального решения. Часто это происходит, когда мне нужно быстро проверять. . .	Работа с ccache/sccache в рамках C++ Loafer 11.06.2025 Утилиты ccache и sccache занимаются тем, что кешируют промежуточные результаты компиляции, таким образом ускоряя последующие компиляции проекта. Это означает, что если проект будет компилироваться. . .