clc; close all; clear all;
t=[0, 4, 8,  12, 16, 20, 24]; %WRES HMERAS POY GINONTAI METRHSEIS 8ERMOKRASIAS
T=[5, 7, 10, 13, 15, 13,  6]; %ESTW OI METRHSEIS 8ERMOKRASIWN AIS8HTHRA;
 
Tmean=mean(T);                %MESH 8ERMOKRASIA
Tstd=std(T);                  %TYPIKH APOKLISH 8ERMOKRASIWN
Tmax=max(T);                  %MEGISTH 8ERMOKRASIA
Tmin=min(T);                  %ELAXISTH 8ERMOKRASIA
Tmedian=median(T);            %MESEA TIMH 8ERMOKRASIAS
Trange=range(T);              %EYROS DIAKYMANSHS 8ERMOKRASIWN
Tvar=var(T);                  %DIASPORA 8ERMOKRASIWN
Tstd=sqrt(var(T));            %Tstd WS ANWTERW
fprintf('Tmean=%0.2f\nTstd=%0.2f\nTmax=%0.2f\nTmin=%0.2f\n',Tmean,Tstd,Tmax,Tmean);
fprintf('Tmedian=%0.2f\nTrange=%0.2f\nTvar=%0.2f\nTstd=%0.2f\n',Tmedian,Trange,Tvar,Tstd);

%SXEDIASH DIAKYMANSHS 8ERMOKRASIWN SE BOXPLOT
figure(1);
subplot(2,1,1);
boxplot(T);
grid on; 
ylabel('хеялойяасиес (╟C)');
title('BOXPLOT диайуламсгс хеялойяасиым');  

%SXEDIASH E3ISWSHS POY NA ANTIPROSWPEYEI TA DEDOMENA
subplot(2,1,2);
[y] = fit(t', T', 'poly3'); %SYNARTHSH MATLAB GIA EYRESH E3ISWSHS POY NA ANTIPROSWPEYEI TA DEDOMENA  
fprintf('\nсумтекестес 3-бахлиоу покуымулоу:\n');
fprintf('p1=%7.5f\np2=%7.5f\np3=%7.5f\np4=%7.5f\n',y.p1,y.p2,y.p3,y.p4)

%YPOLOGISMOS TIMWN 8ERMOKTRASIWN ME BASH THN E3ISWSH
for i=1:length(T)
    z(i)=y.p1*t(i)^3+y.p2*t(i)^2+y.p3*t(i)+y.p4;
end

%SXEDIASH SHMEIWN KAI KAMPYLHS
plot(t,T,'*b','markersize',14,'linewidth',2);
hold on;
plot(t,z,'--k');
grid on;
legend('хеялойяасиес','йалпукг летабокгс хеялойяасиым',2);
xlabel('ыяес');
ylabel('хеялойяасиес (╟C )');
title('покуымулийг енисысг 3ОУ BAхMOY циа летабокг хеялойяасиас')
