clear,clc
global kf1 kf2 kr1 kr2
kf1=1;kf2=2;kr1=2;kr2=0.1;
C0=[1 0 0]';
[t,c]=ode45('kaksistep',[0 10],C0);
figure(1),plot(t,c,'linewidth',2);
set(gca,'linewidth',2,'fontsize',14)
xlabel('T = k_{f1}·t'),ylabel('C_i / [R]_0')