This function calculates the obliquity for a given epoch. It is a
wrapper for function
swe_calc of package swephR.
obliquity(year = cur.year)
Year for which to calculate the obliquity. Defaults to present year as given by Sys.Date()
Laskar, J. et al. (2004), A long-term numerical solution for the insolation quantities of the Earth, Astron. Astroph., 428, 261-285, doi:10.1051/0004-6361:20041335.
#' # Obliquity for year 3999 BC: obliquity(-4000)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.