Implementing a more sophisticated pymc3 model

Thanks @junpenglao. I have posted a new topc here where I dump the output from trying to run one of the simple polynomial integrals in the linked thread.