代码之家  ›  专栏  ›  技术社区  ›  Anthony

R中的Lsoda函数

0

Anthony  · 技术社区  · 9 年前

mo=function(t,x,m){
if(t%%2==0){
   a=m
}else{
   a=0
}
dx=rep(0,2)
dx[1]=-a*x[1]
dx[2]=a*x[1]
res=dx
return(list(dx))
}
xs=c(10,0)
ti=1:100
m=1
data1=as.data.frame(lsoda(xs,ti,mo,m))