function[Zmedreal,Zmedimag]=zmed(Zmedida); % Zmedida pu, zbase Ohm,Zmedreal Ohm, Zmedimag Ohm Zmedreal=real(Zmedida); if Zmedreal>=10000 Zmedreal=exp(1000); elseif Zmedreal<=-10000 Zmedreal=log(0); end Zmedimag=imag(Zmedida); if Zmedimag>=10000 Zmedimag=exp(1000); elseif Zmedimag<=-10000 Zmedimag=log(0); end