А не велико ли подавление 2-ой гармоники считает ваша программа?
При написании использованы формулы приведенные в статье Шульгина по расчету П-контура, в отношении уровня подавления они совпадают с формулами приведенными в справочнике Бунина Яйленко. В Вашем расчете указано, что подавление 2-й гармоники составляет 47,7 дБ, это при нагруженной добротности 17,2. Здесь принимается подавление 6Q для П-контура с переводом в дБ. 20log(6*17,2)=40,27, это подавление самим контуром. К этому значению прибавляется 7,45 дБ, эта величина вычисляется через коэффициенты Берга для заданного угла отсечки, в данном случае 90 градусов, и представляет собой отношение величин токов 1-й и 2-й гармоник, 20log(0,5/0,212)=7,45 дБ. Итого получили 47,7 дБ. Далее считается отношение мощности 1-й и 2-й гармоник с учётом данного значения, 656/10^(0,1х47,7)=11 мВт. В общем рассчитываемые этой программой данные при тестировании сверялись с данными которые давала оригинальная программа Шульгина, и они совпадают, если не брать в расчет несущественные различия. Хотя вообще в самой программе Шульгина расчета мощности 2-й гармоники не было, но формулы выбора необходимой добротности контура (это было в Радио 5 и 7 за 1985 год) заложены именно эти. Ну а в целом добротность 17,2 все же достаточно высока. Только у Вас в расчете оптимальное значение C1 совсем маленькое, всего 12 пф, скорее всего Вы задали только ёмкость анодов ламп, без ёмкости монтажа.