Groups | Search | Server Info | Keyboard shortcuts | Login | Register [http] [https] [nntp] [nntps]

Groups > comp.soft-sys.math.mathematica > #1483

Re: Integation Anomoly

Show all headers | View raw

Sorry, I forgot to include the definitions of wCarStepR and wCarStepL
(thanks to Daniel).  With those definitions included in the proper
place, all the code should be as follows (you don't really need the
definition of wCarStepL to see the behavior I'm seeing, but I left it
in there).  The 'final' definition of wCarStep would be a piecewise
using 0 for t < 0, wCarStepR[t] from 0 <= t <= twTr, and wCarStepL[t]
from t > twTr.


row:=Sqrt[dwCar^2+1]*twTr/dwCar;
tow:=twTr+twTr/(dwCar^2)-row/dwCar;
wCarStepR[t_]=row-Sqrt[row^2-t^2];
wCarStepL[t_]=dwCar*(t-tow);
Assuming[twTr>0 && dwCar>0,(wCarStepR[twTr]==wCarStepL[twTr])//
FullSimplify]
thCA1s[t_]=Assuming[0<=t<=twTr && twTr>0 &&
dwCar>0,Integrate[wCarStepR[tau],{tau,0,t}]]
thCA2s[t_]=Assuming[0<t<=twTr && twTr>0 &&
dwCar>0,Integrate[wCarStepR[tau],{tau,0,t}]]
thCA1s[t]==thCA2s[t]
thCA1sd[t_]:=Assuming[0<=t<=twTr && twTr>0 &&
dwCar>0,Integrate[wCarStepR[tau],{tau,0,t}]]
thCA2sd[t_]:=Assuming[0<t<=twTr && twTr>0 &&
dwCar>0,Integrate[wCarStepR[tau],{tau,0,t}]]
thCA1sd[t]==thCA2sd[t]
Number Trials
dwCar=16;twTr=1/500;N[{row,tow},10]
thCA1s[t]
thCA1sd[t]
(thCA1s[t]==thCA1sd[t])//FullSimplify
(thCA2s[t]==thCA2sd[t])//FullSimplify
{thCA1sd[t]==thCA2sd[t],thCA1s[t]==thCA2s[t]}
Clear[dwCar,twTr]
dwCar=15;twTr=1/500;N[{row,tow},10]
thCA1s[t]
thCA1sd[t]
(thCA1s[t]==thCA1sd[t])//FullSimplify
(thCA2s[t]==thCA2sd[t])//FullSimplify
{thCA1sd[t]==thCA2sd[t],thCA1s[t]==thCA2s[t]}
(thCA1sd[t]==thCA2sd[t])//FullSimplify
a=thCA1sd[t]//FullSimplify;b=thCA2sd[t]//FullSimplify;a==b

-Brian L.

Back to comp.soft-sys.math.mathematica | Previous | Next — Previous in thread | Find similar | Unroll thread

Thread

Integation Anomoly blamm64 <blamm64@charter.net> - 2011-04-04 10:30 +0000
  Re: Integation Anomoly blamm64 <blamm64@charter.net> - 2011-04-05 10:43 +0000

csiph-web