Displaying 1 result from an estimated 1 matches for "lvcomp2".
Did you mean:
lcomp2
2001 May 11
1
lsoda
...cond, "a vector
(possibly with a `names' attribute) of global values that are
required at each point in `times'."
I didn't understand from the example why the massbalance = sum(y) was needed
for that system of equations. Why is it needed?
For my function below ('lvcomp2'), I don't know what global values lsoda
needs at each point in time (see???), other than the previous time step's
values of y and the parms. Any help?
Also, any rule of thumb advice regarding rtol and atol would be appreciated.
p <- c(r1=.01, r2=.01, a12=0.5, a21=0.6, k1=1,k2=1)
t...