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()