Improved discover_equation function by avoiding using the known
valuation of qE.

Fixed also an example from README.rst.