function Mat_G=FormG(Volt,PD,QD,Loadi,Vbi,mVolt,RealVolt,bigM,sigma,PDbi,RealPD,RealQD)

Mat_G=[
    sparse(PD(Loadi))-bigM*PDbi-PD(Loadi)*sigma;
    sparse(PD(Loadi))+bigM*PDbi+PD(Loadi)*sigma;
    sparse(QD(Loadi));
    Volt'-bigM*Vbi-mVolt'-sigma*RealVolt';
    Volt'+bigM*Vbi+mVolt'+sigma*RealVolt';
    Vbi;
    PDbi;
    ];
end