Sean $f_1, f_2, ......., f_n$ funciones linealmente independientes. Construimos la función de prueba $\varphi$ combiando linealmente la $f_i$ funciones.

\begin{equation}\label{ec5}
\varphi=c_1 f_1 + c_2 f_2 +........+c_n f_n =\sum_{j}c_i f_i
\end{equation}

Donde $\varphi$ es la función variacional de prueba, $c_i$ son parámetros variacionales que deben hacer mínima la integral variacional. $c_i$ y $f_i$ son reales.

Planteamos la integral variacional:

\begin{equation}\label{ecu6}
W=\frac{\int\varphi^{\ast}\hat{H}\varphi dq}{\int\varphi^{\ast}\varphi dq}
\end{equation}

Calculamos por separado numerador y denominador de la integral variacional

\begin{equation}\label{ecu7}
\int\varphi^{\ast}\varphi dq=\int\sum_{i}c_{i}^{\ast}f_{i}^{\ast}\sum_{j}c_j f_j dq=\sum_{i}c_{i}^{\ast}\sum_{j}c_{j}\int f_{i}^{\ast}f_j dq=\sum_{i}\sum_{j}c_i c_j S_{ij}
\end{equation}

Donde $S_{ij}$ es la integral de solapamiento.

\begin{equation}\label{ecu8}
\int\varphi^{\ast}\hat{H}\varphi dq=\int\sum_{i}c_{i}^{\ast}f_{i}^{\ast}\hat{H}\sum_{j}c_{j}f_{j}=\sum_{i}\sum_{j}c_{i}^{\ast}c_j\int f_{i}^{\ast}\hat{H}f_j dq=\sum_{i}\sum_{j}c_{i}^{\ast}c_j H_{ij}
\end{equation}

Sustituyendo las ecuaciones (\ref{ecu7}) y (\ref{ecu8}) en (\ref{ecu6}):

\begin{equation}\label{ecu9}
W=\frac{\sum_{i}\sum_{j}c_{i}c_{j}H_{ij}}{\sum_{i}\sum_{j}c_{i}c_{j}S_{ij}}
\end{equation}

Vamos a minimizar W, para aproximarnos tanto como sea posible a $E_0$. Para ello, derivamos $W$ respecto a los coeficientes $c_k$ de la combinación lineal, e igualamos a cero las derivadas. $\frac{\partial W}{\partial c_k}=0$ con $k=1,2.....n$, da como resultado: $\sum_{i=1}^{n}c_{i}\left(H_{ki}-WS_{ki}\right)=0$

Para una función de prueba con dos funciones base $\varphi=c_1f_1+c_2f_2$, obtenemos el siguiente sistema:
\begin{equation}
c_1(H_{11}-WS_{11})+c_2(H_{12}-WS_{12})=0\\
c_1(H_{21}-WS_{21})+c_2(H_{22}-WS_{22})=0
\end{equation}
Para que el sistema tenga solución única (compatible determinado), el determinante de los coeficientes debe anularse.
\begin{equation}
\left| \begin{array}{cc}
H_{11}-WS_{11} & H_{12}-WS_{12}\\
H_{21}-WS_{21} & H_{22}-WS_{22}
\end{array}
\right|=0
\end{equation}
La resolución del determinante nos da los valores propios de $\varphi$, es decir, las energías, W. A partir de las energías se obtienen los coeficientes variacionales $c_1$ y $c_2$.