Ax=b solver