x=[0.01  0.05   0.25  1/3   1   ];
f=[100    20       4    3    1  ];
Da=0;Db=-1;
kubicky_spline(x,f,Da,Db);
