Fixed issue with imaginary parts in RkSolver.
Boundary function is taken into account while calculating k1, ..., k4.