C++ version of DUMKA3 is a template function dumka3.cpp. Idea to use templates was taken from IML++ library. In order to use templates you need to have classes VectorFunction and Vector

This is the first C++ version of dumka3, and as everything new it will be changed and modified in a future, and your suggestions/complains/questions would be very welcome.

Both fortran version of dumka3.for and C++ version of dumka3.cpp require additional function which determine maximal step of explicit Euler method. It can be evaluated by value cou=2/M, where M is upper boundary of maximal eigen value. For example if you solve ODE y'=-5*y, then take cou=2/5 or something smaller, for example 2/5.1. You can use also norm of jacobian of the right hand side function.