function Lz=FormLz(Mat_G,Init_L,Loadi,rPD3P,rQD3P,mVoltABCV,PDL,QDL)
% KK=999;
% 
% VoltL=(0.9)*ones(1,Busnum);
% %VoltL=-10*ones(1,Busnum);
% PDL=PD0(Loadi);
% PDL(PDL>0)=0.800*PDL(PDL>0);
% PDL(PDL<0)=1.200*PDL(PDL<0);
% PDL(PDL==0)=-0.400;
% realPD=PD0(Loadi);
% indPD=find(realPD>0);
% PDL(indPD(3:12:end))=0.95*realPD(indPD(3:12:end));
% PDL(indPD(9:12:end))=0.45*realPD(indPD(9:12:end));
% %PDL=-10*ones(length(Loadi),1);
% QDL=QD0(Loadi);
% QDL(QDL>0)=0.800*QDL(QDL>0);
% QDL(QDL<0)=1.200*QDL(QDL<0);
% QDL(QDL==0)=-0.400;
% realQD=QD0(Loadi);
% indQD=find(realQD>0);
% QDL(indQD(3:12:end))=0.95*realQD(indQD(3:12:end));
% QDL(indQD(9:12:end))=0.95*realQD(indQD(9:12:end));
% QDL=0*PD(Loadi).*sqrt((1-PF.^2))./PF;
% PDL=0*ones(length(Loadi),1);
% QDL=0*ones(length(Loadi),1);

% PDL=0.70*rPD3P.*(1+unifrnd(-0.15,0.15,length(rPD3P),1));
% QDL=0.70*rQD3P.*(1+unifrnd(-0.15,0.15,length(rQD3P),1));

% t1=([PDL',QDL',0.88*ones(length(Volt),1)'])';
t1=([PDL',QDL'])';
t2=Mat_G-Init_L'-t1;
Lz=t2;

end