clear
t=0:10:60;
T=[30 51 68 79 85.5 89 90];
tvork=0:1e-3:60;

%lähendamine
%1) kuupsplain
Tsplain=interp1(t,T,tvork,'spline');
%2) vähimruutude mõttes kuupfunktsioon
kordajad=polyfit(t,T,3);
Tkuup=polyval(kordajad,tvork);

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

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