conc.m

Incorporating Prior Knowledge in Cubic Spline Approximation - Application to the Identification of R

% CONC

global b1 b2 b3 b4
b1 = 1;
b2 = 0.5;
b3 = 10;
b4 = 5;
[t,c] = ode45('freact1',[0:0.02:2],[1 0 0]);
plot(t,c(:,1),'r',t,c(:,2),'g',t,c(:,3),'b');

%12-points
vec=round(linspace(1,1001,12));
t12 = t(vec);
c12 = c(vec,:);

%31-points
vec=round(linspace(1,1001,31));
t31 = t(vec);
c31 = c(vec,:);

save data_conc t12 c12 t31 c31
save data_conc_big t c