jl> p = big(3.6e13)
3.6e+13
jl> q = 1/big(0.563483398)+1
2.774675178628776489804929213602311938188089959295862682297487944939957397047059
jl> beta(p, q)
NaN
??
jl> p = big(3.6e13)
3.6e+13
jl> q = 1/big(0.563483398)+1
2.774675178628776489804929213602311938188089959295862682297487944939957397047059
jl> beta(p, q)
NaN
??