1

虽然存在相关的问题和回复，但我无法将它们转化为我自己的问题，所以这里是：

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)

谢谢你！

4

0 回答 0