%Start
clear all
addpath '.\Libary'
load data_conc
load data_conc_big
cold = c;

figure(1);
hold on
plot(t,c,'k-');

%12, 10%, yes
tin = t12;
cin = c12;
[spa1,spb1,spc1] = spline_no_2(cold,tin,cin,0.10,0,[0 1 2]);

te = [0:0.02:2]';
drawspline(spa1,te,'k--');
drawspline(spb1,te,'k--');
drawspline(spc1,te,'k--');

%
tin = t12;
cin = c12;
[spa1,spb1,spc1] = spline_h_2(cold,tin,cin,0.10,0,[0 1 2]);

te = [0:0.02:2]';
drawspline(spa1,te,'k:');
drawspline(spb1,te,'k:');
drawspline(spc1,te,'k:');

axis([0 2 0 1]);