PROGRAM pulsate

    asymp
    pulse

SUBROUTINES

    asymp
      dosim
      modetrap
    bump
      rkf
      splintl
    consl
    dosim
    fehl
      rkfcow
      rkfliq
    grind
      rkf
      splintl
    init
    modetrap
      pdisc()
      zbrent()
    perg
    pulse
      init
      bump
      perg
      grind
    rkf
      rkfs
    rkfcow
      splintl
    rkfliq
      splintl
    rkfs
      rkfcow
      rkfliq
      fehl
    splintl
      consl

FUNCTIONS

    pdisc()
    zbrent()