!
! sddfn.f90
!
! Vector field function for: sddfn
! These functions are to be used with DDE_SOLVER_M.
!
! This file was generated by the program VFGEN (Version:2.4.0)
! Generated on 10-Jul-2008 at 13:08
!
!
! The lags are: { t-beta}

MODULE DEFINE_sddfn_DDEs

    IMPLICIT NONE
    INTEGER, PARAMETER :: NEQN=1, NLAGS=1

CONTAINS

    SUBROUTINE sddfn_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 beta, delayedy
    DOUBLE PRECISION y
    ! State variables
    y          = x_(1)
    ! Expressions
    beta = -2.0**(1.0/2.0)+y
    delayedy = Zlags_(1,1)
    ! The vector field
    vf_(1) = (1.0/2.0)*( t+1.0)**-(1.0/2.0)*delayedy

    RETURN
    END SUBROUTINE sddfn_ddes

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

    x_(1) = 1.0

    RETURN
    END SUBROUTINE sddfn_history

    SUBROUTINE sddfn_beta(t,x_,bval_)
    ! Arguments
    DOUBLE PRECISION :: t
    DOUBLE PRECISION, DIMENSION(NEQN) :: x_
    DOUBLE PRECISION, DIMENSION(NLAGS) :: bval_
    ! Local variables
    DOUBLE PRECISION beta, delayedy
    DOUBLE PRECISION y

    ! State variables
    y          = x_(1)
    ! Expressions
    beta = -2.0**(1.0/2.0)+y
    bval_(1) = beta
    END SUBROUTINE sddfn_beta

END MODULE DEFINE_sddfn_DDEs