虽然存在相关的问题和答复,但我无法将它们翻译成我自己的问题,所以在这里:
ratedegrad<-0.012
releaser<-0.0016666667
times2<-seq(from=1, to=240, by=1)
state2<-c(degrad=300,titre=101,released=13.7)
dtitre2<-function(y,t,parms) {
with(as.list(y), {
dX2=-ratedegrad*titre
dY2=(-releaser*titre)-(ratedegrad*titre)
dZ2= releaser*titre
list(c(abs(dX2),dY2,dZ2)) } )
}
outtitre2<-ode(y=state2,times=times2,func=dtitre2,parms=NULL)
plot(outtitre2)
我得到的错误信息是: 由func()(153)返回的导数的数量必须等于初始条件向量(3)的长度
谢谢你!