Re: Using "If" and "NDSolve" together
- To: mathgroup at smc.vnet.net
- Subject: [mg106859] Re: Using "If" and "NDSolve" together
- From: RBenf <rui.benfeitas at gmail.com>
- Date: Mon, 25 Jan 2010 05:07:15 -0500 (EST)
Hello Bobby, thanks for all the help! It was good to see how you can simplify what I wrote into just a few coding lines. I'm still learning how to work with Mathematica, and from what I've seen, it will still be a long way until I can use it effectively.
Although the code that you gave me does not do exactly what I want to do, I will surely use it in order to solve my problem.
Thanks a lot!
> Oh, what a tangled web you weave.
>
> Clear[funcao]
> funcao[first_, last_] :=
> Module[{k0 = 10^-2, k1 = 10^-5, k2 = 10^-6, k3 =
> = 10, a, b, c},
> a = a /. First@
> NDSolve[{a'[t] == k0 Boole[first <= t <= last]
> last] - k1 a[t],
> b'[t] == k1 a[t] - k2 b[t], c'[t] == k2 b[t]
> k2 b[t] - k3 c[t],
> a[0] == 2 10^-6, b[0] == 0, c[0] == 0}, {a,
> 0}, {a, b, c}, {t, 0,
> 1000}];
> Plot[a@t, {t, 0, 1000}, PlotRange -> All]
> ]
> funcao[0, 100]
> funcao[0, 900]
> funcao[100, 900]
>
> Bobby