Integrate[y^a WhittakerM[k, m, y] Exp[y/2], y]