У меня только общие соображения "в лоб". Два последовательных CR звена дадут цепь второго порядка. В этом можно сходу убедиться составив характеритсическое уравнение, а также по виду комплексной частотной характеристики (КЧХ). Последнюю всё равно придётся получить. КЧХ представим в виде:

где

- коэффициенты, определяемые структурой цепи и параметрами её элементов (как именно - придётся найти),

- частота среза. Рассматривая КЧХ в нуле, запишем, что

, откуда

. Рассматривая предел на бесконечности, получим

. Тогда выражение для КЧХ должно приводиться к виду:

и имеет три параметра. На частоте среза АЧХ должна соответствовать уровню среза

, откуда получаем связь между параметрами:

Находим АЧХ

, находим её производную на частоте среза

. Используя связь между тремя параметрами, выражаем любой из них через два другие. Ищем максимум полученной функции двух переменных.
Не исключено, что при выводе КЧХ какой-то параметров

сам по себе окажется равным нулю. Тогда всё будет очень красиво. Но это надо проверять.