Kubinės lygties sprendimas be kompleksinių skaičių

Iš testwiki.
Pereiti į navigaciją Jump to search
https://www.journals.vu.lt/LMR/article/view/14970/13988


Yra įrodyta, kad kubinės ir ketvirto laipsnio lygtys gali būti išspręstos be kompleksinių skaičių ir išvestinių panaudojimo (randami tik realieji sprendiniai). Bus pateikti atitinkami algoritmai.

Kubinės lygtys

Bendra forma kubinės lygties yra
ax3+bx2+cx+d=0(a0).(1)
Padalinus šią lygtį iš a0, gaunama lygtis, taip pat vadinama bendros formos:
x3+bx2+cx+d=0.(2)
(2) lygtis gali būti perdaryta/suprastina įvedus naują kintamąjį y, kad x=y+k. Skaičius k pasirenkamas taip, kad (2) lygtis neturėtų antro laipsnio nario (bx2). Iš tikro, po įstatymo x=y+k, tiktai iš x3 ir bx2 gaunamas y kvadratas (padaugintas iš konstantos).
x3+bx2+cx+d=0.(2)
(y+k)3+b(y+k)2+c(y+k)+d=0,
(y3+3y2k+3yk2+k3)+b(y2+2yk+k2)+c(y+k)+d=0,
(y3+3yk2+k3)+[3y2k+by2]+b(2yk+k2)+c(y+k)+d=0.
Igriko kvadratas išnyks, kai [3y2k+by2]=0 arba 3k+b=0. Tada k=b/3.
Tai galima gauti ir taip:
x3+bx2=x2(x+b)=(y+k)2(y+k+b)=(y2+2yk+k2)(y+k+b);
sandaugoje (y2+2yk+k2)(y+k+b) tik ky2+by2+2y2k=3ky2+by2 turi y2. Todėl y2 išnyks, jeigu 3k+b=0, t. y. jeigu k=b/3.
Gauta lygtis vadinama redukuota (nežinomasis vėl yra x):
x3+px+q=0.(3)

Sprendinių skaičius

Kubinė lygtis (3) visada turi bent vieną sprendinį. Funkcija f(x)=x3+px+q yra teigiamai su pakankamai dideliais x ir neigiama su dideliais neigiamais x. Todėl f(x) funkcijos grafikas kerta Ox ašį.
Kubinė lygtis gali turėti vieną, du ar trys sprendinius – pavyzdžiui, lygtys
(x1)3=0,(x1)2(x2)=0,(x1)(x2)(x3)=0
turi atitinkamai sprendinius {1},{1,2},{1,2,3}.
Daugiau nei 3 sprendinius kubinė lygtis negali turėti. Tarkim, kad lygtis (3) turi sprendinį α. Tai reiškia, kad α3+pα+q=0. Dabar (3) lygtį lengva išskaidyti (faktorizuoti)
x3+px+q=x3+px+qα3pαq=x3α3+pxpα=(xα)(x2+xα+α2+p).
Taigi kitus sprendinius lygties (3) rasime išsprenę kvadratinę lygtį (nuo nežinomojo x)
x2+xα+α2+p=0.(4)
Lygtis (4) gali turėti 0 arba 1, arba 2 sprendinius, kurie nėra α. Šie sprendiniai taip pat bus sprendiniais (3) lygties.
[x3+px+q=0.(3)]
Matome, kad žinant bent vieną (3) lygties spreninį, galima rasti visus kitus (3) lygties sprendinius iš (4) lygties. Štai kodėl labai lengva spręsti kubinę lygtį, kuri turi racionalų sprendinį. Beje, galima lengvai surasti bet kokio laipsnio lygties sprendinius, kai jos koeficientai yra racionalieji skaičiai. Pavyzdžiui, jeigu lygties koeficientai yra sveikieji skaičiai ir pirmas koeficientas yra 1 (tokią formą galima suteikti bet kokiai lygčiai su racionaliais koeficientais), tada racionalieji sprendiniai gali būti tik sveikieji skaičiai – teigiami ir neigiami dalikliai laisvojo nario (kubinėje lygtyje (3) laisvasis narys yra q; žr. čia). Todėl apsimoka sprendžiant bet kokią lygtį su racionaliaisiais koeficientais, pradėti lygties sprendimą nuo racionaliųjų sprendinių ieškojimo.

Atvejis, kai p>0.

Kubinė lygtis (3) gali būti toliau suprastinta, padarius koeficiento p0 modulį lygų 3 (kai p=0, tada lygtis (3) tampa x3=q, turinti sprendinį x=q3).
Vėl mes panaudojame paprastą, tiesinį keitinį x=ky ir tinkamai parinksime k. Mūsų (3) lygtis konvertuojama į
k3y3+pky+q=0,y3+py/k2+q/k3=0.
Atveju p>0, pasirenkame k taip:
p/k2=3,p/3=k2,k=p/3.
Todėl su tokiu k, koeficientas prie y tampa lygus 3.
x=ky=p/3y,
q/k3=q/(p/3)3=qp3/(39)=3qpp/3=33qpp.
Pažymėkime
m=33q2pp.
Tada mes turime lygtį (gautą iš lygties y3+py/k2+q/k3=0):
y3+py/k2+q/k3=0,
y3+3y2m=0,(5)
kurioje p/k2=3,33qpp=2m.
Dabar mes bandysime rasti sprendinius, parinkę keitinį y=z1/z:
(z1/z)3+3(z1/z)2m=0,
(z33z21/z+3z1/z21/z3)+3(z1/z)2m=0,
(z33z+3/z1/z3)+3z3/z2m=0,
(z31/z3)2m=0,
z31/z32m=0,|z3
z612mz3=0,
z62mz3=1,
z62mz3+m2=1+m2,
(z3m)2=1+m2,
z3m=±1+m2,
z3=m±1+m2,
z=m±1+m23.
Abi šios z reikšmės duoda tą patį sprendinį (5) lygties.
y=z1/z=m+1+m231m+1+m23=
=m+1+m23m1+m23m+1+m23m1+m23=m+1+m23m1+m23m2(1+m2)23=
=m+1+m23m1+m23m2(1+m2)3=m+1+m23m1+m231=
=m+1+m23+m1+m23;


y=z1/z=m1+m231m1+m23=
=m1+m23m+1+m23m+1+m23m1+m23=m1+m23m+1+m23m2(1+m2)23=
=m1+m23m+1+m23m2(1+m2)3=m1+m23m+1+m231=
=m1+m23+m+1+m23.


Taigi, abiais atvejais (su dviem skirtingom z reikšmėm) gavome tą patį sprendinį
y=m+1+m23+m1+m23.
Tai skiriasi nuo sprendinio iš PDF failo, kuris yra toks:
y=m2+1+m3+m2+1m3.
Įsitikinkime, kad (5) lygtis neturiu daugiau sprendinių. Pagal (4) lygtį (iš (5) lygties), dabar turime:
[x2+xα+α2+p=0.(4)]
[y3+3y2m=0,(5)]
y2+yα+α2+3=0,(5.1)
kur
α=m+1+m23+m1+m23.
Diskriminantas (5.1) lygties lygus
D=α24(α2+3)=3α212<0.
Kadangi diskriminantas neigiamas, (5) lygtis daugiau sprendinių neturi. Todėl (3) lygtis, atveju p>0, turi tik vieną sprendinį. Mes rasime jį perėję nuo y prie x=yp/3.
Rasime kam lygus x (kuris yra sprendinys (3) lygties).
m=33q2pp.
y=m+1+m23+m1+m23.
x=ky=p/3y=p/3(m+1+m23+m1+m23)=
=p/3(33q2pp+1+(33q2pp)23+33q2pp1+(33q2pp)23)=
=p/3(33q2pp+1+93q24p2p3+33q2pp1+93q24p2p3)=
=p/3(33q2pp+1+27q24p33+33q2pp1+27q24p33)=
=p/3(33q2pp+4p3+27q24p33+33q2pp4p3+27q24p33)=
=(p33333q2pp+p3334p3+27q24p33+p33333q2ppp3334p3+27q24p33)=
=pp3333q2pp+pp334p3+27q24p33+pp3333q2pppp334p3+27q24p33=
=q2+1334p3+27q243+q21334p3+27q243=
=q2+4p3+27q24273+q24p3+27q24273=
=q2+p327+q243+q2p327+q243.

Atvejis, kai p<0.

Sunkiau išspręsti (3) lygtį, kai p yra neigiamas. Ir vėl panaudojame keitinį x=ky. Dabar lygtyje
y3+py/k2+q/k3=0
pasirenkame k, kad p/k2=3,
p/3=k2,
k=p/3.
Pošaknyje yra teigiamas skaičius, nes p<0.
q/k3=q/(p/3)3=qp3/(39)=3qpp/3=33qpp;
33qpp=q/k3=2m;
m=33q2pp.
Pažymėję laisvajį narį per -2m, gauname lygtį
y33y2m=0.(6)
Panaudojame keitinį y=z+1/z:
(z+1/z)33(z+1/z)2m=0,
(z3+3z21/z+3z1/z2+1/z3)3(z+1/z)2m=0,
(z3+3z+3/z+1/z3)3(z+1/z)2m=0,
(z3+3z+3/z+1/z3)3z3/z2m=0,
(z3+1/z3)2m=0,
z3+1/z32m=0,|z3
z6+12mz3=0,
z62mz3=1,
z62mz3+m2=m21,
(z3m)2=m21.(7)
Lygtis (7) nevisada turi sprendinius ir reikia išnagrinėti trys atvejus: m2>1, m2=1 ir m2<1.

Atvejis p<0, m*m>1 (unikalus sprendinys).

Jeigu m2>1, tada
[(z3m)2=m21.(7)]
z3m=±m21,
z3=m±m21,
z=m±m213,
ir abi šios z reikšmės duoda tą patį sprendinį:
y=z+1/z=m+m213+1m+m213=
=m+m213+mm213mm213m+m213=m+m213+mm213m2(m21)23=
=m+m213+mm213m2(m21)3=m+m213+mm2131=
=m+m213+mm213;
y=z+1/z=mm213+1mm213=
=mm213+m+m213m+m213mm213=mm213+m+m213m2(m21)23=
=mm213+m+m213m2(m21)3=mm213+m+m2131=
=mm213+m+m213.
Gavome tas pačias y reikšmes.
Gausime x reikšmę ((3) lygties sprendinį):
k=p/3;
m=33q2pp;
x=ky=p/3y=p/3(mm213+mm213)=
=p/3(33q2pp+(33q2pp)213+33q2pp(33q2pp)213)=
=p/3(33q2pp+93q24p2(p)13+33q2pp93q24p2p13)=
=p/3(33q2pp+27q24p313+33q2pp27q24p313)=
=p/3(33q2pp+4p3+27q24p33+33q2pp4p3+27q24p33)=
=(p33333q2pp+p3334p3+27q24p33+p33333q2ppp3334p3+27q24p33)=
=pp3333q2pp+pp334p3+27q24p33+pp3333q2pppp334p3+27q24p33=
=q2+1334p3+27q243+q21334p3+27q243=
=q2+4p3+27q24273+q24p3+27q24273=
=q2+p327+q243+q2p327+q243.


Įsitikinkime, kad šis sprendinys yra unikalus (vienintelis). Šio sprendinio modulis yra didesnis už 2:
y2=(z+1/z)2=z2+2+1/z2=z22+1/z2+4=(z1/z)2+4>4.
Vadinasi, |y|>2 ir iš (6) lygties
2m=y33y=y(y23)>y(43)=y,
2m>y;
2m>y>2(kaiy>0),
m>1(kaiy>0).
Jeigu pvz., y=-3, tai
2m>3(kaiy<0),
2m>2>y(kaiy<0),
m>1(kaiy<0).
Bet mūsų atveju turi būti m2>1, tačiau to nėra, kai m>1, nes jeigu pavyzdžiui, m=-0.9, tai m2=(0.9)2=0.81<1, kas neatitinka sąlygos m2>1. Todėl teigiamos y reikšmės patenka į sąlygą m2>1. Kitaip tariant, jeigu y yra neigiamas, tai jo reikšmės mažesnės už -2 (y<-2, kai y<0) ir 2m gali būti lygu -2 (nes 2m>y). Tada gaunasi, kad 2m=-2, m=-1. Ir m2=(1)2=1, bet tai neatitinka atvejo m2>1.
Bet jeigu, pavyzdžiui, y=-3, tai iš nelygybės 2m>y, gauname
2m>3,
m>1.5.
Ir gaunasi, kad sąlyga m2>1 tenkinama, nes (1.49)2=2.2201>1.
Vadinasi, sąlyga m2>1 netenkinama tik kai 2m=2,m=1. Kalbant dar tiksliau, su bet kokiom y reikšmėm, m gali įgyti vienintelę reikšmę m=1, kuri nepriklauso šiam atvejui p<0,m2>1.
Pataisymas. Jeigu m=-0.5, tai tada m>-1 (kai y<0). Arba 2m=-1>-2>y (kai y<0). Todėl, kai y neigiamas, m reikšmės gali būti segmente [-1; 0]. Arba 1m0, kai y<0. O kai y>0, tada m reikšmės gali būti intervale (1;). Tokios neigiamos m reikšmės (1m0) tenkina visus y. Su neigiamom didelėm absoliučiu dydžiu (modulis didelis) y reikšmėm, m neigiamos reikšmės bus irgi didesnės absoliučiu dydžiu (nelygybė 2m>y bus teisinga, kai y<0). O neigiamos a reikšmės iš segmento [-1; 0] visada tenkins nelygybę 2a>y (kai y<0), bet m reikšmės iš šio segmento nepriklauso šiam atvejui (p<1,m2>1).


Kiti sprendiniai galėtų būti gauti iš (4) lygties, kuri dabar atrodo taip:
[x2+xα+α2+p=0.(4)]
[y33y2m=0.(6)]
y2+yα+α23=0(kur|α|>2).(4.6)
Jos diskriminantas
D=α24(α23)=3α2+12=3(4α2)<0
yra neigiamas ir (4.6) lygtis neturi sprendinių (su α=mm213+m+m213). Todėl (6) lygtis, atveju p<0,m2>1 turi vienintelį sprendinį α (α=m+m213+mm213).

Atvejis p<0, m*m=1 (du sprendiniai).

Jeigu m2=1, tada iš (7) lygties
[(z3m)2=m21.(7)]
turime z3=m.
Kai m=1, tai z=1 ir
y = z + 1/z = 2.
(6) lygtis tranformuojasi į
y33y2=0,
kuri gali būti išskaidyta nepasinaudojant (4) lygtim:
y33y2=0,(y2)(y2+2y+1)=0,(y2)(y+1)2=0.
Taigi lygtis (6) (ir (3)) šiuo atveju turi du sprendinius ((6) lygties sprendiniai yra 2 ir -1, kai p<0, m*m=1). Grįžtame prie x.
x=ky=p/3y.
x1=ky1=p/32=2p/3;
x2=ky2=p/3(1)=p/3.
Kai m=-1, tai z=-1 ir
y = z + 1/z = -1 + 1/(-1) = -2.
(6) lygtis tampa
y33y+2=0,(y+2)(y22y+1)=0,(y+2)(y1)2=0.
Ir taip pat turi du sprendinius (kaip ir (3)). Grįžtame prie x.
x=ky=p/3y.
x1=ky1=p/3(2)=2p/3;
x2=ky2=p/31=p/3.
p<0.
x1 ir x2 yra (3) lygties sprendiniai.


x312x+16=0
pavyzdys, kurios p=-12<0.
Tada
x1=2p/3=2(12)/3=212/3=24=22=4;
x2=p/3=(12)/3=12/3=4=2.
Šie sprendiniai yra tokie patys kaip ir pateiktoje nuorodoje.
Patikriname:
x1312x1+16=0,
(4)312(4)+16=0,
64+48+16=0;
x2312x2+16=0,
23122+16=0,
824+16=0.
Vadinasi sprendiniai teisingi.
m=33q2pp=33162(12)(12)=33162412=33162434=
=316244=316384=1684=1616=1.
Kažkodėl šitame pavyzdyje m=1, bet x1 ir x2 paimti iš atvejo, kai m=-1.


[y3+py/k2+q/k3=0
pasirenkame k, kad p/k2=3,
p/3=k2,
k=p/3.
Pošaknyje yra teigiamas skaičius, nes p<0.
q/k3=q/(p/3)3=qp3/(39)=3qpp/3=33qpp;
33qpp=q/k3=2m;
m=33q2pp.]
Gali būti, kad gavome m=1, o ne m=-1, nes
q/k3=q/(p/3)3=qp3/(39)=3qp2p/3=33qp2p;
tada
33qp2p=q/k3=2m;
m=33q2p2p=
=33162(12)2(12)=
=3316214412=
=331621234=
=3162124=
=3162122=
=316434=1.
Dabar ir gavome m=1.

Atvejis p<0, m*m<1 (trys sprendiniai).

Jeigu m2<1, tai lygtis (7) neturi sprendinių ir keitinys y=z+1/z neras sprendinius. Bet čia padeda trigonometrija: pakeitimas y=2z pakeičia (6) lygtį į tokią:
[y33y2m=0.(6)]
(2z)332z2m=0,
8z36z2m=0,
4z33zm=0.
Darome keitinį z=cosφ:
4cos3φ3cosφm=0,cos(3φ)m=0,cos3φ=m.
Ši lygtis turi sprendinius, kai |m|<1 (arba m2<1). Kosinusas yra lyginė funkcija, todėl su neigiamom ir teigiamom φ reikšmėm gaunamas tas pats kosinuso rezultatas. Todėl
cos3φ=m,
cos(±3φ)=m,
±3φ=arccosm,
±(±3φ)=3φ=±arccosm,
3φ=±arccosm,
3φ2kπ=±arccosm,
3φ=2kπ±arccosm,
φ=2kπ3±13arccosm;
z=cosφ=cos(23kπ±13arccosm),
y=2z=2cosφ=2cos(23kπ±13arccosm),
k - sveikasis skaičius.
Remdamiesi redukcijos formulėmis, suprantame, kad su daugiau nei trim skaičiaus k reikšmėm z reikšmės kartosis, nes kubinė lygtis negali turėti daugiau nei trys sprendinius. Todėl pasirenkame k=0,1,2. Tada
y1=2z1=2cos(13arccosm),
y2=2z2=2cos(23π±13arccosm),
y3=2z3=2cos(43π±13arccosm).
Grįžtame prie x.
x=ky=p/3y;
x1=y1p/3=2p/3cos(13arccosm),
x2=y2p/3=2p/3cos(23π±13arccosm),
x3=y3p/3=2p/3cos(43π±13arccosm).
Kubinė lygtis (3) yra išspręsta.
m=33q2pp.


x319x+30=0,
kurios p=19, q=30.
Jos sprendiniai yra 2, 3 ir 5.
Apskaičiuosime juos pagal ką tik išvestas formules.
m=33q2pp=33302(19)(19)=9032(19)19=453(19)19=
= 45*3^0.5 /(19*19^0.5) = 0.9411150958093732309574814856304.
Arba kitas variantas:
m=33q2p2p=33302(19)2(19)=903236119=4531919=
= -45*3^0.5 /(19*19^0.5) = -0.9411150958093732309574814856304.
Tada
arccosm=arccos(0.9411150958093732309574814856304)= 0.34488276150211935210875728756448.
Arba
arccosm=arccos(0.9411150958093732309574814856304)= 2.796709892087673886353886095715.


Apskaičiuosime x1.
Su pirma m reikšme:
x1=2p/3cos(13arccosm)=2(19)/3cos(13arccos(0.94111509580937323))=
=219/3cos(130.3448827615021193521)=219/3cos(0.11496092050070645)=
= 2*(19/3)^0.5 * cos(0.11496092050070645) = 2*(19/3)^0.5 *0.99339926779878285489956379038765 = 5 (Windows 10 kalkuliatoriaus duoda tokią tikslią reikšmę "5", įstačius pajuodintą reikšmę).
Su antra m reikšme:
x1=2p/3cos(13arccosm)=2(19)/3cos(13arccos(0.94111509580937323))=
=219/3cos(132.796709892087673886353886095715)=219/3cos(0.93223663069589129545)=
= 2*(19/3)^0.5 * cos(0.93223663069589129545) = 2*(19/3)^0.5 *0.59603956067926971293973827423259 = 3 (Windows 10 kalkuliatoriaus duoda tokią tikslią reikšmę "3", įstačius pajuodintą reikšmę).
Su pirma m reikšme, atsakymas 5 yra be minuso ženklo. O su antra m reikšme, atsakymas 3 yra duotos lygties sprendinys.


Toliau apskaičiuosime x2.
Su pirma m reikšme:
x2=2p/3cos(23π±13arccosm)=2(19)/3cos(23π±13arccos(0.3448827615021193521))=
=219/3cos(23π±0.11496092050070645)=
= 2*(19/3)^0.5 * cos(2*pi/3 +0.11496092050070645070291909585483) =
= 2*(19/3)^0.5 * cos(2.0943951023931954923084289221863 +0.11496092050070645070291909585483) =
= 2*(19/3)^0.5 * cos(2.2093560228939019430113480180411) = 2*(19/3)^0.5 * (-0.59603956067926971293973827423256) =
= -2.9999999999999999999999999999999.
O jeigu po pi paimsime ne pliuso, bet minuso ženklą, tai gausime:
= 2*(19/3)^0.5 * cos(2.0943951023931954923084289221863 -0.11496092050070645070291909585483) =
= 2*(19/3)^0.5 * cos(1.9794341818924890416055098263315) = 2*(19/3)^0.5 * (-0.39735970711951314195982551615502) =
= -1.9999999999999999999999999999998.
Čia kosinusas skaičiuojamas su reikšmėm didesnėm nei pi/2 =~ 1.57, todėl galimos kai kurios klaidos. Bendrai, negalimos klaidos, nes čia kosinusas, o ne arkkosinusas.
Su antra m reikšme:
x2=2p/3cos(23π±13arccosm)=2(19)/3cos(23π±13arccos(0.94111509580937323))=
=219/3cos(23π±0.932236630695891295451)=
= 2*(19/3)^0.5 * cos(2*pi/3 +0.93223663069589129545129536523833) =
= 2*(19/3)^0.5 * cos(2.0943951023931954923084289221863 + 0.93223663069589129545129536523833) =
= 2*(19/3)^0.5 * cos(3.0266317330890867877597242874246) = 2*(19/3)^0.5 * (-0.99339926779878285489956379038764) =
= -5 (įstačius į Windows 10 kalkuliatorių pajuodintą tekstą).
O jeigu po pi paimsime ne pliuso, bet minuso ženklą, tai gausime:
= 2*(19/3)^0.5 * cos(2.0943951023931954923084289221863 - 0.93223663069589129545129536523833) =
= 2*(19/3)^0.5 * cos(1.162158471697304196857133556948) = 2*(19/3)^0.5 * 0.39735970711951314195982551615505 =
= 2 (įstačius į Windows 10 kalkuliatorių pajuodintą tekstą).
Su antra m reikšme iškart gavome du teisingus tos lygties sprendinius.


Toliau apskaičiuosime x3.
Su pirma m reikšme:
x3=2p/3cos(43π±13arccosm)=2(19)/3cos(43π±13arccos(0.3448827615021193521))=
=219/3cos(43π±0.11496092050070645)=
= 2*(19/3)^0.5 * cos(4*pi/3 +0.11496092050070645070291909585483) =
= 2*(19/3)^0.5 * cos(4.1887902047863909846168578443727 +0.11496092050070645070291909585483) =
= 2*(19/3)^0.5 * cos(4.3037511252870974353197769402275) = 2*(19/3)^0.5 * (-0.39735970711951314195982551615506) =
= -2.
O jeigu po pi paimsime ne pliuso, bet minuso ženklą, tai gausime:
= 2*(19/3)^0.5 * cos(4.1887902047863909846168578443727 - 0.11496092050070645070291909585483) =
= 2*(19/3)^0.5 * cos(4.0738292842856845339139387485179) = 2*(19/3)^0.5 * (-0.59603956067926971293973827423257) =
= -2.9999999999999999999999999999999.
Su antra m reikšme:
x3=2p/3cos(43π±13arccosm)=2(19)/3cos(43π±13arccos(0.94111509580937323))=
=219/3cos(43π±0.932236630695891295451)=
= 2*(19/3)^0.5 * cos(4*pi/3 +0.93223663069589129545129536523833) =
= 2*(19/3)^0.5 * cos(4.1887902047863909846168578443727 + 0.93223663069589129545129536523833) =
= 2*(19/3)^0.5 * cos(5.121026835482282280068153209611) = 2*(19/3)^0.5 * 0.39735970711951314195982551615505 =
= 2 (įstačius į Windows 10 kalkuliatorių pajuodintą tekstą).
O jeigu po pi paimsime ne pliuso, bet minuso ženklą, tai gausime:
= 2*(19/3)^0.5 * cos(4.1887902047863909846168578443727 - 0.93223663069589129545129536523833) =
= 2*(19/3)^0.5 * cos(3.2565535740904996891655624791344) = 2*(19/3)^0.5 * (-0.99339926779878285489956379038764) =
= -5 (įstačius į Windows 10 kalkuliatorių pajuodintą tekstą).
Su antra m reikšme iškart gavome du teisingus tos lygties sprendinius.


Su antra m reikšmę gavome visus teisingus lygties x319x+30=0 sprendinius (2, 3. -5). O su pirma m reikšme gavome sprendinius su priešingais ženklais (-2, -3, 5; dėl to visi jie neteisingi).

Pavyzdys.

Imkim konkrečią lygtį formoje (3):
x33x4=0.(8)
Bus spręsta iš naujo be teorijos (teorijos gimsta iš konkrečių pavyzdžių sakoma toje nuorodoje). Reikia galvoti apie formulę tipo x33x. Tada ateina į mintis formulė
(a+b)3=a3+3a2b+3ab2+b3=a3+b3+3ab(a+b)
ir dalinis atvejis
(a±1/a)=a3±1/a3±3(a±1/a).
Paskutinė formulė sako, kad (8) lygtyje verta panaudoti keitinį x=y+1/y. Tada lygtis pavirtsta į
(y+1/y)33(y+1/y)4=0,
y3+1/y3+3(y+1/y)3(y+1/y)4=0,
y3+1/y34=0.
Tai yra kvadratinė lygtis y3 atžvilgiu (y6+14y3=0), todėl ją galima išspręsti (galima net nenaudoti kvadratinės lygties sprendimo formulės):
y64y3+1=0,
(y32)23=0,
(y32)2=3,
y32=±3,
y3=2±3,
y=2±32.
Iš skyrelio Atvejis p<0, m*m>1 (unikalus sprendinys). panaudojus formulę
y=z+1/z=m+m213+mm213;
x=y+1/y=2+2213+22213=2+33+233
(abu sprendiniai vienodi su skirtingais y).
Vis tiek reikia įsitikinti, kad nėra kitų sprendinių (su x reikšme x=y+1/y negalima rasti sprendinius intervale (-2; 2), nes iš to, kad y visada teigiamas, |y+1/y|=|y|+1/|y|=(|y|1/|y|)2+22).
(8) lygtis gali būti išskaidyta padalinus iš tiesinio polinomo x(2+33+233), bet tai būtų sunkus darbas. Lengviau pasinaudoti (4) formule, kuri tampa tokia:
[x2+xα+α2+p=0.(4)]
x2+xα+α23=0,
kur α=2+33+233.
Šios lygties diskriminantas
d=α24(α23)=3α2+12
yra neigiamas, nes
α2=(y+1/y)2=(y1/y)2+4>4.
Todėl (8) lygtis turi vienintelį sprendinį (x=2+33+233).
x=2+33+233=
= (2 + 3^0.5)^(1/3) + (2 - 3^0.5)^(1/3) = 2.1958233454456471528327992055497.
2.1958233454456471528327992055497^3 - 3*2.1958233454456471528327992055497 - 4 = -1.1291287212984656298036039312783e-31 =
1.129128721298465629803603931278310310.

Kardano formulė

Dabar aišku, kad nagrinėtus atvejus galima taiktyti lygčiai (3).
Atvejis 1: p327+q24>0.
Taikomas keitinys
x=y3p/(3y3).
Surasti y, tada grįžti prie x:
x=q2+p327+q243+q2p327+q243.
Tada
x3+px+q=0.(3)
(y3p/(3y3))3+p(y3p/(3y3))+q=0,
y3y23p/(3y3)+3y3p2/(9y23)p3/(27y)+p(y3p/(3y3))+q=0,
yy3p+p2/(3y3)p3/(27y)+p(y3p/(3y3))+q=0,
yy3p+p2/(3y3)p3/(27y)+py3p2/(3y3)+q=0,
yp3/(27y)+q=0,
y2p3/27+qy=0,
y2+qyp3/27=0;
D=q24(p3/27);
y1=q+q24(p3/27)2=q2+q24+p3/27,
y2=qq24(p3/27)2=q2q24+p327.
x=y3p/(3y3),
x1=y13p/(3y13)=
=q2+q24+p3273p/(3q2+q24+p3273)=
=q2+q24+p3273p3q2+q24+p3273=
=q2+q24+p3273pq2q24+p32733q2+q24+p3273q2q24+p3273=
=q2+q24+p3273pq2q24+p32733(q2)2(q24+p327)23=
=q2+q24+p3273pq2q24+p32733q24(q24+p327)3=
=q2+q24+p3273pq2q24+p32733p3273=
=q2+q24+p3273pq2q24+p32733(p3)=
=q2+q24+p3273pq2q24+p3273p=
=q2+q24+p3273+q2q24+p3273.
Nesunku matyti, kad įstačius y2 į y3p/(3y3) bus gautas toks pats atsakymas.


Atvejis 2: p327+q24=0.
Kardano formulė duoda sprendinį
α=q2+p327+q243+q2p327+q243=
=q23+q23=2q23=8q23=4q3.
Tada (4) lygtis duoda sprendinį
[x2+xα+α2+p=0.(4)]
x2+xα+α2+p=x24q3x+16q23+p;
D=16q234(16q23+p)=316q234p;
x1=(4q3)+316q234p2=4q/83+(316q234p)/4=
=q/23+316q2/643p=q/23+3q2/43p=
Čia https://lt.wikibooks.org/wiki/Matematika/Kubinė_lygtis pagal (22) formulę
q23=3q2p.(22)
2p=3qq23,
2p=3q2/23,
p=3q2/(28)3=3q2/163.
Tada
x1=q/23+3q2/43p=q/23+3q2/43+3q2/163=
=q/23+3q2/43+343q2/43=q/23+433q2/43+3q2/4343=
=q/23+(43+1)3q2/4343=q/23+3(43+1)q23163=
=q/23+3(43+1)q2343=q/23+q33(43+1)43=
=q/23+q33(43+1)43=
Kvadratinėje šaknyje neigiamas skaičius (3(43+1)), todėl toliau skaičiuot beprasmiška.
Kai q=0 (tai tada ir p=0) ir tada sprendinys yra x=0. Kai q0, gaunami du sprendiniai, tik ne šituo budu.
  • x33x+2=0.
Šios lygties sprendiniai yra -2 ir 1.
x1=(4q3)+316q234p2=4q3+316q234p2=
=423+316(3)234(3)2=
=83+31693+122=
=2+31443+122=
= 1 + 0.5* (-3 * 144^(1/3) +12)^0.5 = 1 + 0.5* (-3 * 5.241482788417793214283322560884 +12)^0.5 = Invalid input, nes skliaustuose neigiamas skaičius.
x2=(4q3)316q234p2=4q/83(316q234p)/4=
= 1 - 0.5* (-3 * 144^(1/3) +12)^0.5 = 1 - 0.5* (-3 * 5.241482788417793214283322560884 +12)^0.5 = Invalid input, nes skliaustuose neigiamas skaičius.
x3=α=4q3=423=2.


Atvejis 3: p327+q24<0.
Suprantama, p<0. Taikant keitį x=4|p|/3cosϕ, gaunasi lygtis
4cos3ϕ3cosϕ=q27/(2p|p|),
cos3ϕ=q27/(2p|p|).(A)
Kadangi
(q27)2<(2p|p|)2,
27q2<4p3,
tai (A) lygties dešinėje pusėje modulio reikšmė mažesnė už 1. Radus φ, randami 3 sprendiniai lygties (3):
α1=x=4|p|/3cos(13arccos(3q3/2p|p|),
α2=4|p|/3cos(23π+13arccos(3q3/2p|p|),
α3=4|p|/3cos(43π+13arccos(3q3/2p|p|).
Juos surasti galima taip:
x=ky,
k=p/3,
y=2z,
x=2kz=2p/3z.
x3+px+q=0,
x3+px+q=8(p/3)3z3+p2p/3z+q=8p3p/3z3+2pp/3z+q=0,
4p3p/3z3+pp/3z=q/2,
4p/3z3+3p/3z=3q2p,
4z3+3z=33q2pp,
4z33z=33q2pp.


Tolesnis skaičiavimas vyksta taip:
z=cosφ,
x=2p/3z=2p/3cosφ.
p<0, todėl
x=2|p|/3z=2|p|/3cosφ.
4z33z=33q2pp,
4cos3φ3cosφ=33q2pp,
cos3φ=33q2pp;
3φ=arccos(33q2pp),
φ=13arccos(33q2pp).
Ir p<0.
y=2z=2cosφ,
x=2kz=2p/3z=2p/3cosφ=2p/3cos(13arccos(33q2pp)).

Ketvirto laipsnio lygtis

Bendra ketvirto laipsnio lygtis yra
ax4+bx3+cx2+dx+e=0(a0).
Ji gali buti padalinta iš a, todėl užtenka nagrinėti lygtį
x4+ax3+bx2+cx+d=0.
Po pakeitimo x=y-a/4, eliminuojamas narys trečio laipsnio. Todėl užtenka mokėti spręsti lygtį
x4+bx2+cx+d=0.(9)
Jeigu ketvirto laipsnio lygtis turi racionalius sprendinius (arba bent vieną racionalų sprendinį α), ją lengva išspręsti: kairioji pusė gali būti išskaidyta vėl, atskirus xα ir belieka išspręst kūbinę lygtį. Bus nagrinėjama lygtis (9), kai ji neturi racionaliųjų sprendinių, nes taip sunkiau.
Jeigu c=0, lygtis (9) yra bikvadratinė, ir pasirinkus naują kintamąjį x2=y, ji tampa kvadratine. Beje, lengva išskaidyti bikvadratinį trinarį x4+bx2+d. Jeigu diskriminantas b2/4d0, tada
x4+bx2+d=(x2+b/2)2(b2/4d)=
=(x2+b/2+b2/4d)(x2+b/2b2/4d).
Jeigu b2/4d<0, tada d>b2/4 (d teigiamas), d>b/2,2db/2>0, ir
x4+bx2+d=(x2+d)2x2(2db/2)=(x2+d+x2db/2)(x2+dx2db/2).
Jeigu c0, tai kairę pusę (9) lygties visada įmanoma užrašyti kaip kvadratų skirtumą:
(x2+m)2p(x+n)2=0(p>0).(9.1)
(x2+m)2p(x+n)2=x4+2mx2+m2p(x2+2xn+n2)=
=x4+2mx2+m2px22xpnpn2=x4+(2mp)x22xpn+m2pn2.
Kairioji pusė (9.1) lygties gali būti išskaidyta į du kvadratinius trinarius,
(x2+m)2p(x+n)2=(x2+m+p(x+n))(x2+mp(x+n))=
=(x2+px+m+pn))(x2px+mpn)),
ir išsprendžiamos dvi kvadratinės lygtys.
Bus nustatyta kaip randamos m, p, n reikšmės. Lygtyje
(x2+m)2p(x+n)2=x4+bx2+cx+d(9.5)
nuėmus skliaustus ir palyginus koeficientus prie x2,x1,x0, gaunama lygčių sistema m, p, n reikšmėms rasti.
x4+2mx2+m2p(x2+2xn+n2)=x4+bx2+cx+d,
x4+2mx2+m2px22xpnpn2=x4+bx2+cx+d,
x4+(2mp)x22xpn+m2pn2=x4+bx2+cx+d;
(2mp)x2=bx2,
2mp=b,
2pn=c,
m2pn2=d;
2mp=b,2pn=c,m2pn2=d.(10)
Kadangi c0, iš antros (10) lygties p0 ir tada n=c/(2p). Iš pirmos (10) lygties m=(b+p)/2. Po įstatymo į trečią lygtį n ir m reikšmes, trečia lygtis tampa
m2pn2=(b+p)2/4p(c/(2p))2=(b+p)2/4pc2/(4p2)=(b+p)2/4c2/(4p)=d.
Po padauginimo iš 4p, gaunama kubinė lygtis:
(b+p)2/4c2/(4p)=d,
p(b+p)2c2=4dp,
p(b2+2pb+p2)c24dp=0,
pb2+2p2b+p3c24dp=0,
p3+2p2b4dp+pb2c2=0.
f(p)=p3+2p2b4dp+pb2c2.
Ši lygtis turi teigiamą sprendinį (tai yra raktas ketvirto laipsnio lygties sprendimui). Tikrai, funkcija f(p), kai p=0, yra lygi f(0)=c2, t. y. neigiama (kai c0), ir su pakankamai didele p reikšme f(p) funkcija yra teigiama (nes funkcijos ženklą su didelėm p reikšmėm apsprendžia narys p3). Todėl funkcijos f(p) grafikas kerta p (arba x) ašį iš intervalo (0; ), nes f(0)=-c*c<0, o f()>0.
Kadangi galima išspręsti kūbinę lygtį (nuo p), galima surasti teigiamą sprendinį p0, gaunant sprendinį (10) lygties
(m,p,n)=(b/2+p0/2,p0,12c/p0).
Lygtyje-lygybėje (9.5) ikso reikšmės turi būti tokios, kad (x2+m)2=0 ir p(x+n)2=0. Tada su šia ar šiom x reikšmėm ir x4+bx2+cx+d=0. Tai vienas variantas. Kito būdo rodos nėra. Todėl galimas variantas, kad randami ne visi ketvirto laipsnio lygties (9) sprendiniai.
Nes lygtyje
(x2+m)2p(x+n)2=x4+bx2+cx+d,(9.5)
(x2+m)2p(x+n)2 gali būti daugiau už 0 su tam tikrom x reikšmėm ir gali būti mažiau už nulį su kitom x reikšmėm. Ir su kai kuriom x reikšmėm reiškinys (x2+m)2p(x+n)2 gali būti lygus nuliui, nors nei (x2+m)2, nei p(x+n)2 neprivalo būti lygūs nuliui.
Ketvirto laipsnio lygtis gali turėti nuo 1 iki 4 sprendinių arba nei vieno sprendinio (pavyzdžiui, lygtis x4+1=0 neturi [realių] sprendinių).
Išspręsti penkto laipsnio lygtį bendru atveju yra neįmanoma (many try, but now all lieing dead).
Tačiau, svarbiausias teiginys pagrindinės algebros teoremos išlieka teisingas: kiekvienas n-to laipsnio polinomas išskaidomas į tiesinius ir kvadratinius daugiklius.
Papildymas. Jeigu (x2+m)2p(x+n)2=0, tai
(x2+m)2=p(x+n)2,
x2+m=p(x+n),
x2+m=px+pn,
x2px+mpn=0.
Išsprendus šią kvadratinę lygtį bus rastos dvi lygties (9) šaknys.
Be to, dėmuo (x2+m)2 visada yra teigiamas skačius su visom x reikšmėm, o dėmuo p(x+n)2 su neigiamom x reikšmėm yra mažesnis teigiamas skaičius (jei n yra teigiamas skaičius), nei su teigiamom x reikšmėm (jei n yra teigiamas skaičius). Todėl gal taip ir prisideda ketvirtas sprendinys. Kai trečias sprendinys ateina iš sąlygos, kad (x2+m)2=0 ir p(x+n)2=0.