//计算分城乡、性别生育年龄模式
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;
       }
}