迭代循环不一致，完全迭代

问：

我正在尝试编写一个步长为 175 的迭代循环，但每次运行程序时，循环都会以不同的长度迭代。每次我运行时，它都会按步长增加最终迭代。不知道发生了什么。

对于上下文，我正在将 ODE 与 Runge-Kutta 子程序集成。

do i= 175,17500,175
call rkf45(real(i,wp),ts,te,y,esp1,esp2,h,hmin,175._wp,dy,flag,neqn)
q(i)=y(1)
ts=te
te=te+175.0_wp
hmax=te-ts
write(*,*)i,q(i)
enddo

连续 3 次运行的最终输出

11200   257478.60158947381
13650   431868.79246271809
17500   808062.17349709466

不知道要尝试什么，尝试在前后将 i 初始化为 0。