При выводе числа выводится #INF

@Frip · Регистрация: 04.11.2013

Студворк — интернет-сервис помощи студентам

В общем, задали написать программу вычисления кубического уравнения. В сети нашёл алгоритм решения (ссыль, который можно спокойно использовать СИ++ без гиперболических функций (Привет формуле Виета).

Все написал, ошибок не выдаёт при компиляции, но вот при вводе любых значений выводится #INF. Поискал, написано , что в большинстве связано с делением на ноль, но таких ошибок, мной допущено не было в написании, вроде.

Посмотрите, может где увидите ошибку (И если нашли и справили, поставьте комментарий, хочу посмотреть, где ошибся)

Кликните здесь для просмотра всего текста

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
#include <iostream>
#include <math.h>
 
using namespace std;
 
int main()
{
 double A,B;
 double a,b,c,d;
 int r;
 cout<<"ax^3+bx^2+cx+d=0"<<endl;
     cout<<"Vvedite a: ";
           cin>>a;
     cout<<"Vvedite b: ";
           cin>>b;                         
     cout<<"Vvedite c: ";
           cin>>c;                                      
     cout<<"Vvedite d: ";
           cin>>d;
     cout<<"\n";      
     //Избавляемся от коэф. перед x^3
         a=b/a;
         b=c/a;
         c=d/a;
         
     //Вычисляем Q и R для метода Виета-Кордано
         double Q = (a*a-3.0*b)/9.0;
         double R = (2.0*a*a*a-9.0*a*b+27.0*c)/54;
         
     /*Сравните полученные значения R и Q. 
Если верно выражение Q³ >R² , следовательно, 
в исходном уравнении присутствуют 3 действительных корня. 
Вычислите их по формулам Виета.*/
 
         if ((Q*Q*Q)>(R*R))
            {
              cout<<"Q^3 > R^2, znachit imeem 3 deistvitelnih kornya"<<endl;
              double t = (acos(R/sqrt(Q*Q*Q)))/3.0;
                  double x1 = -2.0 * sqrt(Q) * cos(t) - a/3.0;
                  double x2 = -2.0 * sqrt(Q) * cos(t+(2*M_PI/3.0))-a/3.0;
                  double x3 = -2.0 * sqrt(Q) * cos(t-(2*M_PI/3.0))-a/3.0;
              
              cout<<"x1= "<<x1<<endl;
              cout<<"x2= "<<x2<<endl;
              cout<<"x3= "<<x3<<endl;             
                           }
     /* При значениях Q³ <= R² , в решении находится один действительный корень х1 
и два комплексно-сопряженных корня. Для их определения нужно 
найти промежуточные значения А и В. Вычислите их по формулам Кардано.*/
           if ((Q*Q*Q)<=(R*R))
              {
               if (R>0) r=1;
               if (R<0) r=-1;
               if (R=0) r=0;
             A =-r*pow((fabs(R)+sqrt(R*R-Q)), 1/3);
                    if (A=0){B=0;} else B=Q/A;
              
              if (A!=B)
              {      
                double x1 = (B+A)-a/3.0;
                double x23_1 = -(A+B)/2.0+a/3.0;
                double x23_2 = sqrt(3.0)*(A-B)/2.0;
                       cout<<"x1= "<<x1<<endl;
                       cout<<"x2= "<<x23_1<<"+i*"<<x23_2<<endl;
                       cout<<"x3= "<<x23_1<<"-i*"<<x23_2<<endl;
           
                } else
                      {
                       double x1 = (B+A)-a/3.0;
                       double x2 = -A-a/3.0;
                              cout<<"x1= "<<x1<<endl;
                              cout<<"x2= "<<x2<<endl;           
                      }
                                     
                              }
      system ("pause");     
}

@olgashat · 09.12.2014, 20:24

Сообщение от olgashat

phi = log(q / (r*r*r) - sqrt(pow(q / (r*r*r),2)-1));

а должно было быть

C++
1
phi = log(q / (r*r*r) + sqrt(pow(q / (r*r*r),2)-1));

и то же самое в 5 ветке

@Frip · 09.12.2014, 20:59 **[ТС]**

ааа, вы про это? Я это изменил, результат я огласил уже вроде, ответы всё-равно не те

@olgashat · 09.12.2014, 22:45

Frip, а что за книга со скриншотами? дело в том, что я получила промежуточные ответы, но у меня в них не сходятся мнимые части... похожей матчасти не могу найти, подскажите автора

все, не надо, все сошлось)

Добавлено через 15 минут

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
#include <iostream>
#include <math.h>
int sign(double x)
{
    if (x > 0) return 1;
    else
        if (x < 0) return -1;
        else return 0;
 
}
using namespace std;
 
int main()
{
    const double PI = 3.14159265358979323846;
    double r, phi;
    double a, b, c, d;
    
    /*cout << "ax^3+bx^2+cx+d=0" << endl;
    cout << "Vvedite a: ";
    cin >> a;
    cout << "Vvedite b: ";
    cin >> b;
    cout << "Vvedite c: ";
    cin >> c;
    cout << "Vvedite d: ";
    cin >> d;
    cout << "\n";*/
    a = -1;
    b = -2;
    c = 3;
    d = 1;
 
    double p = (3.0*a*c - b*b) / (9.0*a*a);
    double q = (b*b*b) / (27.0*a*a*a) - (b*c) / (6.0*a*a) + d / (2*a); //
    double q2 = q*q;
    double p3 = p*p*p;
    
    cout << "p=" << p << endl;
    cout << "q=" << q << endl;
    cout << "Q=q*q + p*p*p=" << q*q + p*p*p << endl;
    /* 1 */ if (p == 0)
    {
        cout << "Vetv 1" << endl;
        double y123 = pow(-2.0*q, 1.0 / 3.0);
        cout << "x1= " << y123 - b / (3 * a) << endl;
        cout << "x2= " << y123 - b / (3 * a) << endl;
        cout << "x3= " << y123 - b / (3 * a) << endl;
    }
 
 
    /* 2 */   if (q == 0)
    {
        if (p <= 0)
        {
            cout << "Vetv 2" << endl;
            double y1 = 0;
            double y2 = sqrt(-3.0*p);
            double y3 = -1.0*sqrt(-3.0*p);
            cout << "x1 = " << y1 - b / (3 * a) << endl;
            cout << "x2 = " << y2 - b / (3 * a) << endl;
            cout << "x3 = " << y3 - b / (3 * a) << endl;
        }
        else
        {
            cout << "Vetv 2.1" << endl;
            double y1 = 0;
            double y23 = sqrt(3.0*p);
            cout << "x1 = " << y1 - b / (3 * a) << endl;
            cout << "x2 = i * " << y23 - b / (3 * a) << endl;
            cout << "x3 = -i * " << y23 - b / (3 * a) << endl;
        }
    }
 
    /* 3 */ if ((p<0) && (q != 0) && ((q*q + p*p*p) <= 0))
    {
        cout << "Vetv 3" << endl;
        r = sign(q)*sqrt(fabs(p));
        phi = acos(q / pow(r, 3.0));
 
        double y1 = -2.0*r*cos(phi / 3.0);
        double y2 = 2.0*r*cos((PI - phi) / 3.0);
        double y3 = 2.0*r*cos((PI + phi) / 3.0);
        cout << "x1 = " << y1 - b / (3 * a) << endl;
        cout << "x2 = " << y2 - b / (3 * a) << endl;
        cout << "x3 = " << y3 - b / (3 * a) << endl;
    }
 
    /* 4 */if ((p<0) && (q != 0) && ((q*q + p*p*p) >0))
    {
        cout << "Vetv 4" << endl;
        r = sign(q)*sqrt(fabs(p));
        phi = acosh(q / pow(r, 3.0));
 
        double y1 = -2.0*r*cosh(phi / 3.0);
        double y23_1 = r*cosh(phi / 3.0);
        double y23_2 = fabs(sqrt(3.0)*r*sinh(phi / 3.0));
        cout << "x1 = " << y1 - b / (3 * a) << endl;
        cout << "x2 = " << y23_1 - b / (3 * a) << " + i * (" << y23_2 << ")" << endl;
        cout << "x3 = " << y23_1 - b / (3 * a) << " - i * (" << y23_2 << ")" << endl;
    }
 
    if ((p>0) && (q != 0))
    {
        cout << "Vetv 5" << endl;
        r = sign(q)*sqrt(fabs(p));
        phi = asinh(q / (r*r*r) );
        double y1 = -2.0*r*sinh(phi / 3.0);
        double y23_1 = r*sinh(phi / 3.0);
        double y23_2 = fabs(sqrt(3.0)*r*cosh(phi / 3.0));
        cout << "x1 = " << y1-b/(3*a) << endl;
        cout << "x2 = " << y23_1 - b / (3 * a) << " + i * (" << y23_2 << ")" << endl;
        cout << "x3 = " << y23_1 - b / (3 * a) << " - i * (" << y23_2 << ")" << endl;
    }
 
    system("pause");
}

@Frip · 10.12.2014, 20:29 **[ТС]**

Сошлись, при ваших данных, я вводил, к примеру a=1, b=2, c=3, d=5. Ответы не сходятся и всё тут, отличаются в десятой части, может погрешности где какие-нибудь выходят.

Не по теме:

Автор А.В. Крячков, И.В. Сухинина, В.К. Томшин. Алгоритм решения рассматривается в 6 главе сие учебника. Как вложение прикрепить не могу, размер не позволяется, вот ссылка на файлообменник http://rghost.ru/59548055

@olgashat · 10.12.2014, 23:14

x1 = −1.84373427789807
x2 = −0.0781328610509654 + i ∙ 1.64492637759997
x3 = −0.0781328610509654 − i ∙ 1.64492637759997
с сайта решение при Ваших коэффициентах и мое решение

@Frip · 12.12.2014, 17:00 **[ТС]**

ДА! Спасибо вам большое! У меня просто стояли коэф. a=1, b=2, c=3, d=4 , и из-за этого все не сходилось, моя вина

Добавлено через 9 минут
Причём очень глупая ошибка*

Новые блоги и статьи Все статьи Все блоги /
ИИ не может найти нужный язык в списке Supersumestria 05.10.2026 Я ему даю вот такое изображение и прошу найти и подчеркнуть немецкий язык. Возвращает он вот это: https:/ / i. **********/ vqBWLe2. png Нужную строчку в 3й колонке просто выдумал. . Это. . .	Новая последняя моя музыка в SUNO zorxor 05.10.2026 Здравствуйте, дорогие мои друзья! С большой радостью я хотел бы представить вам свою новую последнею музыку, которую сгенерировала мне по моей просьбе нейросеть SUNO. С уважением, zorxor. Это. . .	Программный домашний кинотеатр russiannick 27.09.2026 Сподобился на программный домашний кинотеатр. В качестве ЯВУ по традиции выбрал js. В помощники взял Яндекс-Алису. Было создано три зала на разные интересы. исторические и ретро сериал Хичкок. . .	Беседа с ИИ о программистах, недопускающих к созданию и правке кода генеративные ИИ и причины этого zorxor 21.09.2026 Раньше я радовался или получал некоторые эмоции, пусть небольшие, но всё же, от самого процесса написания кода, рекомпиляции и запуска, видя постепенное развитие программы и прочее. А теперь лень. . .
Мобильное приложение ColorStep pavlinmavlin 17.09.2026 Реализовал приложение Красный, Зеленый, Синий в Unity3d + c#. Название изменил на ColorStep. Приложение прошло модерацию и теперь доступно для скачивания. Делал его сам, шаг за шагом — и вот,. . .	Запрет дублирования строк в табличной части Maks 13.09.2026 Реализация из решения ниже выполнена на нетиповом справочнике "Нормы ТО" с табличной часть "Виды ТО", разработанного в КА2, со следующими реквизитами: - ВидТО (СправочникСсылка. ВидыТО); - ВидГСМ. . .	Скрипты Tampermonkey для CyberForum, ChatGPT, Claude и пр. Jin X 06.09.2026 Скрипты Tampermonkey для CyberForum, ChatGPT, Claude и пр. Работая с форумом и нейросетями в браузере часто хочется что-то подкорректировать или добавить какого-то функционала. Ниже прикреплён. . .	Программа опроса у.з. расходомера SLS-720F Argus19 02.09.2026 Программа опроса у. з. расходомера SLS-720F Программа опрашивает один раз в минуту три ультразвуковых расходомера SLS-720F через интерфейс RS-485 по протоколу Modbus RTU. Опрашиваются регистры. . .

@Frip 3 / 3 / 1 Регистрация: 04.11.2013 Сообщений: 285
	09.12.2014, 20:59 [ТС]
	ааа, вы про это? Я это изменил, результат я огласил уже вроде, ответы всё-равно не те 0

@Frip 3 / 3 / 1 Регистрация: 04.11.2013 Сообщений: 285
	10.12.2014, 20:29 [ТС]
	Сошлись, при ваших данных, я вводил, к примеру a=1, b=2, c=3, d=5. Ответы не сходятся и всё тут, отличаются в десятой части, может погрешности где какие-нибудь выходят. Не по теме: Автор А.В. Крячков, И.В. Сухинина, В.К. Томшин. Алгоритм решения рассматривается в 6 главе сие учебника. Как вложение прикрепить не могу, размер не позволяется, вот ссылка на файлообменник http://rghost.ru/59548055 0

@olgashat 91 / 74 / 81 Регистрация: 07.12.2014 Сообщений: 303
	10.12.2014, 23:14
	x1 = −1.84373427789807 x2 = −0.0781328610509654 + i ∙ 1.64492637759997 x3 = −0.0781328610509654 − i ∙ 1.64492637759997 с сайта решение при Ваших коэффициентах и мое решение Изображения 1

@Frip 3 / 3 / 1 Регистрация: 04.11.2013 Сообщений: 285
	12.12.2014, 17:00 [ТС]
	ДА! Спасибо вам большое! У меня просто стояли коэф. a=1, b=2, c=3, d=4 , и из-за этого все не сходилось, моя вина Добавлено через 9 минут Причём очень глупая ошибка* 0

При выводе числа выводится #INF

Решение