Escribir las ecuaciones de malla y nodal para circuito grande y supermalla grande

Aquí está mi circuito:

^{simular este circuito : esquema creado con CircuitLab}

He definido los 'voltajes' del nodo y las direcciones de corriente de malla y las corrientes mismas. Lo que quiero hacer es calcular las ecuaciones nodales y de malla para este circuito dada la información anterior.

Aquí está mi trabajo en el análisis nodal:

$$KCL\;\;V_1:\;\;{ {{V_1-V_2}\over R_2} +{{V_1-V_3}\over R_1} + {V_1\over R_3}= 0}$$ $$KCL\;\;V_2:\;\;{ {-I_b} +{{V_2-V_1}\over R_2} -I_a = 0}$$ $$KCL\;\;V_3:\;\;{ {{V_3-V_1}\over R_1}+I_a +{{V_3-V_4}\over R_4} = 0}$$ $$KCL\;\;V_4:\;\;{ {{V_4-V_3}\over R_4} +{{V_4}\over R_5} - I_c = 0}$$ Por lo que puedo decir, esos deberían ser correctos.

Mi próximo problema es cómo obtener las ecuaciones para las corrientes de malla. Mi problema es que está en la(s) supermalla(s) formada(s) por Ia e Ib. No estoy seguro de cómo obtener las corrientes de malla si están ahí.

Para resumir, estoy bastante seguro de que mis ecuaciones de voltaje/nodo son correctas, pero no estoy 100% seguro de eso. No sé cómo obtener las ecuaciones de corriente de malla si esas fuentes independientes están en esta configuración particular.

Si alguien pudiera ayudar, se lo agradecería mucho. ¡Gracias!