and r is the residual vector of m settlement markers in each ground
array. The Levenberg–Marquardt (LM) method [9] uses a search direction
pk to iteratively minimize F and to find an optimum estimation of
the unknown parameters (i.e. s1, a1, s2 and a2) in Eq. (2). pk is a solution
of the linear set of equations as follows: