Added basic matrix vector product examples
See #1531
Hello :)
This is actually the most embarrasing commit to do, but here
it goes.
BUG: allow .moment to take shapes as keywords.
BUG: allow also rv_discrete.moment to take shape as keywords.

Also fix a test issue for moment() in the continuous tests.
Merge branch 'moment' into master.

Contains relevant commit of https://github.com/scipy/scipy/pull/2753 and an
extra bugfix.
MAINT: small error message change in linalg.solve_triangular. Closes gh-2664.
Merge branch 'pr/2749' into master.

https://github.com/scipy/scipy/pull/2749
DEP: linalg: deprecate expm2 and expm3

The Pade approximant provided by expm is a superior alternative, so
providing two worse methods just clutters the documentation.
DOC: linalg: explain the Taylor order used in expm3

Better to fix this in the documentation than to change the definition.

Fixes gh-2212
DOC: update release notes on deprecations
DOC: release: update release notes vs. LinearOperator changes
Merge pull request #2738 from pv/deprecate-expm23

 DEP: linalg: deprecate expm2 and expm3
TST: sparse: more zero-size matrix tests
Merge pull request #2755 from cowlicks/size-0-matrices

ENH: sparse: Size zero sparse matrices

This allows sparse matrices to be size zero. Meaning one or both of their shape
components is == 0.  This is needed because fancy indexing and other weird edge
cases occasionally need to return a size-0 sparse matrix.
Merge pull request #2745 from Juanlu001/sparse-docs

DOC: sparse: Added basic sparse matrix vector product examples
ENH sparse: Add Ellipsis indexing to sparse matrices.
TST sparse: Add tests Ellipsis indexing with sparse matrices.