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

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Web (WebAssembly): Обработчик клика мыши в браузере ПК и касания экрана в браузере на мобильном устройстве 8Observer8 02.02.2026 Содержание блога Для начала пошагово создадим рабочий пример для подготовки к экспериментам в браузере ПК и в браузере мобильного устройства. Потом напишем обработчик клика мыши и обработчик. . .	Философия технологии iceja 01.02.2026 На мой взгляд у человека в технических проектах остается роль генерального директора. Все остальное нейронки делают уже лучше человека. Они не могут нести предпринимательские риски, не могут. . .	SDL3 для Web (WebAssembly): Вывод текста со шрифтом TTF с помощью SDL3_ttf 8Observer8 01.02.2026 Содержание блога В этой пошаговой инструкции создадим с нуля веб-приложение, которое выводит текст в окне браузера. Запустим на Android на локальном сервере. Загрузим Release на бесплатный. . .	SDL3 для Web (WebAssembly): Сборка C/C++ проекта из консоли 8Observer8 30.01.2026 Содержание блога Если вы откроете примеры для начинающих на официальном репозитории SDL3 в папке: examples, то вы увидите, что все примеры используют следующие четыре обязательные функции, а. . .
SDL3 для Web (WebAssembly): Установка Emscripten SDK (emsdk) и CMake для сборки C и C++ приложений в Wasm 8Observer8 30.01.2026 Содержание блога Для того чтобы скачать Emscripten SDK (emsdk) необходимо сначало скачать и уставить Git: Install for Windows. Следуйте стандартной процедуре установки Git через установщик. . . .	SDL3 для Android: Подключение Box2D v3, физика и отрисовка коллайдеров 8Observer8 29.01.2026 Содержание блога Box2D - это библиотека для 2D физики для анимаций и игр. С её помощью можно определять были ли коллизии между конкретными объектами. Версия v3 была полностью переписана на Си, в. . .	Инструменты COM: Сохранение данный из VARIANT в файл и загрузка из файла в VARIANT bedvit 28.01.2026 Сохранение базовых типов COM и массивов (одномерных или двухмерных) любой вложенности (деревья) в файл, с возможностью выбора алгоритмов сжатия и шифрования. Часть библиотеки BedvitCOM Использованы. . .	SDL3 для Android: Загрузка PNG с альфа-каналом с помощью SDL_LoadPNG (без SDL3_image) 8Observer8 28.01.2026 Содержание блога SDL3 имеет собственные средства для загрузки и отображения PNG-файлов с альфа-каналом и базовой работы с ними. В этой инструкции используется функция SDL_LoadPNG(), которая. . .