Maxima 中也提供了相应的求解函数 desolve(),desolve()函数既可以求解ODE 方程,也可以求解ODE方程组。函数的基本形式如下。 desolve (eqn, y) desolve ([eqn_1, ..., eqn_n], [y_1, ..., y_n]) 这里待解函数不能只写变量名(例如y),而需要明确写出对自变量的依赖关系 如果提供了足够的初值条件,再用的desolve()函数求解时积分常数自然就可以确定了。 (%i3) atvalue(y(x),x=0,1); desolve(eq5,y(x)); atvalue('diff(y(x),x),x=0,-1); desolve(eq5,y(x)); atvalue 说明 desolve 函数提供的初值必须是x=0 处的。 ic1 不能用来直接求解 desolve 函数的结果。必须要人为的处理一下结果的形式。这一点上确实不方便。
求解微分方程 desolve函数 实例1 实例2 实例3 实例4 求解有条件的微分方程 微分方程显示隐式解 未找到显式解决方案时查找隐式解决方案 求微分方程级数解 为具有不同单边限制的函数指定初始条件 (特解) 练习题 desolve函数 S = dsolve(eqn)求解微分方程eqn,其中eqn是符号方程。
I - gamma * I dR <- gamma * I list(c(dS, dI, dR)) }) } image.png 接下来可以就交给R做参数估计和模型求解 library(deSolve
我们来使用最新的数据计算 R0: library(deSolve) R0 <- function(days){ historylist <- suppressMessages( df$jsondata