function prolozeni
t=[1:5];
pr=[20,24,30,31,37];
plot(t,pr,'ro')
xj=[0.5:0.1:5.5];
k=polyfit(t,pr,3);
disp('koeficienty polynomu')
disp(k)
y3=polyval(k,xj);
hold on
plot(xj,y3,'k')
hold off
end

