powered by
sw.beta(S, t, p, is.theta=FALSE)
t
p
ctd
S
sw.beta(40, 10, 4000, is.theta=TRUE) # 0.72088e-3 (from p1964 of McDougall 1987)
Run the code above in your browser using DataLab