Dear all,

There is a bug in getPowerDeriv. In rare cases such as entering x^(alppha-1) with x=0 in steady state and alpha=2 the derivative is wrong. The reason is that the first condition in getPowerDeriv checks whether x=0 and then sets the derivative to 0 as the derivative is p*0^(p-1) which is 0 if p~=1. But this is wrong if the exponent is p=1 as it then is 1*0^0=1.

I pushed a fix to git where the first condition checks whether p=1. This also affects the 4.2 branch.

Best,

Johannes

--------

Johannes Pfeifer

Haußerstr. 29

72076 Tübingen

Tel.: +49-(0)7071-6396184

Mobil: +49-(0)170-6936820

jpfeifer@gmx.de