dvonmises(x, mu, kappa)
pvonmises(q, mu, kappa, tol = 1e-20)
rvonmises(n, mu, kappa, units = c("radians", "degrees"), ...)circular.circular.dvonmises gives the density, pvonmises gives the
distribution function and rvonmises generates random deviates.