function [v,lambda] = geig(A,B)

% Find the first generalized eigenvecotr and eigenvalue of A, B:
% A*v = lambda*B*v
% A and B are square matrices, v is a vector and lambda is scalar.
    
[AA, BB, Q, Z, V, W, lambda] = qz (A, B);

[m,idx] = max(lambda);
lambda = m;
v = V(:,idx);