You should use numerical_integral:
sage: f(x) = 1/(3.33*sqrt(2*pi))*exp(-1/2*((x-9.486)/3.33)^2)

The second number in the output is an estimate of the numerical error (type numerical_intergal? for more details, in particular for options to control the numerical error).