Solving a radioactive decay system using different numerical methods in MatLab

This is a project I wrote during the course of Numerical Methods in Earth Sciences  during my Erasmus in Budapest.

You may find it useful.

download

Here is the content of the read.txt file:

main script:
- radioactive_decay_system_script.m
  this is the main script, that uses all the following functions
  [ JUST RUN IT! ]

functions:
- eulero_method_system.m
  function to implement the simple eulero method
- eulero_method_modif_system.m
  function to implement the modified eulero method
  (better then the simple eulero)
- runge_kutta_method_system.m
  function to implement the 4-order runge-kutta method
- adams_bashforth_method_system.m
  function to implement the adams-bashforth method for s=2,3,4,5
- adams_bashforth_method_system_gener.m
  generalized (s>=2) function to implement the adams-bashforth method
- coefficients.m
  function used to calculate the coefficients for the adams-bashforth method
  this function is utilized by the function adams_bashforth_method_system_gener.m 
- rad_decay_system.m
  function for the system of ode we want to solve

Different methods provides slighty different precisions.

Enjoy it and ask questions if you wish!