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