Решатель Гаусса для БИХ-фильтра

arduino-ide arduino-due

Я пишу фрагмент кода в онлайн-GDB, используя C++, код вычисляет некоторые коэффициенты для БИХ-фильтра, который пытается аппроксимировать дробную производную и интеграл для дробного ПИД-регулятора. Код работает, как и ожидалось, в онлайн-IDE, но когда я пытаюсь выполнить его в Arduino Due, результаты только переполняются и нан, и я действительно не знаю, в чем проблема, поэтому,

Вот код в GDB онлайн https://onlinegdb.com/rycwt62XL, который дает ожидаемые результаты:

Массив А, Массив Б

1,41421,1

-0,707107,0,5

-2,12132,-1,5

0,883883,-0,625

0,883883,0,625

-0,265165,0,1875

-0,0883883,-0,0625

0,0110485,-0,0078125

А вот код для Arduino Due, который только изменяет метод печати с помощью последовательных команд, https://onlinegdb.com/ByAs5xyVU, и результаты:

Массив А, Массив Б

ovf,1.00

овф,овф

ovf,0.00

inf,ovf

ovf,0.00

inf,ovf

ovf,3,89

inf,ovf

И действительно не знаю, с чего начать, чтобы исправить эти результаты, спасибо. Метод точно такой же.

    double* h = new double[L]; 
    for(int n = 0; n <= L; n++){
        for(int k = 0; k <= n; k++){
            h[n] += pow(operador,-alpha-n+k)*pow(1.0-operador,n-k)*Binomial_E(alpha,k)*Binomial(-alpha,n-k);
        }
    }

double h[L+1];