Chapter 10

Introduction to LSODE

Note: This chapter will be included or not depending on whether LSODE is selected. Within the chapter, the subsections using IDL and MATLAB for graphical display are optional; none, either, or both can be selected.

10.1	Radioactive Decay
	10.1.1 Defining the Equations
	10.1.2 Specifying Initial Conditions
	10.1.3 Controlling Output Values
	10.1.4 Specifying Remaining Parameters
	10.1.5 Working Arrays
	10.1.6 Invoking LSODE
	10.1.7 Double Precision
	10.1.8 The Final Program
	10.1.9 Compiling and Running the Program
	10.1.10 Displaying the Solution Graphically with IDL
	10.1.11 Displaying the Solution Graphically with MATLAB
10.2	The Lorenz Attractor
	10.2.1 The Basic Equations
	10.2.2 The Solution with LSODE
	10.2.3 Displaying the Solution Graphically with IDL
	10.2.4 Displaying the Solution Graphically with MATLAB
10.3	Diffusion in One Dimension
	10.3.1 The Basic Equations
	10.3.2 Discretizing the Basic Equations
	10.3.3 The Solution with LSODE
	10.3.4 Displaying the Solution Graphically with IDL
	10.3.5 Displaying the Solution Graphically with MATLAB
10.4	Diffusion in Two Dimensions
	10.4.1 The Basic Equations
	10.4.2 Discretizing the Basic Equations
	10.4.3 The Solution with LSODE
	10.4.4 Displaying the Solution Graphically with IDL
	10.4.5 Displaying the Solution Graphically with MATLAB
10.5	Waves in One Dimension
	10.5.1 The Basic Equations
	10.5.2 Discretizing the Basic Equations
	10.5.3 The Solution with LSODE
	10.5.4 Displaying the Solution Graphically with IDL
	10.5.5 Displaying the Solution Graphically with MATLAB
10.6	Waves in Two Dimensions
	10.6.1 The Basic Equations
	10.6.2 Discretizing the Basic Equations
	10.6.3 The Solution with LSODE
	10.6.4 Displaying the Solution Graphically with IDL
	10.6.5 Displaying the Solution Graphically with MATLAB
10.7	Two-Dimensional Diffusion in Irregular Regions
	10.7.1 The Basic Equations
	10.7.2 Discretizing the Basic Equations
	10.7.3 The Solution with LSODE
	10.7.4 Displaying the Solution Graphically with IDL
	10.7.5 Displaying the Solution Graphically with MATLAB
	10.7.6 A Refined Solution
10.8	References
10.9	Exercises
10.A	Listing of radio.f
10.B	Listing of lorenz.f
10.C	Listing of diffus1d.f
10.D	Listing of diffus2d.f
10.E	Listing of wave1d.f
10.F	Listing of wave2d.f
10.G	Listing of house.f