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

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