Solving Abel's ODEs of the First Kind - Maple Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Mathematics : Differential Equations : Classifying ODEs : First Order : odeadvisor/Abel

Solving Abel's ODEs of the First Kind

 

Description

Examples

Description

• 

The general form of Abel's equation of the first kind is given by:

Abel_ode := diff(y(x),x)=f3(x)*y(x)^3+f2(x)*y(x)^2+f1(x)*y(x)+f0(x);

Abel_ode:=ⅆⅆxyx=f3xyx3+f2xyx2+f1xyx+f0x

(1)
  

where f3(x), f2(x), f1(x) and f0(x) are arbitrary functions.

  

See Differentialgleichungen, by E. Kamke, p. 24. There is as yet no general solution for this ODE. For Abel's equation of the second kind, see Abel2A and Abel2C.

• 

The most general method available at the moment to solve Abel ODEs seems to be the method of "Abel's invariant", described in E. Kamke, p. 26, as sub-method (g) due to M. Chini. The invariant of an Abel equation with f2=0 is the following quantity:

Abel_invariant := -1/27/f3(x)^4*(-diff(f0(x),x)*f3(x)+f0(x)*diff(f3(x),x)+
3*f0(x)*f3(x)*f1(x))^3/f0(x)^5;

Abel_invariant:=127ⅆⅆxf0xf3x+f0xⅆⅆxf3x+3f0xf3xf1x3f3x4f0x5

(2)
  

If the invariant does not depend on x, then the equation can be solved directly.

  

For an Abel equation with f2<>0, the f2 term can be removed by using the following transformation:

y(x)=u(x)-f2(x)/3/f3(x);

yx&equals;ux13f2xf3x

(3)
  

The invariant can then be calculated as in the previous case. Note that if an Abel ODE has a constant invariant, then any other Abel ODE obtained from it by a transformation of the form

{y(x)=G(t)*u(t)+H(t), x=F(t)};

x&equals;Ft&comma;yx&equals;Gtut&plus;Ht

(4)
  

will also have a constant invariant (that is, is also solvable by this method).

  

The method "Chini" (see ?odeadvisor,Chini), also due to Chini, generalizes this method of the constant invariant for Abel ODEs.

Examples

Abel_ode&DifferentialD;&DifferentialD;xyx&equals;f3xyx3&plus;f2xyx2&plus;f1xyx&plus;f0x

Abel_ode:=&DifferentialD;&DifferentialD;xyx&equals;f3xyx3&plus;f2xyx2&plus;f1xyx&plus;f0x

(5)

Abel_invariant1&DifferentialD;&DifferentialD;xf0xf3x&plus;f0x&DifferentialD;&DifferentialD;xf3x&plus;3f0xf3xf1x327f3x4f0x5

Abel_invariant:=127&DifferentialD;&DifferentialD;xf0xf3x&plus;f0x&DifferentialD;&DifferentialD;xf3x&plus;3f0xf3xf1x3f3x4f0x5

(6)

yx&equals;uxf2x3f3x

yx&equals;ux13f2xf3x

(7)

yx&equals;Gtut&plus;Ht&comma;x&equals;Ft

x&equals;Ft&comma;yx&equals;Gtut&plus;Ht

(8)

withDEtools&comma;odeadvisor

odeadvisor

(9)

withPDEtools&comma;dchange

dchange

(10)

odeadvisorAbel_ode

_Abel

(11)

1) An example of an Abel ODE having a constant invariant solved using the related scheme:

ODE&DifferentialD;&DifferentialD;xyx&equals;112x&plus;27x3&plus;27x3yx2&plus;18x2yx&plus;27yx3x3&plus;27x2yx2&plus;9xyx&plus;127x3&colon;

ansdsolveODE

ans:=yx&equals;1929RootOf81&int;_Z1841_a327_a&plus;27&DifferentialD;_a&plus;x&plus;3_C1x3x3x

(12)

Any "linear transformation" of ODE will also be solved by the same method. For example:

TR_LINyx&equals;Gtut&plus;Ht&comma;x&equals;Ft

TR_LIN:=x&equals;Ft&comma;yx&equals;Gtut&plus;Ht

(13)

ODE_pdchangeTR_LIN&comma;ODE&comma;u&comma;t&colon;

ans_pdsolveODE_p&comma;ut

ans_p:=ut&equals;199HtFt29RootOf81&int;_Z1841_a327_a&plus;27&DifferentialD;_a&plus;Ft&plus;3_C1Ft&plus;3Ft&plus;3GtFt

(14)

2) A case for which the solving method is known: f0(x) = f1(x) = 0, and diff(f3(x)/f2(x),x)=a*f2(x).

In this case, one can proceed as follows:

odesubsf1x&equals;0&comma;f0x&equals;0&comma;Abel_ode

ode:=&DifferentialD;&DifferentialD;xyx&equals;f3xyx3&plus;f2xyx2

(15)

First introduce r(t) and t as new variables using:

ITRyx&equals;f2trtf3t&comma;x&equals;t

ITR:=x&equals;t&comma;yx&equals;f2trtf3t

(16)

new_odedchangeITR&comma;ode&comma;rt&comma;t&colon;

Now, introduce the condition on the derivative of f3(t)/f2(t):

constraint&DifferentialD;&DifferentialD;tf3tf2taf2t&equals;0

constraint:=&DifferentialD;&DifferentialD;tf3tf2tf3t&DifferentialD;&DifferentialD;tf2tf2t2af2t&equals;0

(17)

and simplify new_ode with regard to this relation:

new_ode2simplifynew_ode&comma;constraint&comma;&DifferentialD;&DifferentialD;tf3t

new_ode2:=rtf2t3a&plus;f2t&DifferentialD;&DifferentialD;trtf3tf3t2&equals;f2t3rt3&plus;f2t3rt2f3t2

(18)

This ODE is separable.

odeadvisornew_ode2

_separable

(19)

3) Rewrite in "normal form" (no square term in the RHS) when: f3(x)=1/x, f2(x)=1/x, f1(x)=0, f0(x)=4

odeevalsubsf3x&equals;1x&comma;f2x&equals;1x&comma;f1x&equals;0&comma;f0x&equals;4&comma;Abel_ode

ode:=&DifferentialD;&DifferentialD;xyx&equals;yx3x&plus;yx2x&plus;4

(20)

First of all, Abel's ODEs of the first kind can be rewritten in normal form (which is sometimes useful) by making the appropriate change of variables. The transformation is of a general type. After introducing

wx&ExponentialE;&int;f1xf2x23f3x&DifferentialD;x

wx:=&ExponentialE;&int;f1x13f2x2f3x&DifferentialD;x

(21)

the following transformations (where {x,y(x)} = old vars; {t,r(t)} = new vars) will yield the desired normal form:

trt&equals;&int;f3x&apos;w&apos;x2&DifferentialD;x&comma;rt&equals;13yxf3x&plus;f2x3&ExponentialE;&int;f1x1f2x23f3x&DifferentialD;xf3x

tr:=t&equals;&int;f3xwx2&DifferentialD;x&comma;rt&equals;133yxf3x&plus;f2x&ExponentialE;&int;f1x13f2x2f3x&DifferentialD;xf3x

(22)

The transformation equations required for this case are obtained from the general transformation tr (above) as follows:

TRevalsubsf3x&equals;1x&comma;f2x&equals;1x&comma;f1x&equals;0&comma;f0x&equals;4&comma;tr

TR:=t&equals;32x2&sol;3&comma;rt&equals;13x4&sol;33yxx&plus;1x

(23)

ITRsolveTR&comma;yx&comma;x

ITR:=x&equals;32RootOf2_Z2t&plus;3t&comma;yx&equals;13RootOf2_Z2t&plus;33rtRootOf2_Z2t&plus;3

(24)

and the change of variables is implemented as follows:

new_odedchangeITR&comma;ode&comma;t&comma;rt&comma;&apos;known&apos;&equals;indetsode&comma;&apos;unknown&apos;&colon;

new_ode2simplifyop1&comma;mapallvalues&comma;solvenew_ode&comma;&DifferentialD;&DifferentialD;trt

new_ode2:=&DifferentialD;&DifferentialD;trt&equals;11818t3rt3t261t243t3

(25)

Finally, the normal form can be made explicit as follows:

collectnew_ode2&comma;rt&comma;factor

&DifferentialD;&DifferentialD;trt&equals;rt3118t261t&plus;243t3

(26)

See Also

DEtools

odeadvisor

dsolve

quadrature

linear

separable

Bernoulli

exact

homogeneous

homogeneousB

homogeneousC

homogeneousD

homogeneousG

Chini

Riccati

Abel

Abel2A

Abel2C

rational

Clairaut

dAlembert

sym_implicit

patterns

odeadvisor,types

 


Download Help Document

Was this information helpful?



Please add your Comment (Optional)
E-mail Address (Optional)
What is ? This question helps us to combat spam