Lagerstrom ODEs
|
Description
|
|
•
|
The general form of the Lagerstrom ODE is given by the following:
|
>
|
Lagerstrom_ode := diff(y(x),x,x)= -k*diff(y(x),x)/x-epsilon*y(x)*diff(y(x),x);
|
| (1) |
|
See Rosenblat and Shepherd, "On the Asymptotic Solution of the Lagerstrom Model Equation".
|
|
|
Examples
|
|
The second order Lagerstrom ODE can be reduced to a first order ODE of Abel type once the system succeeds in finding one polynomial symmetry for it (see ?symgen):
>
|
|
>
|
|
| (2) |
>
|
|
| (3) |
From which, giving the same indication directly to dsolve, you obtain the reduction of order
>
|
|
| (4) |
For the structure of the solution above see ?ODESolStruc. Reductions of order can also be tested with odetest
>
|
|
| (5) |
The reduced ODE is of Abel type and can be selected using the mouse, or as follows
>
|
|
| (6) |
>
|
|
| (7) |
|
|
See Also
|
|
DEtools, odeadvisor, dsolve, and ?odeadvisor,<TYPE> where <TYPE> is one of: quadrature, missing, reducible, linear_ODEs, exact_linear, exact_nonlinear, sym_Fx, linear_sym, Bessel, Painleve, Halm, Gegenbauer, Duffing, ellipsoidal, elliptic, erf, Emden, Jacobi, Hermite, Lagerstrom, Laguerre, Liouville, Lienard, Van_der_Pol, Titchmarsh; for other differential orders see odeadvisor,types.
|
|