Mathematica 9 is now available
Student Support Forum
Student Support Forum: 'DSolve and NDSolve' topicStudent Support Forum > General > Archives > "DSolve and NDSolve"

< Previous CommentHelp | Reply To Comment | Reply To Topic
Author Comment/Response
Bill Simpson
02/20/13 11:06pm

If you evaluate

With[{s = 0.1, h = 0.5, r0 = 0.01, r1 = 1},
odes2 = {...YourOdes2CodeUnChanged...}

and inspect the output very very very carefully, using Zoom may help, you will see of examples of

0. YT0A[t]

and similar constructs.

That is saying this function is being multiplied by "approximately zero" or floating point zero.

Those "worry me."

If you then kill your kernel and restart evaluation of

With[{s = 1/10, h = 1/2, r0 = 1/100, r1 = 1},
odes2 = {...YourOdes2CodeWithALLDecimalsReplacedByExactFractions...}

then the result is a little smaller, all decimal points have disappeared and the functions being multiplied by "approximately zero" are gone.

If you can duplicate that then try this next.

In[3]:= solnApproximate=NDSolve[odes2,{YT0A[t],YT0a[t],YT1A[t],MTA[t]},{t,0,100}]

Out[3]= {{

In[4]:= Plot[YT0A[t]/.solnApproximate,{t,0,100}]

Out[4]= ...NiceIncreasingPlotSnipped...

In[5]:= Plot[YT0a[t]/.solnApproximate,{t,0,100}]

Out[5]= ...NiceDecreasingPlotSnipped...

In[6]:= Plot[YT1A[t]/.solnApproximate,{t,0,100}]

Out[6]= ...NiceIncreasingPlotSnipped...

In[7]:= Plot[MTA[t]/.solnApproximate,{t,0,100}]

Out[7]= ..NiceIncreasingPlotSnipped...

This worries me a little less now.

If you can duplicate all this and verify that the plots make sense for your problem then perhaps we can think what the next step might be.

There is a great line in a 25 year old book by Donald Norman, called I believe "The Psychology of Everyday Things", something like "If you had 1/10 the trouble with a toaster that you have with a computer you would kick it down the stairs..." and now toasters have computers in them. I apologize for your frustration, it isn't even my fault, but I still apologize.

URL: ,

Subject (listing for 'DSolve and NDSolve')
Author Date Posted
DSolve and NDSolve Brian 02/20/13 3:15pm
Re: DSolve and NDSolve Bill Simpson 02/20/13 11:06pm
< Previous CommentHelp | Reply To Comment | Reply To Topic