Circuit Theory/2Source Excitement/Node and Mesh/mcode1

syms va vb vc
v1 = 5*(sqrt(3) + i);
i1 = (1-i)/(2 * sqrt(2));
Zc1 = -10i;
Zc2 = -5i;
ZL = i;
R1 = 1000;
R2 = 500;
R3=2000;
S = solve((v1-va)/(R1 + ZL) + i1 - (va - vb)/Zc1, (va-vb)/Zc2 - vb/R3 - (vb-vc)/Zc2, (vb-vc)/Zc2 - i1 - vc/ZL, va,vb,vc)
S = [S.va S.vb S.vc];
VA = S(1);
VB = S(2);
VC = S(3);
vpa(VA,3)
vpa(VB,3)
vpa(VC,3)
vpa(abs(VA),3)
vpa(abs(VB),3)
vpa(abs(VC),3)
vpa(atan(imag(VA)/real(VA))-pi(),3)
vpa(atan(imag(VB)/real(VB))-pi(),3)
vpa(atan(imag(VC)/real(VC)),3)