Controling the global error
Before the start of the computation process, the user can set
the maximum admitted global error of the approximate solution. The
implicit value is e=10-5, but this value can be changed. Depending
on this value some computations are made in order to establish the
corresponding level of the stepsize and the estimated computation time.
Example for the method DIRK4 applied to the system B1.