Partial Fractions Expansion
for Inverse Laplace Transform


GNU Octave
(Sistema em código aberto compatível com Matlab)
I really like the one-dimensional approach
used in the Wikipedia article about Newton's
method. Then, seeing the Jacobian and thee Hessian
as N-dimensional generalizations of the first
and second derivatives will really make things
easier to swallow.
Section 2.2 "Overview of Algorithms" from
Nocedal & Wright's book provides a smooth
introduction about what quasi-Newton methods
do. Then, chapter 6 dives into the specifics,
but honestly I still didn't take a look at.
Nocedal, J., Wright, S. (2006). Numerical Optimization.
United States: Springer New York.
Some colleagues have been using SciPy's
implementation of L-BFGS-B quite succesfully
for large scale problems (e.g. FWI).
https://docs.scipy.org/doc/scipy/reference/optimize.minimize-lbfgsb.html
Discretization of ODE- and
PDE-based systems
For an undergraduate-level introduction
to discretization, see pages 256-269 from Green, R. A., Lathi, B. P. (2018). Linear Systems and Signals. United Kingdom: Oxford University Press.
Although this video is about solving PDEs,
it provides a nice hint about how they can
be discretized.