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)