int(mu^(K+a-1)*exp((-N+b)*mu),mu,0,inf)