upgrade solving for abs, powers, and real/imag variables

Also:

    use powdenest before factoring in _tsolve

    fix detection of difference of functions
add bivariate solving to solvers
edits after coverage check
modify preprocessing for lambert