Форум программистов, компьютерный форум CyberForum.ru

Ошибка вычисления - C++

Восстановить пароль Регистрация
 
ЗеХель
4 / 4 / 0
Регистрация: 04.11.2014
Сообщений: 125
18.01.2016, 19:35     Ошибка вычисления #1
Доброго времени суток. В результате работы везде в интегральной площади выдаёт -nan(ind). Что это и с чем связано может быть? Пример ввода координат: -1.5 и 1.5
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
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
#include <iostream>
#define _USE_MATH_DEFINES
#include <math.h>
 
using std::cout;
using std::cin;
using std::endl;
 
typedef double(*funcX)(double);
 
double integr_rignt_rect(double, double, funcX, double, double, double);
double func1(double);
double func2(double);
double func3(double);
double func4(double);
double func5(double);
double func6(double);
double func7(double);
double func8(double);
double func9(double);
double func10(double);
double func11(double);
double func12(double);
double func13(double);
 
int main()
{
    setlocale(LC_ALL, "Russian");
    cout << "Добро пожаловать\n";
    cout << "ВВедите координаты точки x, y:\n";
    double x, y;
    cin >> x >> y;
    double eps = 0.000001, step = 0.001;
 
 
    int i;
    if (x > -2 && x < 0 && y > 0 && y < 2 && x * x + (y - 1)*(y - 1) > 1)
    {
        cout << "Точка принадлежит M1\n";
        i = 1;
    }
    else
        if (x > 0 && y < 1 && x * x + (y - 1)*(y - 1) < 1)
        {
            cout << "Точка принадлежит M2\n";
            i = 2;
        }
        else
            if (y < 0 && x * x + y * y < 1 && (x - 1) * (x - 1) + y * y < 1)
            {
                cout << "Точка принадлежит M3\n";
                i = 3;
            }
            else
                if (x * x + (y + 1)*(y + 1) < 1 && x * x + y * y > 1 && (x + 1) * (x + 1) + y * y < 1)
                {
                    cout << "Точка принадлежит M4\n";
                    i = 4;
                }
                else
                    if (x > 0 && x * x + (y + 1)*(y + 1) > 1 && x < 1 && y < -1 && y > -2)
                    {
                        cout << "Точка принадлежит M5\n";
                        i = 5;
                    }
                    else
                    {
                        i = 0;
                    }
 
    cout.precision(4);
    switch (i)
    {
        double M1, M2, M3, M4, M5, sum;
    case 1:
        cout << "Геометрическая площадь = " << (M1 = 2 + 2 * (1 - M_PI/4)) << endl;
        cout << "Интегральная площадь = " << (M1 = fabs((integr_rignt_rect(-2, 0, func1, sum, step, eps) - integr_rignt_rect(-2, -1, func2, sum, step, eps) - integr_rignt_rect(-1, 0, func3, sum, step, eps)) + (integr_rignt_rect(-2, -1, func2, sum, step, eps) + integr_rignt_rect(-1, 0, func5, sum, step, eps) - integr_rignt_rect(-2, 0, func4, sum, step, eps)))) << endl;
        break;
    case 2:
        cout << "Геометрическая площадь =" << (M2 = M_PI/4 - sqrt(3)/4) << endl;
        cout << "Интегральная площадь = " << (M2 = fabs(integr_rignt_rect(0, 1, func2, sum, step, eps) - integr_rignt_rect(0, sqrt(3)/2, func6, sum, step, eps) - integr_rignt_rect(sqrt(3) / 2, 1, func5, sum, step, eps))) << endl;
        break;
    case 3:
        cout << "Геометрическая площадь =" << (M3 = sqrt(3)/4) << endl;
        cout << "Интегральная площадь = " << (M3 = fabs(integr_rignt_rect(0, 1, func4, sum, step, eps) - integr_rignt_rect(0, 0.5 , func8, sum, step, eps) - integr_rignt_rect(0.5 , 1, func7, sum, step, eps))) << endl;
        break;
    case 4:
        cout << "Геометрическая площадь =" << (M4 = sqrt(3)/2 + M_PI/3 - M_PI/4 - 1) << endl;
        cout << "Интегральная площадь =" << (M4 = fabs(integr_rignt_rect(-1, -sqrt(3)/2, func10, sum, step, eps) + integr_rignt_rect(-sqrt(3)/2, -0.5, func9, sum, step, eps) - integr_rignt_rect(-1, 0.5, func11, sum, step, eps))) << endl;
        break;
    case 5:
        cout << "Геометрическая площадь =" << (M5 = 1 - M_PI / 4) << endl;
        cout << "Интегральная площадь =" << (M5 = fabs(integr_rignt_rect(0, 1, func12, sum, step, eps) - integr_rignt_rect(0, 1, func13, sum, step, eps))) << endl;
        break;
    case 0:
        cout << "Не принадлежит областям.\n";
    }
    system("pause");
}
 
double integr_rignt_rect(double a, double b, funcX func, double s1, double step, double eps)
{
    double x;
    double h = (b - a) / step;
    x = a + h;
    s1 = 0;
    for (int i(0); i <= h; i++)
    {
        s1 += func(x);
        x = x + h;
    }
    return s1 * h;
}
 
 
double func1(double x)
{
    return 2;
}
 
double func2(double x)
{
    return 1;
}
 
double func3(double x)
{
    return (sqrt(1 - x * x) + 1);
}
 
double func4(double x)
{
    return 0;
}
 
double func5(double x)
{
    return (-sqrt(1 - x * x) + 1);
}
 
double func6(double x)
{
    return (sqrt(1 - x * x));
}
 
double func7(double x)
{
    return (-sqrt(1 - x * x) + 1);
}
 
double func8(double x)
{
    return (-sqrt(1 - x * x));
}
 
double func9(double x)
{
    return (-sqrt(1 - x * x));
}
 
double func10(double x)
{
    return (sqrt(1 - x * x) - 1);
}
 
double func11(double x)
{
    return (-sqrt(-2 * x - x * x));
}
 
double func12(double x)
{
    return (-sqrt(1 - x * x) - 1);
}
 
double func13(double x)
{
    return -2;
}
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
18.01.2016, 19:35     Ошибка вычисления
Посмотрите здесь:

Ошибка 1 error C2064: результатом вычисления фрагмента не является функция, принимающая 0 аргументов C++
Алгоритм Диффи-Хеллмана. Ошибка вычисления C++
Составить программу для вычисления выражения. В чем ошибка? C++
C++ Программа для вычисления суммы двух неотрицательных чисел. Ошибка: неправильный формат вывода.
C++ Программа для вычисления квадратного корня. Ошибка компиляции в аргументе функции
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
olya7
476 / 293 / 91
Регистрация: 18.02.2013
Сообщений: 684
18.01.2016, 19:51     Ошибка вычисления #2
Цитата Сообщение от ЗеХель Посмотреть сообщение
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
double integr_rignt_rect(double a, double b, funcX func, double s1, double ste
p, double eps)
{
* * double x;
* * [B]double h = (b - a) / step;[/B]
* * x = a + h;
* * s1 = 0;
* * for (int i(0); i <= h; i++)
* * {
* * * * s1 += func(x);
* * * *
         x = x + h;
* * }
* * return s1 * h;
}
и
Цитата Сообщение от ЗеХель Посмотреть сообщение
step = 0.001;
не очень вяжется с
Цитата Сообщение от ЗеХель Посмотреть сообщение
x = x + h;
zss
Модератор
Эксперт С++
 Аватар для zss
5950 / 5555 / 1786
Регистрация: 18.12.2011
Сообщений: 14,195
Завершенные тесты: 1
18.01.2016, 19:55     Ошибка вычисления #3
C++
1
2
3
4
5
6
7
8
9
double integr_rignt_rect(double a, double b, funcX func,double step)
{
    double s1=0;
    for (double x=a+h; x<=b;x+=step)
    {
        s1 += func(x);
    }
    return s1 * step;
}
ЗеХель
4 / 4 / 0
Регистрация: 04.11.2014
Сообщений: 125
18.01.2016, 21:40  [ТС]     Ошибка вычисления #4
zss, Ваша функция не работает

Добавлено через 52 секунды
не считает интеграл методом правого прямогольника.

Добавлено через 5 минут
olya7, а не могли бы мне помочь, пожалуйста, написать функцию для вычисления интеграла методом правого прямоугольника для моего случая.
C++
1
double integr_rignt_rect(double a, double b, funcX func,double step)
a - нижняя граница, b - верхняя, func - передаваемая функция, step - количество шагов.
zss
Модератор
Эксперт С++
 Аватар для zss
5950 / 5555 / 1786
Регистрация: 18.12.2011
Сообщений: 14,195
Завершенные тесты: 1
18.01.2016, 22:01     Ошибка вычисления #5
Цитата Сообщение от ЗеХель Посмотреть сообщение
Ваша функция не работает
В чем заключается неработоспособность этой функции?
step - это значение шага, а не количество шагов.
ЗеХель
4 / 4 / 0
Регистрация: 04.11.2014
Сообщений: 125
18.01.2016, 22:21  [ТС]     Ошибка вычисления #6
zss, она выдает 0 или какие-то астрономические числа.

Вот код, а внизу изображение задания и функции для интегрирования.
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
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
#include <iostream>
#define _USE_MATH_DEFINES
#include <math.h>
 
using std::cout;
using std::cin;
using std::endl;
 
typedef double(*funcX)(double);
 
double integr_rignt_rect(double, double, funcX, double, double, double);
double func1(double);
double func2(double);
double func3(double);
double func4(double);
double func5(double);
double func6(double);
double func7(double);
double func8(double);
double func9(double);
double func10(double);
double func11(double);
double func12(double);
double func13(double);
 
int main()
{
    setlocale(LC_ALL, "Russian");
    cout << "Добро пожаловать\n";
    cout << "ВВедите координаты точки x, y:\n";
    double x, y;
    cin >> x >> y;
    double eps = 0.000001, step = 0.001;
 
 
    int i;
    if (x > -2 && x < 0 && y > 0 && y < 2 && x * x + (y - 1)*(y - 1) > 1)
    {
        cout << "Точка принадлежит M1\n";
        i = 1;
    }
    else
        if (x > 0 && y < 1 && x * x + (y - 1)*(y - 1) < 1)
        {
            cout << "Точка принадлежит M2\n";
            i = 2;
        }
        else
            if (y < 0 && x * x + y * y < 1 && (x - 1) * (x - 1) + y * y < 1)
            {
                cout << "Точка принадлежит M3\n";
                i = 3;
            }
            else
                if (x * x + (y + 1)*(y + 1) < 1 && x * x + y * y > 1 && (x + 1) * (x + 1) + y * y < 1)
                {
                    cout << "Точка принадлежит M4\n";
                    i = 4;
                }
                else
                    if (x > 0 && x * x + (y + 1)*(y + 1) > 1 && x < 1 && y < -1 && y > -2)
                    {
                        cout << "Точка принадлежит M5\n";
                        i = 5;
                    }
                    else
                    {
                        i = 0;
                    }
 
    cout.precision(4);
    switch (i)
    {
        double M1, M2, M3, M4, M5, sum;
    case 1:
        cout << "Геометрическая площадь = " << (M1 = 2 + 2 * (1 - M_PI/4)) << endl;
        cout << "Интегральная площадь = " << (M1 = fabs((integr_rignt_rect(-2, 0, func1, sum, step, eps) - integr_rignt_rect(-2, -1, func2, sum, step, eps) - integr_rignt_rect(-1, 0, func3, sum, step, eps)) + (integr_rignt_rect(-2, -1, func2, sum, step, eps) + integr_rignt_rect(-1, 0, func5, sum, step, eps) - integr_rignt_rect(-2, 0, func4, sum, step, eps)))) << endl;
        break;
    case 2:
        cout << "Геометрическая площадь =" << (M2 = M_PI/4 - sqrt(3)/4) << endl;
        cout << "Интегральная площадь = " << (M2 = fabs(integr_rignt_rect(0, 1, func2, sum, step, eps) - integr_rignt_rect(0, sqrt(3)/2, func6, sum, step, eps) - integr_rignt_rect(sqrt(3) / 2, 1, func5, sum, step, eps))) << endl;
        break;
    case 3:
        cout << "Геометрическая площадь =" << (M3 = sqrt(3)/4) << endl;
        cout << "Интегральная площадь = " << (M3 = fabs(integr_rignt_rect(0, 1, func4, sum, step, eps) - integr_rignt_rect(0, 0.5 , func8, sum, step, eps) - integr_rignt_rect(0.5 , 1, func7, sum, step, eps))) << endl;
        break;
    case 4:
        cout << "Геометрическая площадь =" << (M4 = sqrt(3)/2 + M_PI/3 - M_PI/4 - 1) << endl;
        cout << "Интегральная площадь =" << (M4 = fabs(integr_rignt_rect(-1, -sqrt(3)/2, func10, sum, step, eps) + integr_rignt_rect(-sqrt(3)/2, -0.5, func9, sum, step, eps) - integr_rignt_rect(-1, 0.5, func11, sum, step, eps))) << endl;
        break;
    case 5:
        cout << "Геометрическая площадь =" << (M5 = 1 - M_PI / 4) << endl;
        cout << "Интегральная площадь =" << (M5 = fabs(integr_rignt_rect(0, 1, func12, sum, step, eps) - integr_rignt_rect(0, 1, func13, sum, step, eps))) << endl;
        break;
    case 0:
        cout << "Не принадлежит областям.\n";
    }
    system("pause");
}
 
double integr_rignt_rect(double a, double b, funcX func, double s1, double step, double eps)
{
    for (double i= a; i <= b; i+=step)
    {
        s1 += func(i);
    }
    return s1 * step;
}
 
 
double func1(double x)
{
    return 2;
}
 
double func2(double x)
{
    return 1;
}
 
double func3(double x)
{
    return (sqrt(1 - x * x) + 1);
}
 
double func4(double x)
{
    return 0;
}
 
double func5(double x)
{
    return (-sqrt(1 - x * x) + 1);
}
 
double func6(double x)
{
    return (sqrt(1 - x * x));
}
 
double func7(double x)
{
    return (-sqrt(1 - x * x) + 1);
}
 
double func8(double x)
{
    return (-sqrt(1 - x * x));
}
 
double func9(double x)
{
    return (-sqrt(1 - x * x));
}
 
double func10(double x)
{
    return (sqrt(1 - x * x) - 1);
}
 
double func11(double x)
{
    return (-sqrt(-2 * x - x * x));
}
 
double func12(double x)
{
    return (-sqrt(1 - x * x) - 1);
}
 
double func13(double x)
{
    return -2;
}
Миниатюры
Ошибка вычисления   Ошибка вычисления  
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
19.01.2016, 07:21     Ошибка вычисления
Еще ссылки по теме:

Ошибка в программе вычисления arctg(x) C++
Написание программы вычисления функции. Ошибка C++
C++ Ошибка вычисления суммы второго массива

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
zss
Модератор
Эксперт С++
 Аватар для zss
5950 / 5555 / 1786
Регистрация: 18.12.2011
Сообщений: 14,195
Завершенные тесты: 1
19.01.2016, 07:21     Ошибка вычисления #7
Вы не исправили функцию integr_rignt_rect
привожу ее код ЕЩЕ РАЗ:
C++
1
2
3
4
5
6
7
8
9
double integr_rignt_rect(double a, double b, funcX func,double step)
{
    double s1=0;
    for (double x=a+h; x<=b;x+=step)
    {
        s1 += func(x);
    }
    return s1 * step;
}
Yandex
Объявления
19.01.2016, 07:21     Ошибка вычисления
Ответ Создать тему
Опции темы

Текущее время: 06:38. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru