nday=200
Rdc=((360/365)*nday)*(pi/180)
dec=0.33281-22.984*cos(Rdc)-0.3499*cos(2*Rdc)-0.1398*cos(3*Rdc)+3.7872*sin(Rdc)+0.03205*sin(2*Rdc)+0.07187*sin(3*Rdc)*(pi/180)

print(dec)

if (nday < 21){
    teq = 2.6+0.44*nday
} else if (21 <= nday & nday <= 136){
    teq = 5.2+9.0*cos((nday-43)*0.0357)
} else if (136 <= nday & nday <= 241){
    teq = 1.4-5.0*cos((nday-135)*0.0449)
} else if (241 < nday & nday <= 336){
    teq = (-6.3)-(10*cos((nday-306)*0.036))
} else {
    teq = 0.45*(nday-359)
}

TZ=12
Lw=10*(pi/180)
Tshift=TZ-(Lw/15)

nhour=12
ts=(nhour-teq)/(60-Tshift)

w=((180/12)*(12.5-ts))*(pi/180)

alt=asin((sin(dec)*sin(Lw)+cos(dec)*cos(Lw)*cos(w)))

print(alt*(180/pi))

y=180*(pi/180)
b=90*(pi/180)

inc=acos(sin(dec)*sin(y)*cos(b)-sin(dec)*cos(Lw)*sin(b)*cos(y)*cos(w)+cos(dec)*sin(b)*sin(y)*sin(w))
print(inc)