How to set the NDSolve method to LSODA

Question

I notice that off all the Method options available for NDSolve[...], LSODA is invoked quite automatically for stiff equations when Method -> Automatic. Is there a way to set the Method -> LSODA? I ask because as LSODA switches between Adams' and Gear's method, it should have sub-options to it that I could set (like order etc.)?
Under tutorial/NDSolveStiffnessSwitching I encountered Needs["DifferentialEquations`NDSolveProblems"]. What other example problems are available and where can I find documentation on them?

Edit (03/24/2011):

In light of comments made to this post, I'd also refer readers to page 409 of "The Mathematica Guidebook for Numerics". Springer has this book online and can be found here.

score 12 · Accepted Answer · edited Mar 24 '12 at 12:03

12

If you load the package using Needs as shown in the Tutorial and then evaluate $NDSolveProblems you will list a list of names of all the problems included in the package.

I think I can now also answer the first part of the question. The method has to be called as

NDSolve`LSODA

rather than LSODA (for some reason it is in the NDSolve` context, which suggests that it is not intended for "ordinary users"). To see its options, just evaluate

?NDSolve`LSODA

The whole issue of these "hidden contexts" in Mathematica could be the topic of an interesting discussion ;-)

edited Mar 24 '12 at 12:03

Szabolcs

234,956
30
623
1,263

answered Mar 22 '12 at 21:06

Andrzej Kozlowski

2,839
19
20

?NDSolveLSODA gives me: Options[NDSolveLSODA]={LinearSolveMethod->Automatic,MaxDifferenceOrder->12} This is rather useful. However, what does NDSolve`LSODA do? Did it load the LSODA method for use? – dearN Mar 24 '12 at 15:49
It is itself the name of a method for NDSolve, that is, you use it by evaluating NDSolve[eqn,Method->NDSolve`LSODA]. Whether it will do any good depends of course on eqn. – Andrzej Kozlowski Mar 24 '12 at 16:16
Perhaps it's worth mentioning that the symbol NDSolve`LSODA does not get loaded (and thus cannot be queried) until NDSolve[] is used at least once. – Szabolcs Mar 24 '12 at 16:17
A bit of spelunking shows that NDSolve`LSODA is defined like this: http://pastebin.com/0nayVBHP This confirm what @Andrzej said that it is just a method identifier. You can learn about how methods work here (I didn't have time to read this, but it seems like a relevant resource). – Szabolcs Mar 24 '12 at 16:23
@Szabolcs According to "The mathematica guidebook of numerics" (whose link I've included in the post now) says that LSODA can be loaded.. I am trying to figure out how and what the options are and how to set them. (Thanks to me being a beginner!) – dearN Mar 24 '12 at 16:24
1

More precisely: it can be used as NDSolve`LSODA even the first time but once NDSolve is used, you can call it just LSODA. – Andrzej Kozlowski Mar 24 '12 at 16:25
@AndrzejKozlowski Thanks for the input. – dearN Mar 24 '12 at 16:31

acl · Answer 2 · 2012-03-23T02:01:22.067

5

Does NDSolve[eqn, opts, Method->LSODA] not work? It should, according to this.

The problems seem to be at

FileNameJoin[{$InstallationDirectory, "AddOns", "ExtraPackages", 
  "DifferentialEquations", "NDSolveProblems.m"}]

edited Mar 23 '12 at 02:01

answered Mar 22 '12 at 20:58

acl

19,834
3
66
91

I appreciate the link to the pdf. Thanks! – dearN Mar 22 '12 at 22:47
@acl the linked file does not exist anymore... could you update it, please? – Nobody Nov 29 '21 at 05:54

How to set the NDSolve method to LSODA

2 Answers2

Linked