context('Numerical Integration') test_that('Numerics get generated by antiD()',{ ff = antiD(a*exp(x^2) ~ x, a=100) expect_equal(diff(ff(x=c(-2,2))), 3290.526, tol=0.01) })