This example computes all the eigenvalues of $Az=\lambda Bz$, where Here $B$ is symmetric positive definite and must first be factorized by dpotrf. The program calls dsygst to reduce the problem to the standard form $Cy=\lambda y$; then dsytrd to reduce $C$ to tridiagonal form, and dsterf to compute the eigenvalues.