function [a,b]=myLinearRegression(x,y)
n=length(x);
if length(y)~=n, error('dumm gelaufen'); end
x=x(:); y=y(:);
sx=sum(x); 
y=0:1:100;
xp=min(x):.01:max(x);
yp=a+b.*xp;
astr=num2str(a,'%6.4f');
bstr=num2str(b,'%+6.4f');
fname=['y(x)=' astr bstr 'x'];
plot(x,y,'.k',xp,yp,'g');
legend('data',fname,'Location','NorthWest');
xlabel('x');
ylabel('y');
print(gcf,'-dpng','regressionExample');
end