Rでガウスの超幾何関数を計算する。 hypergeo 関数は一部パラメータでバグがあるのでgslの超幾何関数でおきかえる。 ただし、gslの超幾何関数は複素数に解析接続されていないので、そこはhypergeoを使う。 library("hypergeo") library(gsl) library("BAS") …

RでHurwitz's zeta関数 フルビッツゼータ関数の近似 よいライブラリがないので、とりあえず、つくってみた。 library("Rmpfr") Hurwit_zeta_approx<-function(beta4,x){ #N<-10^7 #k<-0:N #s<-beta4 ans2<-NULL if(beta4==beta4){ for(j in 1:length(x)){ a3…