我有一个微分方程系统。我用数值 dsolve 解决它:
sol := dsolve({First, Fourth, Second, Third, iA(0) = 0, iB(0) = 0, theta(0) = 0, (D(theta))(0) = 0}, numeric, vars,maxfun=100000):
First,Second,Third,Fourth---是微分方程,然后有初始条件,vars---是我想要得到的一组变量(关于时间)。maxfun - dsolve 为计算导数而进行的最大评估次数。
所以,我需要检索计算变量的最大时间值。我尝试设置范围参数,但解决时间太长(我需要调用 dsolve 一百次甚至更多次)。使用 maxfun=1e5 我会在三四秒内得到结果。
有什么技巧可以在这里完成吗?也许有一种方法可以在衍生品评估之间设置时间步长?
在此先感谢并为我的英语不好感到抱歉。