function y1 = linreg(x, y, x1)
% Lineare Regression mit Abbildung

pf = polyfit(x, y, 1);  % y = pf(1)*x + pf(2)

yf = polyval(pf, x);
y1 = polyval(pf, x1);

figure
plot(x, y, 'ro', x1, y1, 'go')
hold on

plot(x, yf, 'b-')
hold off

end
