x=[10,15,20];
y=[227.04,362.78,517.35];
t=16;
n=length(x);
b0=y(1);
b1=(y(2)-y(1))/(x(2)-x(1));
b2=(y(3)-y(2))/(x(3)-x(2));
result=b0+b1*(t-x(1))+b2*(t-x(2))*(t-x(1));
disp(result)