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