The algorithm credited to Richardson and Lucy has the following form:
Note that in this form of operator notation, an operation between two vectors
is done point-by-point. For example, the division by is
point-by-point, and the multiplication by
is point-by-point.