Estás haciendo esto demasiado difícil.
No es necesario insertar todas las corrientes en el diagrama y escribir expresiones para cada una. Mejor es, en cada nodo, escribir KCL directamente usando: suma de corrientes fuera = 0.
Por lo tanto, nodo \ $ \ small V_A \ $ :
$$ \ small 4+ (V_A-V_B) + \ frac {V_A} {j} = 0 $$
nodo \ $ \ small V_B \ $ :
$$ \ small (V_B-V_A) + \ frac {(V_B-V_C)} {- j} +2 (V_C-V_B) = 0 $$
nodo \ $ \ small V_C \ $ :
$$ \ small -4+ \ frac {V_C-V_B} {- j} + \ frac {V_C} {2} = 0 $$
Resuelva estos a la vez para \ $ \ small V_C \ $, luego divida por \ $ \ small 2 \ $ para \ $ \ small V_O \ $:
$$ \ small V_O = \ frac {V_C} {2} = \ frac {4} {5} (3-j) = 2.53 \ large \ angle \ small-18.4 ^ o $$
Veo que la respuesta de tu modelo da un ángulo de fase de \ $ \ small 71.6 ^ o \ $. Revisaré mi análisis nuevamente, ya que su respuesta es arctan (3), mientras que la mía es arctan (-1/3).