integrate exp (p * cos(x)) dx from 0 to pi