Algorithm 407: DIFSUB for solution of ordinary differential equations D2 | Synapse