function T8A2

x1=1;
y1=32;
x2=2;
y2=60;
x=0:3;%Länge der Geraden
x3=1.5;
y=28*x+4;%Geradengleichung
y3=28*x3+4%Geradengeleichung für x3
plot(x1,y1,'^',x2,y2,'^',x,y,x3,y3,'^')

%% Berechnung des quadratischen vertikalen Abstands von x3 zur Geraden
y4=44;%Messwert
r=(y4-y3)^2
figure
plot(x1,y1,'^',x2,y2,'^',x,y,x3,y4,'^')

%% Näherungsgerade
y5=28*x+3;
figure
plot(x1,y1,'^',x2,y2,'^',x,y5,x3,y4,'^')

%% Berechnung des mittleren quadratischen Fehlers
a=[32 60 44];%Messwerte
b=[32 60 46];%Schätzwerte
F=1/3*(a-b).^2;
sum(F)

%% Polynom
a1=[1 1.5 2];%Messwerte
b1=[32 44 60];%Schätzwerte
[p,S]=polyfit(a1,b1,1)
%plot(p)
[z,delta]=polyval(p,a1,S)
s=1:3;
figure
errorbar(s,z,a1,'r')
hold on
plot(z)








