该文章介绍了计算双重非中心Beta分布概率密度函数的鞍点近似方法。代码来源于Marc Paolella的中间概率生成，详见2007年第370页清单10.11。由于鞍点积分，归一化略高于1的近似值。推荐使用四边形标准化。使用方法如下：f = SPncbetapdf(bvec, n1, n2, theta1, theta2)，同时需要下载SPncfpdf.m。