added array_algebra which gives about 20% performance gain for the rk4 stepper however, const correctness is a problem there....