def beta0func(alpha, beta, TH, aTH, bTH, cosTH, tanTH, W):