!
! MackeyGlass.f90
!
! Vector field function for: MackeyGlass
! These functions are to be used with DDE_SOLVER_M.
!
! This file was generated by the program VFGEN (Version:2.1.0)
! Generated on  3-Jun-2007 at 12:24
!
!
! The lags are: {tau}

MODULE DEFINE_MackeyGlass_DDEs

    IMPLICIT NONE
    INTEGER, PARAMETER :: NEQN=1, NLAGS=1
    ! Parameters
    DOUBLE PRECISION a, b, tau

CONTAINS

    SUBROUTINE MackeyGlass_ddes(t,x_,Zlags_,vf_)
    ! Arguments
    DOUBLE PRECISION :: t
    DOUBLE PRECISION, DIMENSION(NEQN) :: x_, vf_
    DOUBLE PRECISION, DIMENSION(NEQN,NLAGS) :: Zlags_
    ! Local variables
    DOUBLE PRECISION delayedx
    DOUBLE PRECISION x
    ! State variables
    x          = x_(1)
    ! Expressions
    delayedx = Zlags_(1,1)
    ! The vector field
    vf_(1) = -b*x+delayedx*a/(delayedx**10.0+1.0)

    RETURN
    END SUBROUTINE MackeyGlass_ddes

    SUBROUTINE MackeyGlass_history(t,x_)
    DOUBLE PRECISION :: t
    DOUBLE PRECISION, DIMENSION(NEQN) :: x_

    x_(1) = 2.0000000000000000e-02*t+5.0000000000000000e-01

    RETURN
    END SUBROUTINE MackeyGlass_history

END MODULE DEFINE_MackeyGlass_DDEs
