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

%võrk
tvork=0:0.001:60;

%lähendamine
%1) kuupsplain
Tsplain=interp1(t,T,tvork,'spline');
%arvutame kuupsplaini väärtused tvork igas punktis

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

%joonis
plot(t,T,'r*',tvork,Tsplain,tvork,Tkuup,'k')
legend('algandmed','kuupsplain','kuupfunktsioon vähimruutude mõttes')
title('Ülesanne 82')
xlabel('aeg')
ylabel('temperatuur')

%arvutame temperatuuri, kui mootor on töötanud 54 minutit
vastus1=interp1(t,T,54,'spline')
vastus2=polyval(kordajad,54)