Mathematica 9 is now available
Student Support Forum
-----
Student Support Forum: 'Varying Parameter' topicStudent Support Forum > General > "Varying Parameter"

< Previous CommentHelp | Reply To Comment | Reply To Topic
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: ,

Subject (listing for 'Varying Parameter')
Author Date Posted
Varying Parameter Jane Shevtsov 09/25/05 00:48am
Re: Varying Parameter Daniele 09/28/05 10:36am
Re: Varying Parameter Jane Shevtsov 09/28/05 7:33pm
< Previous CommentHelp | Reply To Comment | Reply To Topic