function a = Quadrik2D(Punkte)
    x = Punkte(:,1);
	y = Punkte(:,2);

    A1 = [x.^2  y.^2 sqrt(2).*x.*y];
    A2 = [x y ones(size(x))];

    H = A2*inv(A2'*A2)*A2';
    H = eye(size(H))-H;
    H = A1'*H*A1;
    [evec eval] = eig(H);

    posMinEval = find(diag(abs(eval)) == min(diag(abs(eval))));
    
    x1 = evec(:,posMinEval);
    x2 = -inv(A2'*A2)*A2'*A1*x1;
    
    a = [x1; x2];
return;