Esercizio - Corrente Alternata - Calcolo Thevenin in regime sinusoidale con codice OCTAVE

Breve esercizio di calcolo dell'equivalente di Thevenin in sinusoidale con calcolo impedenza equivalente #elettronica #elettrotecnica
Codice OCTAVE:
clc
clear
w=100
E1= 30*e^(i*pi/6)
J2= 6*e^(i*-pi)
E4= 10*e^(i*pi/2)
R1=3
L1= i*w*0.04
R2= 2
L3= i*w*0.05
R3= 2
C3= 1/(i*w*0.0025)
L4= i*w*0.01
R4= 4
Z1= R1+L1
Z2= R2
Z3= R3+C3
Z4= R4+L4
ZL3= Z1+L3
Z13= (ZL3*Z3)/(ZL3+Z3)
ZTH= Z13+Z4
A= [1/Z1+1/L3 -1/L3;-1/L3 1/L3+1/Z3]
B= [(E1/Z1)+J2;0]
C= inv(A)*B
VA= C(2,1)
VB= E4
ETH= VA-VB