function ret = falpha(alpha)

ret = -3*cos(-1/2-(1/8+8802709594530323/72057594037927936*alpha)*(1/2-1253749088136497/4503599627370496*alpha))*(-3787713241984335/144115188075855872+11036389127272785701120029498531/162259276829213363391578010288128*alpha)*cos(1/2+8802709594530323/18014398509481984*alpha)+26408128783590969/18014398509481984*sin(-1/2-(1/8+8802709594530323/72057594037927936*alpha)*(1/2-1253749088136497/4503599627370496*alpha))*sin(1/2+8802709594530323/18014398509481984*alpha);