Su esquema y las direcciones actuales elegí tomar en mi cabeza:
simular este circuito – Esquema creado mediante CircuitLab
Tengo la misma ecuaciones que hacer. Mi trabajo está a la izquierda, y la de los tuyos (o más cerca de la suya) está a la derecha:
$$\pequeño\begin{align*}
-I_1\:R_1 - \left(I_1-I_3\right)R_2+V_1-V_2 &=0\:\text{V} &
I_1\:R_1 + \left(I_1-I_3\right)R_2 &=V_1-V_2\\\\
V_2-\left(I_2-I_3\right)R_4-I_2\left(R_5+R_6\right) &= 0\:\text{V} &\left(I_2-I_3\right)R_4+I_2\left(R_5+R_6\right)&=V_2\\\\
-V_1-\left(I_3-I_1\right)R_2-I_3\:R_3 -\left(I_3-I_2\right)R_4&=0\:\text{V} & \left(I_3-I_1\right)R_2+I_3\:R_3 +\left(I_3-I_2\right)R_4&=-V_1
\end{align*}$$
Estos son el resultado de la misma matriz formada:
$$\left[\begin{array}{ccc}
R_1+R_2 & 0 & -R_2\\\\
0 & R_4+R_5+R_6 & -R_4\\\\
-R_2 & -R_4 & R_2+R_3+R_4
\end{array}\right]\left[\begin{array}{c}I_1\\\\I_2\\\\I_3\end{array}\right]=\left[\begin{array}{c}V_1-V_2\\\\V_2\\\\-V_1\end{array}\right]$$
Si quieres resolver esto por la mano, la regla de Cramer es que a menudo se enseña. Aquí:
$$\tiny\begin{align*}
I_1 &= \frac{\text{det}\left[\begin{array}{ccc}V_1-V_2 & 0 & -R_2\\
V_2 & R_4+R_5+R_6 & -R_4\\
-V_1 & -R_4 & R_2+R_3+R_4\end{array}\right]}{\text{det}\left[\begin{array}{ccc}R_1+R_2 & 0 & -R_2\\
0 & R_4+R_5+R_6 & -R_4\\
-R_2 & -R_4 & R_2+R_3+R_4\end{array}\right]}&
I_2 &= \frac{\text{det}\left[\begin{array}{ccc}R_1+R_2 & V_1-V_2 & -R_2\\
0 & V_2 & -R_4\\
-R_2 & -V_1 & R_2+R_3+R_4\end{array}\right]}{\text{det}\left[\begin{array}{ccc}R_1+R_2 & 0 & -R_2\\
0 & R_4+R_5+R_6 & -R_4\\
-R_2 & -R_4 & R_2+R_3+R_4\end{array}\right]}&
I_3 &= \frac{\text{det}\left[\begin{array}{ccc}R_1+R_2 & 0 & V_1-V_2\\
0 & R_4+R_5+R_6 & V_2\\
-R_2 & -R_4 & -V_1\end{array}\right]}{\text{det}\left[\begin{array}{ccc}R_1+R_2 & 0 & -R_2\\
0 & R_4+R_5+R_6 & -R_4\\
-R_2 & -R_4 & R_2+R_3+R_4\end{array}\right]}
\end{align*}$$
Así que ¿has probado este enfoque?
Un estudio reciente, Una condensación basado en la aplicación de la regla de Cramer para resolver problemas de gran escala de sistemas lineales, Habgood & Arel, Diario de Discreto Algoritmos de 10, 2012, pp 98-109, demuestra que la regla de Cramer puede realizar en el mismo computacional de la orden de los otros métodos, tales como la descomposición LU. Es agradable leer sobre eso.
Suponiendo que se ha trabajado en los tres corrientes anteriores, ahora usted puede fácilmente averiguar la magnitud de la tensión a través de \$R_5\$ (si se me entiende el uso de \$V_x\$, correctamente) y que pueda trabajar fuera de su polaridad con respecto a la presunta actual, \$I_2\$, así.
Creo que usted sabe cómo calcular la corriente en \$V_1\$ como la suma de dos de las tres corrientes que trabajó anteriormente. Y del mismo modo, para \$V_2\$. A partir de esas corrientes, y el conocido de los voltajes de cada uno, que no debería tener problemas para trabajar el poder que contribuyen (la señal que le indica si se están generando (-) o disipe ( + )).
Estos últimos detalles no debería ser difícil, una vez que tenemos las tres corrientes trabajado.