Without dependent (voltage) sources one could normally define a matrix equation to solve for the values in the network:
[Z][I] = [V]
In this case though we must add in dependent voltage sources
[Z][I] = [V] + [k][I]
Apologies, I renamed I3 above as I2 below, per the usual custom of ordering currents left to right as shown in sketch.
From this diagram, we extract the matrix equation suggested above.
This rearranges as
[Z][I] - [k][I] = [V]
Rearranging again [I] = [V] · ([Z]-[k])^-1
Solution of the matrix math is standard matrix math left to the reader. It produces an equation for each of the 3 currents defined in the sketch as a function of the single independent voltage.
This gives the currents as a function of the one independent voltage.
At this point
Va = (I1 - I2)·R2
And Vb = (I2 - I3)·R4