function pinene();

x0=[100 0 0 0 0]; 
tspan = [0,1230,3060,4920,7800,10680,15030,22620,36420];
k=[0.59259,0.29634,0.20471,0.27445,0.39981];
[t,p] = ode15s(@ode_pinene,tspan,x0,[],k);
plot(t,p(:,1),'ro-');
hold on;
plot(t,p(:,2),'bs-');
plot(t,p(:,3),'gd-');
plot(t,p(:,4),'kx-');
plot(t,p(:,5),'mv-');
legend('X_1','X_2','X_3','X_4','X_5');



function dx = ode_pinene(t,x,k);
dx = zeros(5,1);   
dx(1)=-(k(1)+k(2))*x(1);
dx(2)=k(1)*x(1);
dx(3) = k(2)*x(1)-(k(3)+k(4))*x(3)+k(5)*x(5);
dx(4) = k(3)*x(3);
dx(5) = -k(4)*x(3)+k(5)*x(5);




