%
%
%    Beispiel zum Schnittpunkt vieler Geraden
%
%    R. Unger ... 05/2018
%


clear;
close all;


N=5;


rand('seed',333);
a=2*(rand(N,1)-0.5);
b=rand(N,1)*5;

A=[ -a , ones(N,1)]


figure(1);
hold on
for i=1:N
 x1=-10;
 x2= 10;
 y1=a(i)*x1+b(i);
 y2=a(i)*x2+b(i);
 plot( [x1,x2],[y1,y2]);
end


% x= A\b

M=A'*A;
c=A'*b;
x=M\c;

plot(x(1),x(2),'r*')
xlabel('x')
ylabel('y')
t=sprintf('Schnittpunkt-Approximation von N Geraden mit N=%d',N)
title(t)

print -dpng 'x-n-geraden.png'

