function Lz=FormLz(Mat_G,Init_L,GenL,Busnum,PVQL,PD,PD0,QD0,Loadi,KK,PF,PDL,QDL)
KK=999;

VoltL=(0.90)*ones(1,Busnum);
%VoltL=-10*ones(1,Busnum);
% PDL=PD0(Loadi);
% PDL(PDL>0)=0.500*PDL(PDL>0).*unifrnd(0.75,1.25,sum(PDL>0),1);
% PDL(PDL<0)=1.500*PDL(PDL<0).*unifrnd(0.75,1.25,sum(PDL<0),1);
% 
% % 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));
 
% QDL=QD0(Loadi);
% 
% QDL(QDL>0)=0.500*QDL(QDL>0).*unifrnd(0.75,1.25,sum(QDL>0),1);
% QDL(QDL<0)=1.500*QDL(QDL<0).*unifrnd(0.75,1.25,sum(QDL<0),1);
% 
% % 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));
 
t1=([PDL',QDL',VoltL])';
t2=Mat_G-Init_L'-t1;
Lz=t2;

end