%
% Supplementary Matlab files for the paper:
%
% " Deterministic Modelling and Stochastic Simulation of Pathways using Matlab"
% by M.Ullah, H.Schmidt, K-H.Cho and O.Wolkenhauer
%
% This script file simulates the MAPK pathway and generates the temporal
% evolution of three important concentrations
%
% We would appreciate a citation if these files are used by others.

clc;clear all;close all
load mapk       
[t, x] = ode15s(gma(D,k), [0 50], x0); 


[ax, h(1), h(2:3)] = plotyy(t, 1e3*x(:,4), t, x(:,[13 21]));

set(gcf,'PaperPosition',[2 2 8 6.3],'PaperType','A4');
 
set(ax, 'FontName','cmr','FontSize',8,'Position',[0.1264 0.15 0.7879 0.7643]);

set(h,{'LineWidth'},{1.5;1.5;1.2},{'LineStyle'},{'--';':';'-'},{'Color'},{'b';'k';'r'});

xlabel(ax(1),'time  (sec)');
ylabel(ax(1),'[MAPKKK-P]  (nM)');
ylabel(ax(2),'concentration  (nM)');
legend(h,{' MAPKKK-P',' MAPKK-P',' MAPK-P'},'FontName','Times','FontSize',8);
legend boxoff

% print  -depsc 'figures\timehistory.eps'