void CalFertilityPattern()
{
//TODO: Add your control notification handler code here
inti,j;
for(i=0;i<MAX_AGE;i++)
for(j=0;j<2;j++)
m_FertilityRatio[i][j]=0.0;
for(i=m_ufR1;i<=49;i++)
{
m_FertilityRatio[i][0]=KafangFun(i,m_ufR1,m_ufN);
}
for(i=m_sfR1;i<=49;i++)
{
m_FertilityRatio[i][1]=KafangFun(i,m_sfR1,m_sfN);
}
}
double KafangFun(int i,int R1,int N)
{
returnpow((i-R1),((double)N/2-1))*
exp(-(double)(i-R1)/2)/(pow(2,((double)N/2))*calGMfun(N));
}
//计算GM值
double calGMfun(int N)
{
doubleGM=1,SG=(double)N/2;
if((SG-1)<0.000001)
{
returnGM;
}
else
{
while((SG-0.5)>0.000001)
{
SG=SG-1;
GM=GM*SG;
if((SG-1)<0.000001)
return GM;
}
GM=GM*sqrt(3.141592653);
return GM;
}
}
