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

Некорректная работа функции pow - C++

Восстановить пароль Регистрация
 
pro100skif
0 / 0 / 0
Регистрация: 17.03.2014
Сообщений: 41
18.03.2014, 20:36     Некорректная работа функции pow #1
Доброго времени суток. задали в универе написать прогу которая решала б кубическое уравнение.В расчетных формулах использовал функцию pow для извлечения кубического корня.Однако программа выводит неопределенность.Вот код:
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
#include<iostream>
using namespace std;
#include<math.h>
#include<conio.h>
#include<Windows.h>
#define M_PI       3.14159265358979323846
 
void vertification(double& number);
double a, b, c, p, q, d, u, v, r, f, y, y2, y3, i1, i2, x1, x2, x3;
 
 
int main()
{
    int repeat = 1;
    do{
    cout << "x^3+ax^2+bx+c=0" << endl;
    SetConsoleOutputCP(1251);
    cout << "Введіть коефіцієнт а: ";
    vertification(a);
    cout << "Введіть коефіцієнт b: ";
    vertification(b);
    cout << "Введіть коефіцієнт с: ";
    vertification(c);
    cout << "Рівняння матиме вигляд: "<<endl;
    cout << "x^3";
    if (a < 0)
        cout << a<<"x^2";
    else
        cout << "+" << a << "x^2";
    if (b < 0)
        cout << b << "x";
    else
        cout << "+" << b << "x";
    if (c < 0)
        cout << c;
    else
        cout << "+" << c;
    cout << "=0" << endl;
 
    p = b - pow(a,2) / 3;
    q = 2 * pow(a,3) / 27 - a*b / 3 + c;
    d = pow(p, 3) / 27 + pow(q, 2) / 4;
    u = pow(-q / 2 + sqrt(d), 1.0 / 3.0);
    v = -p / (3 * u);
    r = pow(-pow(p, 3) / 27, 0.5);
    f = acos(-q / 2 * r);
 
    if (d == 0)
    {
        y = 3 * p / q - a / 3;
        y2 = -3 * q / (2 * p) - a / 3;
        y3 = y2;
        x1 = y;
        x2 = y2;
        x3 = y3;
    }
    if (d < 0)
    {
        y = 2 * fabs(pow(r, 1.0 / 3.0))*cos(f / 3) - a / 3;
        y2 = 2 * fabs(pow(r, 1.0 / 3.0))*cos((f + 2 * M_PI) / 3) - a / 3;
        y3 = 2 * fabs(pow(r, 1.0 / 3.0))*cos((f + 4 * M_PI) / 3) - a / 3;
        x1 = y;
        x2 = y2;
        x3 = y3;
    }
    if (d>0)
    {
        y = u + v - a / 3;
        y2 = -(u + v) / 2 - a / 3;
        i1 = pow(3, 2)*(u - v) / 2;
    }
    cout << "Корені рівняння:" << endl;
    if (d <= 0)
    {
        cout << x1<<endl;
        cout << x2<<endl;
        cout << x3<<endl;
    }
    if (d > 0)
    {
        cout << y-a / 3 << endl;;
        cout << y2 - a / 3 << "+" << "i*" << i1 << endl;
        cout << y2 - a / 3 << "-" << "i*" << i1 << endl;
    }
    if (d <= 0)
    {
        cout << pow(x1, 3) + a*pow(x1, 2) + b*x1 + c << endl;
    }
    cout << "Для продовження натисніть ENTER, для виходу з програми натисніть ESCAPE." << endl;
    while (repeat != 13 || repeat != 27)
    {
        repeat = _getch();
        if (repeat == 13 || repeat == 27)
            break;
        cout << endl;
    }if (repeat == 27)
        break;
    } while (repeat != 27);
    return 0;
}
 
void vertification(double& number) //определение функции vertification
{
    while ((cin >> number).fail() || cin.get() != '\n')//цикл который повторяется пока не будет присвоено значение переменной number
    {
        cin.clear();
        cin.sync();
        cout << "Помилка, введено некоректні дані. Введіть дані повторно." << endl;
    }
}
Методом проб и ошибок определил что ошибка в строке
C++
1
u = pow(-q / 2 + sqrt(d), 1.0 / 3.0);
скрин с вводимыми значениями прилагается
P.S. с действительными корнями все работает корректно.И при комплексных корнях работает раз через раз.
Заранее спасибо за помощь)
Миниатюры
Некорректная работа функции pow  
Лучшие ответы (1)
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
AndrSlav
44 / 44 / 6
Регистрация: 20.12.2013
Сообщений: 241
18.03.2014, 21:36     Некорректная работа функции pow #2
Что сразу бросается в глаза - числа целые. Здесь они, очевидно, должны быть действительными. Даже если это не приводит к ошибке, то приведет позже - после цифр ставьте десятичные точки.
pro100skif
0 / 0 / 0
Регистрация: 17.03.2014
Сообщений: 41
18.03.2014, 21:38  [ТС]     Некорректная работа функции pow #3
Я проверил все значения параллельно с написание программы проверял в Маткаде.И до переменной u все біло нормально...
zss
Модератор
Эксперт С++
 Аватар для zss
5946 / 5551 / 1784
Регистрация: 18.12.2011
Сообщений: 14,175
Завершенные тесты: 1
18.03.2014, 21:42     Некорректная работа функции pow #4
Сообщение было отмечено автором темы, экспертом или модератором как ответ
Проконтролируйте, чтобы sqrt(d)-q / 2 не было отрицательным,
т.к. pow работает через вычисление логарифма:
pow(x,y)=exp(log(x)*y)
pro100skif
0 / 0 / 0
Регистрация: 17.03.2014
Сообщений: 41
18.03.2014, 21:49  [ТС]     Некорректная работа функции pow #5
Спасибо вроде помогло)
Yandex
Объявления
18.03.2014, 21:49     Некорректная работа функции pow
Ответ Создать тему
Опции темы

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