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

syms i2 i3 vis
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(i1*R2 + (i1+i2)*Zc1 + (i1+i3)*Zc2 - vis, v1 + i2*R1 + (i2 + i1)*R3 + i2*ZL, (i3 + i1)*Zc2 + i3*ZL + (i3 - i2)*R3, i2,i3,vis)
S = [S.i2 S.i3 S.vis];
I2 = S(1);
I3 = S(2);
VIs = S(3);
vpa(I2,3)
vpa(I3,3)
vpa(VIs,3)
vpa(abs(I2),3)
vpa(abs(I3),3)
vpa(abs(VIs),3)
vpa(pi() + atan(imag(I2)/real(I2)),3)
vpa(pi() + atan(imag(I3)/real(I3)),3)
vpa(atan(imag(VIs)/real(VIs)),3)