clear all, close all, clc

T=-15:5:0;
P=[1.24 1.95 3.01 4.58];

P81=P(2)+(P(3)-P(2))/(T(3)-T(2))*(-8-T(2))

plot(T,P,'O',-8,P81,'*',T(2:3),P(2:3))
hold on

aktT=T(2:end)';
A=[aktT.^2 aktT ones(size(aktT))];
aktP=P(2:end)';
c=A\aktP;
kortT=aktT(1):0.01:aktT(end);
kortP=polyval(c,kortT);
P82=polyval(c,-8)
plot(kortT,kortP, -8,P82,'*')

aktT=T(1:end-1)';
A=[aktT.^2 aktT ones(size(aktT))];
aktP=P(1:end-1)';
c=A\aktP;
kortT=aktT(1):0.01:aktT(end);
kortP=polyval(c,kortT);
P83=polyval(c,-8)
plot(kortT,kortP, -8,P83,'*')
