%DANE:
t = [0,5,10,15,20,25,30,35,40]; %czas
T = [90,77,64,56,50,46,43,40,37]; %temperatura
T0 = 22; %temperatura pokojowa

%OBLICZENIA:
W = T - T0; %różnica pomiędzy temp cieczy, a temp pokojową
M = log(W); %zlogarytmowanie różnic temperatur
p = polyfit(t,K,1); %aproksymacja
k = -p(1); %wyznaczenie współczynnika stygnięcia
y = polyval(p,t); %obliczenie wartości wielomianu

%WYKRES:
hold on
plot(t,M,'o'); %rysowanie punktów
plot(t,y); %rysowanie prostej
xlabel('czas[min]') %nazwanie osi OX
ylabel('ln(T-T0)') %nazwanie osi OY