EIG Function
General and symmetric matrices
EIGSYMGEN Function
Real symmetric matrices and B positive definite
GENEIG Procedure
Computes the generalized eigenexpansion of a system Ax = lBx.