Add Gauss-Legendre and Gauss-Laguerre quadratures

Added gauss_legendre() and gauss_laguerre() functions and tests for calculating
the quadrature points and weights. The functions work for any order and any
precision.
Add Gaussian integrals into docs
Fix docstring formatting, added "Returns" section

Now it should follow SymPy's conventions. Also it should make clear that the
function returns numerical values (not symbolic).
Simplify condition to "a < root < b"
Use `` instead of '
Fix bisect's docstring
Fix RootOf.n() for complex roots

The problem was the same as for the real roots that was fixed in
5c58406b3be8fc6d75d0e258a67b6847d5274a4f. This commit provides the same fix for
complex roots by refining the complex interval.
Add a test for RootOf() for complex roots

These tests are correct, but currently fail because they have revealed a few
other bugs, which will be fixed in subsequent commits.  After the other bugs
are fixed, these tests will start passing on their own.
Do not store the refined interval in RootOf

This fixes a caching bug. A test has been added, that failed before.
What happened is that the evalf() method refines the interval and stores it
back into the RootOf() method, which uses global cache... Subsequent calls to
RootOf with the same polynomial and the *same* root index then assigns a wrong
interval for the root.

Ideally, the RootOf() implementation should be improved so that it is robust
against refining the intervals. For the time being, this patch fixes the issue.
Refine the complex interval properly

Now all tests pass. Previously a wrong root was returned due to numerical
inaccuracies.

In particular, the following case was happening: the imaginary part of the root
was supposed to be in [0, 3] and it was a complex root, so it couldn't be 0.
However, due to numerics, the returned imaginary part was 1e-69. Which belongs
into this interval technically, but in reality the returned root was simply a
real root, already returned with some other root index. (The real part happened
to fit into the same interval by accident.)

This patch makes sure that the imaginar part is always distinct from other
roots. That way things start working.
add global_assumptions import to docs
replace global_assumptions with assuming in docs
Travis-CI: add tests for Python 2.7/3.3 + gmpy

This adds two new tests:

Python 2.7 + gmpy
Python 3.3 + gmpy

At the moment, there are some test failures. So the tests are marked as
"Allowed Failures". Travis-CI then tests them, but keeps the results separate
from the main build. That way it will not disturb the SymPy workflow, but we
can work on fixing the failures, as well as manually examining the results to
make sure that no more failures are being introduced.
Make RootOf() work with gmpy

Now it works with gmpy as well as Python ground types.
.travis.yml: force pip to install gmpy 1.16
Merge pull request #1756 from certik/travis_gmpy

Travis-CI: add tests for Python 2.7/3.3 + gmpy
Travis-CI: propagate tests failures back to Travis

The problem was that somehow the nonzero exit code was not always propagated to
Travis. I have found that using just one executable statement fixes it. So now
we run both tests and doctests using one Python session.
Merge pull request #1755 from mrocklin/assuming-sphinx

Sphinx docs for assuming
Merge pull request #1758 from certik/travis_fix

Travis-CI: propagate tests failures back to Travis
Merge pull request #1736 from prasoon2211/evalf-trig-iss3538

Passing a dictionary to .n() now works on cos and sin. Issue 3538