В работе рассматривается плоский биллиард, ограниченный эллипсом, в поле потенциальной силы. Была найдена явная формула полиномиального потенциала, сохраняющего интегрируемость такого биллиарда. Для него была изучена структура слоения Лиувилля на всех неособых уровнях энергии с помощью метода разделения переменных. А именно, был предложен алгоритм, который строит бифуркационную диаграмму, а также инвариант Фоменко-Цишанга, исходя из значений параметров потенциала. Кроме того, была изучена топология изоэнергетического многообразия и обнаружены случаи динамики твердого тела, лиувиллево эквивалентные нашему биллиарду.