Displaying 3 results from an estimated 3 matches for "yd3".
Did you mean:
d3
2004 Jun 10
0
lsoda with arbitrary zero thresholds (with psuedo-solution)
...tol <- c(1e-6, 1e-10, 1e-6)
> times <- seq(0,1000)
> lsexamp <- function(t, y, p)
> {
> if(y[1] < .4) yd1 <- -y[1] ### These if, else statements are new
> else yd1 <- -p["k1"] * y[1] + p["k2"] * y[2]*y[3]
> if(y[3] < .4) yd3 <- -y[3] ### These if,else statements are new
> else yd3 <- p["k3"] * y[2]^2
> list(c(yd1,-yd1-yd3,yd3),c(massbalance=sum(y)))
> }
> out <- lsoda(c(.5,0,.5),times,lsexamp, parms, rtol=1e-4, atol=
my.atol,
> hmax=.1)
> matplot(out[,1],out[,2:5], ty...
2004 Jun 14
1
olesolve: stepsize
...quot;kl"] + p["fx"]*y[4])*y[1]*p["vg"] +
p["k12"]*y[2])/p["vg"] - (p["k21"] + p["k01"])*y[1]
yd2 <- p["k21"]*y[1]*p["vg"] - (p["k12"] + p["k02"] +
(1-p["fx"])*y[4])*y[2]
yd3 <- (1-p["f"])*s(t)*p["kg"]*y[1]/p["vi"] - p["ki"]*y[3]
yd4 <- -p["p2"]*y[4] + p["p3"]*(y[3]-ib)
list(c(yd1, yd2, yd3, yd4))
}
__________________________________________________
First is that The computing time f...
2004 Dec 13
2
Switching to Mac, suggestions? (was switching to linux)
I'm looking to switch to Mac platform. Anyone had any experience
with that? I'm expecting on a power G4 laptop later this week.... hope
R behaves...
-- dok