function ddgzw=func_ddg(PGi,PVi,Busnum,ContrlCount,RestraintCount,Loadi,PD,QD,Init_Z,Init_W)
CCVbiZW=Init_Z(length(Loadi)*4+Busnum*2+1:length(Loadi)*4+Busnum*2+Busnum)+Init_W(length(Loadi)*4+Busnum*2+1:length(Loadi)*4+Busnum*2+Busnum);
ddgzw=[
sparse(length(Loadi)*2+Busnum*2,ContrlCount);
sparse(Busnum,length(Loadi)*2+Busnum*2),sparse(-2*eye(Busnum))*diag(CCVbiZW),sparse(Busnum,length(Loadi)*2);
sparse(length(Loadi)*2,ContrlCount);
];
end