| Author |
Comment/Response |
Jane Shevtsov
|
 09/28/05 7:33pm
When I try to use Piecewise, I get an error saying, "LogicalExpand::elist: List encountered during logical expansion of {#6 == 0 \
&& Solve`Auxiliary[#11 > 0]}."
Here's the actual code. Turn off spell check before trying it!
varsBounce={rBounce[t], eBounce[t],lBounce[t],vBounce[t]};
gBounce=Γ τ -μ rBounce[t]-β rBounce[t] vBounce[t];
hBounce=(1-ρ) β rBounce[t] vBounce[t]+α[t]*lBounce[t]-δ eBounce[t];
jBounce=ρ β rBounce[t] vBounce[t]-α[t]*lBounce[t]-μ lBounce[t];
kBounce=p eBounce[t]-σ vBounce[t];
α[t_]:=Piecewise[{
0.36,{vBounce'[t]\[Equal]0&&vBounce''[t]>0}},{0.036,!{vBounce'[t]\
\[Equal]0&&vBounce''[t]>0}}]
eqnsBounce={rBounce'[
t]\[Equal]gBounce,eBounce'[t]\[Equal]hBounce,lBounce'[t]\[Equal]
jBounce, vBounce'[t]\[Equal]kBounce};
initsBounce={rBounce[2000]\[Equal]24.5338,
eBounce[2000]\[Equal]0.720501, lBounce[
2000]\[Equal]0.638742, vBounce[2000]\[Equal]36.0251,vBounce'[2000]\
\[Equal]0};
Γ=1.36;μ=1.36*10^-3; τ=0.2; ρ=0.1; σ=2; δ=0.33; p=100;
β=0.0001;
solBounce2000=varsBounce2000/.NDSolve[Join[eqnsBounce,initsBounce],\
varsBounce,{t,2000, 4000}][[1]];
URL: , |
|
