sun_declination <- function(day_number){
	result <- 23.45 * sin(2 * pi / 365 * (day_number + 284))
}