def METHOD_NAME(n_chains):
    return multivariate_normal.rvs(mean=[0.015, 500, noise],
                                   cov=np.diag([10, 10000, noise]),
                                   size=nchains)