x=[0 1 2 3  4 5 6 7 8  9 10];
f=[3 5 0 2 -4 0 2 2 -2  2  3];
Da=1;Db=1;
kubicky_spline(x,f,Da,Db);
