integrate exp (a*cos(x)-b*cos(2*x)) dx from 0 to pi