Fixed discover_equation function. Now can discover equations over cubic fields.
Added progress bars.
Also improved integration performance.