Figured out that epsilon is not 1 at N=0, so that's a problem. Still working on a fix.