THE SUBROUTINE Bsub

The form of this subroutine depends on the choices of Info(12) and Info(14). When Info(12)=0 or Info(14)=0, a dummy form of Bsub suffices  

When Info(12) 0 and Info(14) = l, Bsub will be called by DDAPLUS, and may be organized as in the following pseudocode. 

Subroutine Bsub (t, Nstvar, U, fpj, Jspar, Rpar, Ipar, Irange, Imodel, 

This subroutine calculates the nonzero elements of fpj= df /dOj, for use by DDAPLUS in constructing a column of the matrix B t). Bsub is executed only if MAIN sets lnfo(12) 0 and lnfo(14)=l before calling DDAPLUS. If dEjdOj is nonzero, leave lres=0 to activate DDAPLUS s finite-difference algorithm for this column of 

Double precision t, U( ), fpj( ), Rpar(0 ), Integer Nstvar, Jspar, lpar(0 r), Irange, Ires, 

---