Matrix solver for large system