def f の中で途中結果を print しまくってみてください。 まずは N=1 とか 2 とかで試すといいですよ。