x=[500 1000 2000 4000];
y=[3.08 1.56 0.78 0.39];
xi=500:4000;
yi=interp1(x,y,xi,'cubic');
plot(xi,yi,':',x,y,'o');
set(gcf,'color','w');
