Initial Value Problems

The code GAM numerically solves solves first order ordinary differential equations, either stiff or nonstiff in the form y'=f(x,y), with a given initial condition. The code GAMD is a generalization of GAM for the solution of Differential Algebraic Equations of index less than or equal to 3 in the form M y' = f(x,y), with a given initial condition. By Francesca Mazzia.