从日出到正午， 算是光照的上半天， 再乘上2， 就是一天的光照时间

先算赤纬角a23.45sin(360°×(284＋n)/365)，n是从1月1日到要计算那天。查到某地纬度，假设为b(a和b单位是°)算日出日落时角c＝±arccos(-tana×tanb)，±对应出落日照时间t＝2×c/15