Add one-dimensional numerical integration capabilities to Stan. This would compute definite integrals. Do we allow +/- infinity as bounds here?