Custom theano Op to do numerical integration

That is a normal cdf, its from pymc3.distributions.dist_math