Buenas noches.


soy nuevo en esto de los foros y quisiera pedir su ayuda ya que me encuentro realizando los cálculos para el azimut y la altura solar con la siguientes formulas;


Ele=360*((dn-1)/365)*0.0174533;


deltha=(0.006918-0.399912*cos(Ele)+0.070257*sin(Ele)-0.006758*cos(2*Ele)+0.000907*sin(2*Ele)-0.002697*cos(3*Ele)+0.00148*sin(3*Ele));


Et=(0.000075+0.001868*cos(Ele)-0.032077*sin(Ele)-0.014615*cos(2*Ele)-0.04089*sin(2*Ele))*229.19;
LST=(LCT+4*(LM-LH)+Et)/60;
w=15*(LST-12)*(pi/180);
alfa=(asin(sin(deltha)*sin(fi)+cos(deltha)*cos(fi) *cos(w)))*(180/pi);
psi=(acos((sin(alfa*(0.0174533))*sin(fi)-sin(deltha))/(cos(alfa*(0.0174533)*cos(fi)))))*(180/pi);


donde:
ele= angulo diario
deltha= declinacion solar
Et= ecuacion del tiempo en minutos
LST=tiempo solar verdadero
W=angulo diario
alfa=altura solar
psi=azimut


mi problema radica en que a la hora de realizar los cálculos, los grados obtenidos no son los correctos, para comprobar dichos grados ocupo la pagina NOAA SOLAR CALCULATOR, si alguien me pudiera ayudar y decirme que estoy haciendo mal se los agrade seria mucho.


PSD. Estoy trabajando en las formulas en matlab.