Runge-Kutta法求解轧制力