我想请问今年4月4日上海的日出日落时间以及正午太阳高度角
（一）计算4月4日这天太阳高度角h的值
我知道太阳高度角的计算公式为：sin h＝sinφsinδ＋cosφcosδcos15°（t-12）
（公式来源于教材《农业气象学》）.
其中φ为地理纬度,上海的纬度为北纬31°2',δ为赤纬,t为地方时,计算正午的太阳高度角,那么t=12.
而赤纬的计算公式为δ=χ*180/π,而χ为中间变量,计算公式为：
χ=0.006918-0.39912cosθ+0.70257sinθ-0.006758cos2θ+0.000907sin2θ-0.002697cos3θ+0.001480sin3θ
其中,θ=360°d/N,d为一年中日期序数,4月4日即一年中的第93天,故d=93.N为365天.
请根据以上公式帮我求得太阳高度角h的值.
（二）计算4月4日的日出日落时间
首先计算这天的时角ω=tanψtanδ,接着计算这天的日照时数T=（2*|ω|）/15°.
请帮我计算出以上公式,直到日照时数后,我自然能够推算出日出日落时间,先给100分,完成后再另行追加,