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

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
ЗеХель
4 / 4 / 0
Регистрация: 04.11.2014
Сообщений: 171
Завершенные тесты: 1
#1

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

18.01.2016, 19:35. Просмотров 214. Ответов 6
Метки нет (Все метки)

Доброго времени суток. В результате работы везде в интегральной площади выдаёт -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     Ошибка вычисления
Посмотрите здесь:

Большая ошибка вычисления - C++
Реализую класс комплексных чисел. и при вычислении чисел накапливается большая ошибка Например просто выводя на экран число 10exp( 9pi*I...

Ошибка в программе вычисления arctg(x) - C++
Программа открывается, но говорят, что какой-то блок не работает... Посмотрите, пожалуйста, решить данную проблему. #include...

Написание программы вычисления функции. Ошибка - C++
Написал программу, не могу понять, где описать часть условия, и вывести главное. Задание написать программу вычисления функции: ...

Алгоритм Диффи-Хеллмана. Ошибка вычисления - C++
Добрый вечер, я реализую алгоритм Диффи-Хелмана и столкнулся со следующими проблемами 1) в программе неправильно вычисляется значение k2...

Ошибка вычисления суммы второго массива - C++
Всем привет.. помогите исправить ошибку.. задание в 3 одномерных заданых масивах нада обчислить суму отрицательных элементов и вывести... ...

Составить программу для вычисления выражения. В чем ошибка? - C++
Составить программу для вычисления выражения (log2a+logbc)/logb+2(a+c): #include &lt;iostream&gt; #include &lt;cmath&gt; using namespace std; ...

Написать рекурсивную и нерекурсивную функции вычисления полинома (ошибка в цикле) - C++
Здравствуйте,Помогите найти ошибку в цикле. Задание: Написать рекурсивную и не рекурсивную функции вычисления полинома: ...

После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
olya7
477 / 294 / 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
Модератор
Эксперт С++
6282 / 5885 / 1906
Регистрация: 18.12.2011
Сообщений: 15,099
Завершенные тесты: 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
Сообщений: 171
Завершенные тесты: 1
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
Модератор
Эксперт С++
6282 / 5885 / 1906
Регистрация: 18.12.2011
Сообщений: 15,099
Завершенные тесты: 1
18.01.2016, 22:01     Ошибка вычисления #5
Цитата Сообщение от ЗеХель Посмотреть сообщение
Ваша функция не работает
В чем заключается неработоспособность этой функции?
step - это значение шага, а не количество шагов.
ЗеХель
4 / 4 / 0
Регистрация: 04.11.2014
Сообщений: 171
Завершенные тесты: 1
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     Ошибка вычисления
Еще ссылки по теме:

Программа для вычисления квадратного корня. Ошибка компиляции в аргументе функции - C++
Всем привет, программа для вычисления квадратного корня. Не компилируется из за вызова функции вычисления кв корня. Ругается на аргументы. ...

Ошибка 1 error C2064: результатом вычисления фрагмента не является функция, принимающая 0 аргументов - C++
Помогите пожалуйста разобраться не пойму почему выдает 2 ошибки? Ошибка 1 error C2064: результатом вычисления фрагмента не является...

Программа для вычисления суммы двух неотрицательных чисел. Ошибка: неправильный формат вывода. - C++
Написал супер элементарную программу для вычисления суммы двух неотрицательных чисел. #include &lt;stdio.h&gt; int main() { int a, b,...

Arduino - ошибка вычисления - Микроконтроллеры
long V; V = 1/2; Serial.println(V); в порт выводит значение &quot;0&quot;, как так?

Вычисления суммой. Ошибка - MS Access
Доброго времени суток! Прошу помощи. У меня уровень любителя в работе с Access. Проблема такая: 1. Есть состав игроков. 2....


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

Или воспользуйтесь поиском по форуму:
zss
Модератор
Эксперт С++
6282 / 5885 / 1906
Регистрация: 18.12.2011
Сообщений: 15,099
Завершенные тесты: 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     Ошибка вычисления
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru