powered by
ab <- sw.alpha.over.beta(S, t, p, is.theta=FALSE);
# 0.3476 (from p1964 of McDougall 1987) ab <- sw.alpha.over.beta(40, 10, 4000, is.theta=TRUE)
Run the code above in your browser using DataLab