Zdefiniuj . Następnie oryginalne równanie można zapisać jako
gdzie
Chcemy znaleźć analitycznie, rozwiązując (1).α:=r−1
(1)dx1dz=[x1+βr−αz][−1x0dx0dz]
(2)x20=[1+βα]−[βα1(r−αz)3]
x1
Różnicowanie (2) daje
Dlatego też ostatni czynnik w (1) można zapisać jako
Teraz, używając (2),
Możemy teraz zapisać (4) jako
(3)dx0dz=−12x03β(r−αz)4
(4)−1x0dx0dz=32β2x20(r−αz)4
x20(r−αz)4β2=1β2(1+βα)(r−αz)4−βα(r−αz)=r−αzαβ2[(α+β)(r−αz)3−β]
(5)−1x0dx0dz=(3αβ2)(βr−αz)[1(α+β)(r−αz)3−β]
Wtyczka (5) do (1), aby uzyskać
Dokonaj zmiany zmiennych, definiując . Następnie
Możemy teraz napisać (6) jako
(6)dx1dz=3αβ2[x1+βr−αz](βr−αz)[1(α+β)(r−αz)3−β]
y:=r−αzdx1dz=dx1dydydz=−αdx1dy
(7)dx1dy=−3β2[x1+βy][β(α+β)y4−βy]
lub
Możesz teraz rozwiązać ODE przy użyciu standardowych metod, np .
Http://tutorial.math.lamar.edu/Classes/DE/Linear.aspx . Rozwiązanie to złożona mieszanka pierwiastków kwadratowych, kłód i arktanów.
(8)dx1dy=−3β22[x1(α+β)y4−βy+β(α+β)y5−βy2]