New LCurve class to estimate regularization parameters

The new class was added to `fatiando.inversion.regularization`. It behaves like any other solver but under the hood it runs the L-curve and finds the corner value when `fit()` is called. 

The L-curve is built using user-specified values of the regularization parameter. The corner value is found using the "triangle method" of [Longina Castellanos et al (2002)](http://www.sciencedirect.com/science/article/pii/S0168927401001799).

[This gist](http://nbviewer.ipython.org/gist/leouieda/9099363) was the first prototype developed.

Notes:

* So far it's not behaving well for total variation regularization

TODO:

- [ ] Use it in recipes
- [ ] Fill docstrings