next up previous

2.3.4 VODE Example     continued...

This code is followed by two tabulations of the solution at , , , : the first for a supplied Jacobian (MF=21) and the second for an internally generated Jacobian (MF=22). Both methods produce acceptably accurate solutions, but more work is required for the internally-generated Jacobian case. The quantity SUM is also printed out. Since , SUM should have the value ``1.'' Note that this condition is necessary for accuracy, but not for sufficient. The output is followed by a combined plot of , , and vs. t for .

Figure 8: Output from the VODE solution to Robertson's chemical kinetics problem.

(See exercise sets 3 and 3.)