Solving Index-1 DAEs in MATLAB and Simulink. - Lawrence F. Shampine; Mark W. Reichelt; Jacek Kierzenka (SIAM Rev., 1999)

A BVP solver based on residual control and the Maltab PSE. - Jacek Kierzenka; Lawrence F. Shampine (ACM Trans. Math. Softw., 2001)

Non-negative solutions of ODEs. - Lawrence F. Shampine; Skip Thompson; Jacek Kierzenka; George D. Byrne (Appl. Math. Comput., 2005)