Simplify coefficient calculation

The `diff(polyval(..., 0:1))[1]` is used to calculate the integral over `[0, 1]` but is not really necessary since `polyint(::Poly, k::Number=0)` always returns a polynomial with constant term `k` default to `0`. Just calculating the value at `1` is much simpler and returns the same result.