Numerical Solution of Ordinary Differential Equations
Matlab Component

This page contains links to the hand-outs, problem sheets and designed for a Matlab practical on numerical methods for ordinary differential equations held in the 99/00 autumn semester at Brunel University.

Handouts

All handouts are downloadable as post-script files.
Organisational Information
Guide: "How to use Matlab"
Problem Sheet 1
Assignment 1
Problem Sheet 2
Assignment 2

Solutions to Problem Sheets/Assignments

Solutions to Problem Sheet 1


Euler_One_Step.m
Euler.m
test_Euler.m
better_Euler.m
test_better_Euler.m
 

Solutions to Assignment 1

Problem 1

Euler_One_Step2.m
f_example.m
Testing of Euler_One_Step2

Problem 2

new_Euler.m
f_ex2.m
M-file to test new_Euler

Problem 3

The required M-file (also includes instructions for Problem 4)

Problem 4

Runge3_One_Step.m
Runge_Kutta3.m
The required M-file
 

Solutions to Problem Sheet 2


Euler_Step_Sys.m
Euler_Sys.m
Test file for Problem 2: test_Euler_Sys.m
Implementation of function f : f_prob2.m
Implementation of function g : g_prob2.m
 

Solutions to Assignment 2

Problem 1

Runge3_Step_Sys.m, a function called by Runge_Kutta3_Sys
Runge_Kutta3_Sys.m
Implementation of function f : f_prob2.m
Implementation of function g : g_prob2.m
M-file to test Runge_Kutta3_Sys

Problem 2

Implementation of function f : f_a2p2.m
Implementation of function g : g_a2p2.m
M-file to compute and plot the maximum relative errors
The answer to the final question is, that the relative error is approximately divided by eight if the number of mesh points is doubled. This is the expected behaviour of a third order method.

Problem 3

Implementation of function f : f_a2p3.m
Implementation of function g : g_a2p3.m
M-file to compute the height of the spear
plot of the different flight-paths (post-script)
 

This page maintained by Tilo Arens

last changed: 22/3/2001
email me at: arens@numathics.com
visit my home page at: http://www.numathics.com/arens