Switch PyMC3 to PyMC (which is the new name starting from PyMC v4.0.0)

This MR updates the PyMC3 interface to use the new version of PyMC, which is now v4.0.0 (the version number suffix has been dropped). This is mainly a very straightforward port by replacing imports of pymc3 with pymc. The Theano library has also been renamed to aesara, so those are changed too.

The test suite runs and also the two examples, linear_regression_pymc.py and linear_regression_pymc_custom_likelihood.py, both work exactly as expected.

examples/core_examples/gaussian_process_celerite_example.py

@@ -6,7 +6,6 @@ import matplotlib.pyplot as plt
 import numpy as np
 from bilby.core.prior import Uniform
 
-
 # In this example we show how we can use the `celerite` package within `bilby`.
 # We begin by synthesizing some data and then use a simple Gaussian Process
 # model to fit and interpolate the data.