Здравствуйте!
Помогите, пожалуйста разобраться с задачей. Требуется найти

,

в замкнутом виде.
Известно, что

,

,

,

- четное. Ясно, что можно считать

при нечетных

.
Из других соображений известно, что если

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

.
После некоторых преобразований получается уравнение на

:

.
Далее, используя mathematica, находим решение однородного уравнения:

Что делать дальше - неясно. Mathematica зависает, если её просить решать неоднородное уравнение.
А как методом вариации постоянной искать решение, где 2 постоянных - тоже вопрос. Может быть функция Грина поможет, но какие краевые условия тогда брать...
Подскажите, что делать дальше.