The method of lines translates an evolutionary PDE into a large system of ODEs. This system exhibits stiffness which in consequence requires specialized numerical algorithms. Our simulations proved that the Radau IIA method is a good choice.
The source code of Radau IIA method (RADAU5) well as other related software is freely available at.....
