I didn't have a MATLAB, so I tried to running it on Octave.
I was running it in this way:
octave:1> DI_PSEUDOCOUNT_WEIGHT = 0.05
DI_PSEUDOCOUNT_WEIGHT = 0.050000
octave:2> DI_THETA = 0.62
DI_THETA = 0.62000
octave:3> calculate_evolutionary_constraints('PF02142.selected','PUR9_CHICK','salida.txt');
But I get all 0.0 as result.
How can I use it on Octave?
Thanks!!!