3,809 research outputs found
NumGfun: a Package for Numerical and Analytic Computation with D-finite Functions
This article describes the implementation in the software package NumGfun of
classical algorithms that operate on solutions of linear differential equations
or recurrence relations with polynomial coefficients, including what seems to
be the first general implementation of the fast high-precision numerical
evaluation algorithms of Chudnovsky & Chudnovsky. In some cases, our
descriptions contain improvements over existing algorithms. We also provide
references to relevant ideas not currently used in NumGfun
Accurate calculation of the solutions to the Thomas-Fermi equations
We obtain highly accurate solutions to the Thomas-Fermi equations for atoms
and atoms in very strong magnetic fields. We apply the Pad\'e-Hankel method,
numerical integration, power series with Pad\'e and Hermite-Pad\'e approximants
and Chebyshev polynomials. Both the slope at origin and the location of the
right boundary in the magnetic-field case are given with unprecedented
accuracy
Symbolic computation of exact solutions expressible in hyperbolic and elliptic functions for nonlinear PDEs
Algorithms are presented for the tanh- and sech-methods, which lead to
closed-form solutions of nonlinear ordinary and partial differential equations
(ODEs and PDEs). New algorithms are given to find exact polynomial solutions of
ODEs and PDEs in terms of Jacobi's elliptic functions.
For systems with parameters, the algorithms determine the conditions on the
parameters so that the differential equations admit polynomial solutions in
tanh, sech, combinations thereof, Jacobi's sn or cn functions. Examples
illustrate key steps of the algorithms.
The new algorithms are implemented in Mathematica. The package
DDESpecialSolutions.m can be used to automatically compute new special
solutions of nonlinear PDEs. Use of the package, implementation issues, scope,
limitations, and future extensions of the software are addressed.
A survey is given of related algorithms and symbolic software to compute
exact solutions of nonlinear differential equations.Comment: 39 pages. Software available from Willy Hereman's home page at
http://www.mines.edu/fs_home/whereman
Complex oscillations in the delayed Fitzhugh-Nagumo equation
Motivated by the dynamics of neuronal responses, we analyze the dynamics of
the Fitzhugh-Nagumo slow-fast system with delayed self-coupling. This system
provides a canonical example of a canard explosion for sufficiently small
delays. Beyond this regime, delays significantly enrich the dynamics, leading
to mixed-mode oscillations, bursting and chaos. These behaviors emerge from a
delay-induced subcritical Bogdanov-Takens instability arising at the fold
points of the S-shaped critical manifold. Underlying the transition from
canard-induced to delay-induced dynamics is an abrupt switch in the nature of
the Hopf bifurcation
On elliptic solutions of the quintic complex one-dimensional Ginzburg-Landau equation
The Conte-Musette method has been modified for the search of only elliptic
solutions to systems of differential equations. A key idea of this a priory
restriction is to simplify calculations by means of the use of a few Laurent
series solutions instead of one and the use of the residue theorem. The
application of our approach to the quintic complex one-dimensional
Ginzburg-Landau equation (CGLE5) allows to find elliptic solutions in the wave
form. We also find restrictions on coefficients, which are necessary conditions
for the existence of elliptic solutions for the CGLE5. Using the investigation
of the CGLE5 as an example, we demonstrate that to find elliptic solutions the
analysis of a system of differential equations is more preferable than the
analysis of the equivalent single differential equation.Comment: LaTeX, 21 page
Construction of Special Solutions for Nonintegrable Systems
The Painleve test is very useful to construct not only the Laurent series
solutions of systems of nonlinear ordinary differential equations but also the
elliptic and trigonometric ones. The standard methods for constructing the
elliptic solutions consist of two independent steps: transformation of a
nonlinear polynomial differential equation into a nonlinear algebraic system
and a search for solutions of the obtained system. It has been demonstrated by
the example of the generalized Henon-Heiles system that the use of the Laurent
series solutions of the initial differential equation assists to solve the
obtained algebraic system. This procedure has been automatized and generalized
on some type of multivalued solutions. To find solutions of the initial
differential equation in the form of the Laurent or Puiseux series we use the
Painleve test. This test can also assist to solve the inverse problem: to find
the form of a polynomial potential, which corresponds to the required type of
solutions. We consider the five-dimensional gravitational model with a scalar
field to demonstrate this.Comment: LaTeX, 14 pages, the paper has been published in the Journal of
Nonlinear Mathematical Physics (http://www.sm.luth.se/math/JNMP/
Computation of saddle type slow manifolds using iterative methods
This paper presents an alternative approach for the computation of trajectory
segments on slow manifolds of saddle type. This approach is based on iterative
methods rather than collocation-type methods. Compared to collocation methods,
that require mesh refinements to ensure uniform convergence with respect to
, appropriate estimates are directly attainable using the method of
this paper. The method is applied to several examples including: A model for a
pair of neurons coupled by reciprocal inhibition with two slow and two fast
variables and to the computation of homoclinic connections in the
FitzHugh-Nagumo system.Comment: To appear in SIAM Journal of Applied Dynamical System
- …