tb0 = np.vstack((theta00, beta00))