Gauss solver for IIR filter
I can't find the root cause of the problem - it's probably because the Due is possibly the worst board that Arduino ever made, and the compiler / linker / support libraries barely work at the best of times.

However, some things you need to look at in your code:

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


There you allocate L (15) doubles in an array. But then you go on to access array slices 0-15. That's 16 entries in an array of size 15, so you're overflowing it and corrupting other memory.

Also you have no delete to balance the new, so you have a memory leak.  You would be better using:

double h[L+1];


instead of new.