acosd(sind(d)*cosd(eps)-sind(a)*cosd(d)*sind(eps)) where a=(12+(51+26.00/60)/60)*15,d=27+(7+42.00/60)/60,eps=23+(26+21.406/60)/60