Best way Solve linear system of sparse symetric matrix

Hello everyone i have a linear system of size n X n, which i get from c++, then load it to matlab and get the resutls back. To solve the linear system i use lsqr(A,b). My matrix is Sparse and Symetric. Is there a better way to do it?

Thank you all