clear model data parama options

A0 = 0.02090; B0 = A0/3; C0 = 0; D0 = 0; E0 = 0;
data.y0 = [A0;B0;C0;D0;E0];

y0= data.y0;

k(1)=14.402;
k(2)=1.5663;
k(3)=0.29042;

tspan=(0:501/3000:501);
options = odeset('RelTol', 1e-12, 'AbsTol', 1e-12);
[t,z1]=ode45(@(t,y) himmelode(t,y,k),tspan,data.y0,options);

xi=1e-5;
S=[];
colorspec= {[0 0 0]; [0 0 1]; [1 0 1]};

linespec={'-', '-', '-'};

for i=1:3
   pad_k=k;
   pad_k(i)=k(i)*(1+xi);
   [t,z2]=ode45(@(t,y) himmelode(t,y,pad_k),tspan,data.y0,options);
   mz=mean(z1,1);
   B = repmat(mz,3001,1);
   sens=(z2-z1)/(xi*k(i))*k(i)./B;
   
   figure(1); hold on; plot(t,sens(:,1), 'Color', colorspec{i},...
       'LineStyle',linespec{i}); xlabel('Time [d]'); ylabel('Sens.A');         

   figure(2); hold on; plot(t,sens(:,2), 'Color', colorspec{i},...
       'LineStyle',linespec{i}); xlabel('Time [d]'); ylabel('Sens.B');
   figure(3); hold on; plot(t,sens(:,3), 'Color', colorspec{i},...
       'LineStyle',linespec{i}); xlabel('Time [d]'); ylabel('Sens.C');
   
   figure(4); hold on; plot(t,sens(:,4), 'Color', colorspec{i},...
       'LineStyle',linespec{i}); xlabel('Time [d]'); ylabel('Sens.D');
   
   figure(5); hold on; plot(t,sens(:,5), 'Color', colorspec{i},...
       'LineStyle',linespec{i}); xlabel('Time [d]'); ylabel('Sens.E');
   
    S=[S; sens];

end 



for i=1:5
    figure(i)
    legend('k1','k2','k3');
end