function pxi=NcalGCum(cgi, ngi, pgi)

% calculate cumulants of each generator
%cgi - perunit
%ngi- No. of generator
%pgi-availability

afa1=pgi;
afa2=pgi;
afa3=pgi;
afa4=pgi;
afa5=pgi;
afa6=pgi;
afa7=pgi;
afa8=pgi;

px1=afa1;
px2=afa2-px1^2;
px3=afa3-3*px1*px2-px1^3;
px4=afa4-3*px2^2-4*px1*px3-6*px1^2*px2-px1^4;
px5=afa5-5*px1*px4-10*px2*px3-10*px1^2*px3-15*px1*px2^2-10*px1^3*px2-px1^5;
px6=afa6-10*px3^2-15*px2*px4-6*px1*px5-15*px2^3-15*px1^2*px4-60*px1*px2*px3-20*px1^3*px3-45*px1^2*px2^2-15*px1^4*px2-px1^6;
px7=afa7-7*px1*px6-21*px2*px5-35*px3*px4-21*px1^2*px5-105*px2^2*px3-70*px3^2*px1-105*px1*px2*px4-35*px1^3*px4-105*px2^3*px1-210*px1^2*px2*px3-35*px1^4*px3-105*px1^3*px2^2-21*px1^5*px2-px1^7;
px8=afa8-35*px4^2-8*px1*px7-28*px2*px6-56*px3*px5-28*px1^2*px6-210*px2^2*px4-280*px3^2*px2-168*px1*px2*px5-280*px1*px3*px4-105*px2^4-56*px1^3*px5-280*px1^2*px3^2-420*px1^2*px2*px4-840*px1*px2^2*px3-70*px1^4*px4-560*px1^3*px2*px3-420*px1^2*px2^3-56*px1^5*px3-210*px1^4*px2^2-28*px1^6*px2-px1^8;

pxi=zeros(1,8);
pxi(1,1)=ngi*cgi*px1;
pxi(1,2)=ngi*cgi^2*px2;
pxi(1,3)=ngi*cgi^3*px3;
pxi(1,4)=ngi*cgi^4*px4;
pxi(1,5)=ngi*cgi^5*px5;
pxi(1,6)=ngi*cgi^6*px6;
pxi(1,7)=ngi*cgi^7*px7;
pxi(1,8)=ngi*cgi^8*px8;

return