clear
global K1 K2
K1=10;K2=0.1;
y0=[1 0 0 0]';
[t,y]=ode45('perakkain',[0 10],y0);
figure(1),plot(t,y,'linewidth',2)
set(gca,'linewidth',2,'fontsize',14)
xlabel('k_A·t'),ylabel('[X]/[A]_0')