Numerical integration using interval arithmetic gets you the same thing but in a completely rigorous way.