clear
%algandmed
t=0:10:60; %aeg 10minutilise sammuga
T=[30 51 68 79 85.5 89 90]; %temperatuur 

%võrk
tvork=0:1e-3:60;

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

%2) kuupfunktsioon vähimruutude mõttes
kordajad=polyfit(t,T,3); %otsime kuupfunktsiooni, sestap 3
Tkuup=polyval(kordajad,tvork);

%joonis
plot(t,T,'r*',tvork,Tsplain,'k',tvork,Tkuup)
legend('algandmed','kuupsplain','kuupfunktsioon')
xlabel('aeg')
ylabel('temperatuur')
title('Töötava mootori mähise temperatuur ajast sõltuvana')
grid on

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