C语言实现计算太阳高度角程序

/**********************************************************************************************
*文件名：  精确太阳角度计算 
*程序功能：此程序计算规定规定日期的规定时间的太阳能电池位地平方位角和视高度角
*编制时间：  2011年05月13日
*程序版本： 针对太阳角度计算做了修改 
*改动部分：  对太阳赤纬角、时角等算法做了修正，加入时差等补偿
**********************************************************************************************
*输入参数：
*所测时间的年、月、日              y，m，d
*所测时间的时分秒              h,m,s
*所测地地理经度                  LONG
*所测地地理纬度                  LAT
*程序输出：                         屏幕打印
说明：  程序中输入输出借用C++输入输出库 （头文件<iostream>）
*************************************************************************************/
#include<iostream>
#include<math.h>
using namespace std;


#define  PI   3.14159265358979
/******************************************************************************
 *声明积日计算函数n_ca
 *输入参数为月：a、日：b
 *输出为当天为当年的第n天
 ******************************************************************************/
unsigned int  n_ca(unsigned int a,unsigned int b);
/*******************************************************************************
 *声明日角计算函数
 *输入参数为：年：y，积日：n

 *输出为日角 angle_day 
******************************************************************************/
double angle_day_ca(unsigned int n,unsigned int y);

/*******************************************************************************
*声明赤纬角的计算函数delt_sun_ca
*输入参数为：日角angle_day
*输出为赤纬角 delt_sun
********************************************************************************/
double delt_sun_ca(double a);

/******************************************************************************
 *声明时差计算函数
 *输入参数为日角angle_day
 *输出为时差et
*******************************************************************************/
double et_ca(double a);

/*****************************************************************************
 *声明时角计算函数
 *输入参数为:时分秒：hour，minute，sec；时差et,地理经度LONG
 *输出为时角angle_hour
*****************************************************************************/
double angle_hour_ca(unsigned int h,unsigned int m,unsigned s,double e,double l);


/****************主函数计算太阳高度角与方位角************/
int main()
{