ee.m

可以方便的调用这个函数

year=[0.125 0.250 0.375 0.500 0.625 0.750]';
price=[0.79618 0.77334 0.74371 0.70413 0.65632 0.60228]';
y=linspace(min(year),max(year),10);
p=newtint(year,price,y);

f(0.1581)=0.79618+(-0.18272)*(0.1581-0.125)+(-0.21728)*(0.1581-0.125)*(0.1581-0.250)+(-0.26965)*(0.1581-0.125)*(0.1581-0.250)*(0.1581-0.375)

f(0.6367)=0.74371+(-0.31664)*(0.6367-0.375)+(-0.26336)*(0.6367-0.375)*(0.6367-0.500)+(0.17066)*(0.6367-0.375)*(0.6367-0.500)*(0.6367-0.625)