improve integrator efficiency

eliminate a numpy array allocation inside rhs
update example fig2a script to fix inconsistency