求解方法：

1、用三次样条函数interp1（）插值xi=1:0.1:15;yi=interp1(x,y,xi,'spline');