Mark Dickinson added the comment: Sorry: those lines should have been: double temp = 1-z.real; printf("temp is %.19e\n", temp); r.real = log1p(4.*z.real/(temp*temp + ay*ay))/4.; __________________________________ Tracker <report at bugs.python.org> <http://bugs.python.org/issue1640> __________________________________