Circuit Theory/2Source Excitement/Thevenin and Norton/code1

v1 = 5*(sqrt(3) + i);
i1 = (1-i)/(2 * sqrt(2));
Zc1 = -10i;
Zc2 = -5i;
ZL = i;
R1 = 1000;
R2 = 500;
R3=2000;
iT = v1/(2*ZL + R1 + Zc1 + Zc2);
VAB1 = iT * (Zc2 + ZL);
vpa(VAB1,3)
VN = i1*(1/(1/(R1 + 2*ZL)) + 1/(Zc1 + Zc2));
VAB2 = VN*Zc2/(Zc1 + Zc2) - VN*ZL/(R1 + 2*ZL);
vpa(VAB2,3)
VTH = VAB1 + VAB2;
vpa(VTH,3)