Z1=6+j*2;
Z2=-j*4;
Z3=10;
[X,Y]=pol2cart(20/180*pi,100);
VS=X+Y*j
Zth=Z3+Z1*Z2/(Z1+Z2) % V th
[theta1,rho1]=cart2pol(real(Zth),imag(Zth))
% voltage
% Vab=Volage of -j4, -j4 and Z1 is in serial,
%using voltage division
Vth=(-j*4)/(Z1-j*4)*VS
[abs(Vth),angle(Vth)*180/pi]
![answer](https://i-blog.csdnimg.cn/blog_migrate/3ff17e2e7e8738799bc37092483b5219.png)
Zth= Z3+Z1||Z2;
Vth=voltage of Z2