#Cuanta energía llega la la parte superior de la atmósfera dependiendo del día, a una superficie perpendicular a la radiación Solar el día (n) por unidad de area por unidad de tiempo.

library(ggplot2)

DIA_n<- 1

ALPHA<- 2*pi*DIA_n/365

CONSTANTE_SOLAR<- 1353

ENERGIA_LLEGA_DIA<- CONSTANTE_SOLAR*(1+ 0.034*cos(ALPHA))
ENERGIA_LLEGA_DIA
## [1] 1398.995
#Vamos a graficar

DIA<- 1:365
ALPHA_n<- 2*pi*DIA/365

Grafico_ENERGIA_LLEGA_DIA<- CONSTANTE_SOLAR*(1+ 0.034*cos(ALPHA_n))

Data<-data.frame(DIA,Grafico_ENERGIA_LLEGA_DIA)

ggplot(data= Data,aes(DIA,Grafico_ENERGIA_LLEGA_DIA)) + geom_line()