function value = rayleigh2 ( n, a, x, y )

%*****************************************************************************80
%
%% rayleigh2() returns the generalized Rayleigh quotient.
%
%  Formula:
%
%    RAYLEIGH2 = X' * A * Y / ( X' * Y )
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    20 October 2007
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    integer N, the order of A.
%
%    real A(N,N), the matrix.
%
%    real X(N), Y(N), the vectors used in the 
%    Rayleigh quotient.
%
%  Output:
%
%    real VALUE, the Rayleigh quotient of A and X.
%
  value = ( x(1:n) * a(1:n,1:n) * y(1:n)' ) / ( x(1:n) * y(1:n)' );

  return
end
