3,809 research outputs found

    NumGfun: a Package for Numerical and Analytic Computation with D-finite Functions

    Get PDF
    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

    Get PDF
    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

    Get PDF
    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

    Full text link
    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

    Full text link
    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

    Full text link
    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

    Get PDF
    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 ϵ\epsilon, 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
    • …
    corecore