function[Vpol,Iop,AngDir]=PDire(V0,I0) Vpol=-3*V0; Iop=3*I0; AngDir=angle(Vpol/Iop); Vpol=abs(Vpol); Iop=abs(Iop); AngDir=-(AngDir/pi)*180;