Displaying 2 results from an estimated 2 matches for "pileggi".
Did you mean:
leggi
2011 Oct 03
0
deSolve - Function daspk on DAE system - Error (Vince)
...need to look at the formulation itself.
By the way, there is a special R-mailinglist that deals with this type
of problems:
r-sig-dynamic-models at r-project.org
Hope this helps,
Karline
---------------------
Original message:
Date: Sat, 1 Oct 2011 20:20:10 -0700 (PDT)
From: Vince <vince.pileggi at ontario.ca>
To: r-help at r-project.org
Subject: [R] deSolve - Function daspk on DAE system - Error
Message-ID: <1317525610060-3864298.post at n4.nabble.com>
Content-Type: text/plain; charset=us-ascii
I'm getting this error on the attached code and breaking my head but
can't
fi...
2011 Oct 02
0
deSolve - Function daspk on DAE system - Error
I'm getting this error on the attached code and breaking my head but can't
figure it out. Any help is much appreciated. Thanks, Vince
CODE:
library(deSolve)
Res_DAE=function(t, y, dy, pars) {
with(as.list(c(y, dy, pars)), {
res1 = -dS -dES-k2*ES
res2 = -dP + k2*ES
eq1 = Eo-E -ES
eq2 = So-S -ES -P
return(list(c(res1, res2, eq1, eq2)))
})
}
pars <- c(Eo=0.02,