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