Solving.Differential.Equations.in.R笔记1

例1：y'=ry(1-y/K)

          y(0)=2

对应的R代码为：

library(deSolve)
#parameters and initial values
r<-1
K<-10
yini<-2
#the function
derivs<-function(t,y,parms){
                 return(list(r*y*(1-y/K)))
        }
#the independent variables
times<-seq(from=0,to=20,by=0.2)
#to calculate y values
out<-ode(y=yini,times=times,func=derivs,parms=NULL)
#the results
out

参考文献：

Karline S, Thomas P, Setzer R W. Solving Differential Equations in R[M]. Springer Publishing Company, Incorporated, 2012.

转载于:https://www.cnblogs.com/lww1993/p/7489826.html