Hello list,

I am using linalg.solve to solve a system of linear equations.
As I have to solve multiple systems with the same matrix, but different right-hand sides, I tried to make the right-hand side a matrix and that works fine.
So the docstring should say:

    Solve the equation ``a x = b`` for ``x``.

    Parameters
    ----------
    a : array_like, shape (M, M)
        Input equation coefficients.
    b : array_like, shape (M,) or array_like, shape (M,N)
        N can be arbitrary size
        Equation target values.

    Returns
    -------
    x : array, shape (M,) or array, shape (M,N)

Thanks,

Mark