clear
t=0:10:60;
T=[30 51 68 79 85.5 89 90];

tvork=0:1e-2:60;

%lähendamine
Tsplain=interp1(t,T,tvork,'spline');

kordajad=polyfit(t,T,3);
Tkuup=polyval(kordajad,tvork);

%joonis
plot(t,T,'b*',tvork,Tsplain,tvork,Tkuup)
legend('algandmed','kuupsplain','kuupfunktsioon vähimruutude mõttes')

%arvutamine
vastus1=interp1(t,T,54,'spline')
vastus2=polyval(kordajad,54)