cancel
Showing results for 
Search instead for 
Did you mean: 
cancel
Showing results for 
Search instead for 
Did you mean: 

ODE solve tolerance

SOLVED
Highlighted
Regular Member

ODE solve tolerance

Greetings!

I am solving a pair of odes using mathcad built in odesolve function. I want to know if there is a way to alter the default tolerance and number of iterations of this function.

Because when i solve it its gives a statement that solution is not converging to a solution refer to Help.

 

Thanks.

1 ACCEPTED SOLUTION

Accepted Solutions
Highlighted

Re: ODE solve tolerance

You can change the tolerance via the menu or by typing TOL:=... at the top.

But the problem is that the numerical derivative of v is not defined at the end values (0 and 0.34) - its wrongly set to 0!.

So you may extend the ends a little bit like shown in the picture to make it work.

B.JPG

 

View solution in original post

22 REPLIES 22
Highlighted

Re: ODE solve tolerance

You can change the value of TOL.  Set TOL:=10^-6 say (the default value is 10^-3). 

 

However, this might not be the root cause of your problem.  There are many reasons why the solution might not converge.  It would be a good idea to upload your worksheet so the community can investigate.

 

Alan

Highlighted

Re: ODE solve tolerance


@Asad wrote:

Greetings!

I am solving a pair of odes using mathcad built in odesolve function. I want to know if there is a way to alter the number of iterations of this function.

 

 

Thanks.


May by so

Odesolve(y, t, [number of iterations])

Highlighted

Re: ODE solve tolerance

I am using mathcad prime 2.0 to solve the ODEs. In the solve block I have typed x:=odesolve(x(z),0.47,1000000000) and now it gives a new error that there is not enough memory for this calculation. I am uploading the work sheet, help would be highly appreciated! A little about the ODEs: there are two varibles velocity (v) and concentration (x) that are varying along the z-axial direction. Taking an initial guess for dx/dz in the first equation to obtain and answer for v and dv/dz. Then using the answers obtained in last step I put them in the second equation to solve for x and dx/dz. After successfully getting an answer I will then make a mathcad while loop with these solve blocks and iterate them until a desired solution is obtained. Thanks
Highlighted

Re: ODE solve tolerance

and if i put another zero

x:=odesolve(x,0.47, 10000000000)

 

it says found a number with a magniture greater than 10^307 while trying to evaluate this expression

Highlighted

Re: ODE solve tolerance

Is there a way to change the tolerance?

Highlighted

Re: ODE solve tolerance

You can change the tolerance via the menu or by typing TOL:=... at the top.

But the problem is that the numerical derivative of v is not defined at the end values (0 and 0.34) - its wrongly set to 0!.

So you may extend the ends a little bit like shown in the picture to make it work.

B.JPG

 

View solution in original post

Highlighted

Re: ODE solve tolerance

Or, just change the end time for the v ode to something significantly larger:

 

Alan

 

odes.PNG

Highlighted

Re: ODE solve tolerance

Hmm, thats interesting, that the end value of the first odesolve block has such an impact on the acceptance  of the start value in the second solve block.

Its also a matter of the precision settings. With TOL:=10^-9 the solve block still fails with end value 0.5 in the first one.

Highlighted

Re: ODE solve tolerance

I've just noticed the following taken from the Prime 2 help file:

 

tol.png

 

Alan

Highlighted

Re: ODE solve tolerance

After your first post I played around and found that your way of working around the problem failed as long as I had TOL:=10^-9 in front. When I changed it to the default 10^-3 the second odesolve block worked OK. So in contrary to what the help says it look like the setting of TOL has some effect on the solve block.

 

When I open the sheet you posted (in Prime 5) the second odesolveblock fails - even when I delete the last line and retype it.

 

Strange effects, indeed.

Highlighted

Re: ODE solve tolerance

A similar problem occurs in mathcad 15.  However, if Adams/BDF is replaced by Radau in the second odesolve it works ok with 0.34 as the end time in both.  Of course, this isn't an option in Prime.

 

Alan

Highlighted

Re: ODE solve tolerance

Can both of you please send me your worksheet file because my answer is still not converging.

Highlighted

Re: ODE solve tolerance

Strangely, when I reloaded my worksheet it also failed!  However, by deleting the last odesolve command and simply re-entering it the calculation converged!

 

Worksheet attached.

 

Alan

 

Highlighted

Re: ODE solve tolerance

Hi Asad,

Have you tried to change algorithm?

Ode options.jpg

Highlighted

Re: ODE solve tolerance


@-MFra- wrote:

Hi Asad,

Have you tried to change algorithm?

 


 Only in Mathcad 15!

Highlighted

Re: ODE solve tolerance

RIght! So a combination of both of your solutions worked for me and if the solution does not converge I just alter the starting and ending points of the boundary conditions.

 

I worked on it and:

v:=odesolve(v(z),2,100000000) with v(0)=0.4

x:=odesolve(x(z),1.2,1000000) with x(10^-307)=x_in

 

then for the second iteration

v:=odesolve(v(z),0.9) with v(10^-306)=0.4

x:=odesolve(x(z),0.47) with x(10^-305)=x_in

 

Thanks a lot for the help!

Highlighted

Re: ODE solve tolerance

Also, is there any other function besides odesolve to solve the differential equations like the ones stated in comments such as radau or any other. For better accuracy..

Highlighted

Re: ODE solve tolerance

As was already written in contrary to Mathcad 15 and below we can not manually chose the algorithm used by odesolve - we have to rely on the automatic selection.

For other ways to solve ODEs numerically you may look up the various function in the help.

See the attached picture.

Highlighted

Re: ODE solve tolerance

Can you send me your worksheet file in which you have set the initial and final boundary conditions for the derivative of v.

The worksheet of the picture you posted.

 

Thanks!

Highlighted

Re: ODE solve tolerance


@Asad wrote:

Can you send me your worksheet file in which you have set the initial and final boundary conditions for the derivative of v.

The worksheet of the picture you posted.

 

Thanks!


I did not save and store that worksheet.

I could change the few values again and post but it won't help as I am using a newer version of Prime and you would not be able to open the worksheet as you wrote that you are using Prime 2. Newer versions of prime are not able to save a sheet in an older format.

Highlighted

Re: ODE solve tolerance

The other functions are still there and usable, but they require setting up and they don't use a solve block.  They also don't improve the accuracy.  (If you have accuracy issues the old wisdom suggests rescaling your problem.)

 

Spoiler
Spoiler
Radau, Rkfixed, Rkadapt and the others do not return a function like odesolve; they return a matrix of points with the ordinate points in the first column and solutions in the subsequent  columns.
Highlighted

Re: ODE solve tolerance


@FredKohlhepp wrote:

The other functions are still there and usable, but they require setting up and they don't use a solve block.  They also don't improve the accuracy.  (If you have accuracy issues the old wisdom suggests rescaling your problem.)

 

Spoiler
Spoiler
Radau, Rkfixed, Rkadapt and the others do not return a function like odesolve; they return a matrix of points with the ordinate points in the first column and solutions in the subsequent  columns.

When I open this (the "Quote function" did not replicate it) I saw the first paragraph and 

>Spoiler

 

Two clicks later I saw this:

Capture.PNG

So, what's a "Spoiler"?

Announcements