Т.е. надо решить уравнение

. Ну можно методом Ньютона:

На каждом шаге придётся считать интеграл численно, но это не страшно, т.к. шагов понадобится и всего-то несколько штук.
Конечно, нужно знать подходящее начальное приближение

, но тут уж никаких универсальных рекомендаций дать невозможно, тут уж ничего не поделаешь. Можно, например, постепенно меленькими шажками накапливать интеграл по формуле просто прямоугольников, пока не перешагнёшь уровень

-- это и будет начальным приближением. А степень меленькости выбирать, исходя из специфики задачи и из соображений здравого смысла.