def _toeplitz_real_sym(a, b): """ Solve linear system Ax=b for real symmetric Toeplitz matrix A. :param a: first row of Toeplitz matrix A :param b: vector b :return: x=A^-1*b """ return sto_sl(np.hstack((a, a[1:])), b, job=0)