Выдаёт -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

Новые блоги и статьи Все статьи Все блоги /
Модель ЗдрввоСохранения 7: больше работников, больше ресурсов. anaschu 08.04.2026 работников и заданий может быть сколько угодно, но настроено всё так, что используется пока что только 20%	Дальние перспективы сервера - слоя сети с космологическим дизайном интефейса карты и логики. Hrethgir 07.04.2026 Дальнейшее ближайшее планирование вывело к размышлениям над дальними перспективами. И вот тут может быть даже будут нужны оценки специалистов, так как в дальних перспективах всё может очень сильно. . .	Горе от ума kumehtar 07.04.2026 Эта мне ментальная установка, что вот прямо сейчас, мол, мне для полного счастья не хватает (нужное вписать), и когда я этого достигну - тогда и полный кайф. Одна из самых сильных ловушек на пути. . . .	Использование значений реквизитов справочника в документе, с определенными условиями и правами Maks 07.04.2026 1. Контроль срока действия договора Алгоритм из решения ниже реализован на примере нетипового документа "ЗаявкаНаРаботу", разработанного в конфигурации КА2. Задача: уведомлять пользователя, если. . .
Доступность команды формы по условию Maks 07.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: сделать доступной кнопку (команда формы "ЗавершитьСписание") при. . .	Уведомление о неверно выбранном значении справочника Maks 06.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "НарядПутевка", разработанного в конфигурации КА2. Задача: уведомлять пользователя, если в документе выбран неверный склад. . .	Установка Qt Creator для C и C++: ставим среду, CMake и MinGW без фреймворка Qt 8Observer8 05.04.2026 Среду разработки Qt Creator можно установить без фреймворка Qt. Есть отдельный репозиторий для этой среды: https:/ / github. com/ qt-creator/ qt-creator, где можно скачать установщик, на вкладке Releases:. . .	AkelPad-скрипты, структуры, и немного лирики.. testuser2 05.04.2026 Такая программа, как AkelPad существует уже давно, и также давно существуют скрипты под нее. Тем не менее, прога живет, периодически что-то не спеша дополняется, улучшается. Что меня в первую очередь. . .