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

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

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

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

18.01.2016, 19:35. Просмотров 234. Ответов 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;
}
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
18.01.2016, 19:35
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Ошибка вычисления (C++):

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

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

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

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

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

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

6
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;
0
zss
Модератор
Эксперт С++
6403 / 5969 / 1945
Регистрация: 18.12.2011
Сообщений: 15,387
Завершенные тесты: 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;
}
0
ЗеХель
4 / 4 / 0
Регистрация: 04.11.2014
Сообщений: 182
Завершенные тесты: 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 - количество шагов.
0
zss
Модератор
Эксперт С++
6403 / 5969 / 1945
Регистрация: 18.12.2011
Сообщений: 15,387
Завершенные тесты: 1
18.01.2016, 22:01 #5
Цитата Сообщение от ЗеХель Посмотреть сообщение
Ваша функция не работает
В чем заключается неработоспособность этой функции?
step - это значение шага, а не количество шагов.
0
ЗеХель
4 / 4 / 0
Регистрация: 04.11.2014
Сообщений: 182
Завершенные тесты: 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;
}
0
Миниатюры
Ошибка вычисления   Ошибка вычисления  
zss
Модератор
Эксперт С++
6403 / 5969 / 1945
Регистрация: 18.12.2011
Сообщений: 15,387
Завершенные тесты: 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;
}
0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
19.01.2016, 07:21
Привет! Вот еще темы с ответами:

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

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

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

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


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

Или воспользуйтесь поиском по форуму:
7
Yandex
Объявления
19.01.2016, 07:21
Ответ Создать тему
Опции темы

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