JuliaLang/ODE.jl — pull request #52 by yuyichao (2014-12-26 01:34:01+00:00)
Title: 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 return a polynomial with constant term `k` default to `0`. Just calculating the value at `1` is much simpler and returns the same result.