不完全ベータ関数(1-alpha,0,x)の漸近展開の実装の例 後半は部分分数分解やテーラー展開で積分 bx_real<-function(x,alpha){ z<-x if(length(x[abs(x)>1])>=1){ n<--log(0.00001,base=min(x[abs(x)>1])) str1<-0 for(i in 0:n){ str1<-str1+x[abs(x)>1]^-i/(…