include working RK4 code and local potential calculation