Solving large system of linear equations using MATLAB

I am using MatLab to solve a system of linear regressions: Ax=b. But my matrix A is really big. Its shape is (4098,3628800). I used A\b to solve this problem. It takes about 1 hour to do. I was wondering if there are more efficient ways to do this in Matlab?