IGEM:IMPERIAL/2009/M1/Modelling/Matlabcode1