import math
from math import floor
import time
from datetime import date, timedelta
import datetime


LAT=0
MONTH=3
DAY=21

DA = (datetime.date(2015,MONTH,DAY).timetuple().tm_yday*(360/365.25))*(math.pi/180)
DEC = math.asin(0.3978*math.sin(DA-1.4+0.0355*math.sin(DA-0.0489)))
ALT = math.asin((math.sin(LAT*(math.pi/180))*math.sin(DEC))+((math.cos(LAT*(math.pi/180))*math.cos(DEC))))

print('DEC = %.2f' % ((DEC)*(180/math.pi)))
print('ALT = %.3f' % ((ALT)*(180/math.pi)))