This function uses the following argument for output:
b - On output, b is overwritten by the solution to the matrix equation ax=b.
ier (Fortran only) - When called from a Fortran program, this function returns error status in ier.