ut <- daspk(y = state, dy = dval, times = z, res = reaction, parms = 0)
head(out)
plot(out)
-Message d'origine-
De : Berend Hasselman [mailto:b...@xs4all.nl]
Envoyé : jeudi 11 juillet 2013 14:13
À : Raphaëlle Carraud
Cc : r-help@r-project.org
Objet : Re: [R] Differential problem
On 11-0
ms = 0))
Is it wrong? When I modify the order, I obtain again that object dA is not
found, so I guessed the doc was right.
-----Message d'origine-----
De : Berend Hasselman [mailto:b...@xs4all.nl]
Envoyé : jeudi 11 juillet 2013 12:33
À : Raphaëlle Carraud
Cc : r-help@r-project.org
Objet : Re:
On 11-07-2013, at 13:53, Raphaëlle Carraud
wrote:
> Sorry for the bug, I had eliminated some lines to avoid making the program
> too big. Here is the version that works :
>
> reaction<-function(z, state, dval, parameters) {
> with(as.list(c(state)),{
># rate of change
>
>Tr <- 273+9
deSolve)
out <- daspk(y = state, dy = dval, times = z, res = reaction, parms = 0)
head(out)
plot(out)
-Message d'origine-
De : Berend Hasselman [mailto:b...@xs4all.nl]
Envoyé : jeudi 11 juillet 2013 11:18
À : Raphaëlle Carraud
Cc : r-help@r-project.org
Objet : Re: [R] Differential
nd shouldn't the first line of the reaction function be this
with(as.list(c(state,dval,parameters)),{
in stead of this
with(as.list(c(state)),{
The call of daspk also seems incorrect; shouldn't it be
out <- daspk(y = state, dy = dval, times = z, res = reaction, parms =
parameter
On 11-07-2013, at 09:13, Raphaëlle Carraud
wrote:
> Hello,
>
> I have the following differential equation system to solve, the variables I
> wish to obtain being A, B, C, D, E, I, G and H.
>
>0 = -dA + dB + 2*dC - 2*r1 - 2*r5
>0 = dA + dD + r1 + r4
>0 = K2 - C/B^2
>0 = K3 - D
6 matches
Mail list logo
