A set of demo programs for illustrating the use and behavior of the Matlab odexx integrators.

Codes and data for two problems

Radioactive Decay Chains:
	This problem solves for the activities of radioactive decay chains where nuclide 1 decays to nuclide 2 etc.
	
	to run
	load decaydata; % loads half-lives, initial conditions, integration times
	options=odeset; % sets defaults ode options (see odeset)
	figure; plotDecayChain(tHalf,a0,tRange,options,@ode45);

	to run stiff problem use decadatastiff;
	to try other ode solvers change @ode45 to @ode15s etc.
	
	Uses: decayChain.m to calculate RHS of the linear system da/dt = Aa

VanDerPol oscillator:
	
	This problems solves for the time dependent behavior of the vanderpol non-linear oscillator
	Equations and RHS are given in vanDerPol.m

	to visualize a run for a given value of mu for an integration time of 30, using ode45 with default options
	mu=-6
	tRange = [ 0 30];
	options=odeset;
	figure; plotVanDerPol(mu, tRange, options, @ode45)

	seeVDPtol runs this twice with different relative and absolute tolerances