% Vstupní hodnoty a základní výpoèty parametrù
l=200
Rho=0.0178
S=10
Un=0.4
Pnm=0.03
Eta=0.92
CosFi=0.85
Xp1=0.302

Rp=Rho*l/S
Xp=Xp1*l/1000
Pn=Pnm/Eta
Sn=Pn/CosFi
In=Sn/(sqrt(3)*Un)
Ufn=Un/sqrt(3)
Fi=acos(CosFi)

% Základní metoda výpoètu
dUf=Rp*In*CosFi+Xp*In*sin(Fi)
dUp=dUf/Ufn*100

% Alternativnì s využitím sdružených hodnot napìtí
dU=dUf*sqrt(3)
dUp=dU/Un*100

% Alternativnì s využitím výkonových parametrù
Qn=Sn*sin(Fi)
dUf=(Rp*Pn+Xp*Qn)/(3*Ufn)
dUp=(Rp*Pn+Xp*Qn)/Un^2*100

% Pøesný výpoèet
Delta=(Xp*In*CosFi-Rp*In*sin(Fi))^2/(2*Ufn)
dUf=Rp*In*CosFi+Xp*In*sin(Fi)+Delta
dUp=dUf/Ufn*100

% Alternativní pøesný výpoèet s využitím výkonových parametrù
Delta=(Xp*Pn-Rp*Qn)^2/(18*Ufn^3)
dUf=(Rp*Pn+Xp*Qn)/(3*Ufn)+Delta
dUp=((Rp*Pn+Xp*Qn)/Un^2 + (Xp*Pn-Rp*Qn)^2/(2*Un^4))*100

% Procentní chyba nepøesného výpoètu
dDeltaProc=Delta/dUf*100
