A variable time-step-size code for advection-diffusion-reaction PDEs

The main point for the solution of the large-dimensional ODEs is not to solve for the stage values of the Radau method until convergence until convergence, but only giving a very few iterations and take as advancing solution the latter stage value computed.