tb = np.hstack((theta0, beta0))