Совпадение с измерениями достаточно хорошее
Измерения       Расчёт
340                  341
205                  224

Оказывается частота не та, для которой считал, а выше - 8.6 Мгц, тогда совпадение ещё ближе Q=214, а не 224.

Лампа в момент положительной полуволны полностью заперта. А реактивности скомпенсированы П контуром. И вряд ли кому потребуется сунуть туда дополнительный резистор, что бы тупо сажать входную мощность на корпус. Ну а посчитать, проблем наверное не составит.

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

Есть много разных реализаций функции поиска минимума (или максимума, если приписать минус ). Можно поискать для языка, которым пользуетесь.

Игорь моделировал в Мультисим,где то есть на форуме.

Практически совпало

Это просто замечательно.
Я считаю прямо в лоб, прохожу с мелким шагом по частотному диапазону и ищу две точки, где спад от максимума 0.707. Потом делю частоту на полосу. Можно сделать две итерации, сначала относительно крупным шагом пройти, найдя приблизительно максимум, потом в окрестности пройтись мелким шагом. Тогда и задавать приблизительную частоту резонанса вручную не придётся.