Searched refs:sinsun (Results 1 – 1 of 1) sorted by relevance
2162 double sinsun, sinpos, sinapp, num; local2168 sinsun = sin(torad(sundec));2173 num = 1 - sinsun*sinsun - sinpos*sinpos - sinapp*sinapp2174 - 2*sinsun*sinpos*sinapp;2176 if (sinsun*sinpos>0) 2181 duration = 12.0 + 24.0*atan((sinsun*sinpos+sinapp)/sqrt(num))/PI;