quadrature1d_modified.m evaluates line integral in 2d