Question

Mathematica中利用NDsolve怎么打印它所选的步长或时间点

当用NDsolve求解一个以时间t为积分常数的常微分方程时，怎样能够输出它用龙格库塔法求解时所选的步长，或者时间点。最好是选定了下个步长，就能把选定的步长或者下一个时间点输出

Answer 1

values = {}; Monitor[
 NIntegrate[Sqrt[x (1 - x)], {x, 0, 1}, 
  EvaluationMonitor :> (Pause[0.025]; AppendTo[values, x];)], 
 Quiet@ListPlot[values]]

上面是一个看NIntegrate函数计算过程的例子，修改一下就可以看NDSolve函数了。NDSolve也有EvaluationMonitor选项。