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