Input , , , MaxIter, and . Initialize
For k = 1: MaxIter
 compute using the explicit scheme (18).
 update the parameter according to (19).
end for
Algorithm 1: The time-marching gradient descent algorithm for solving (13).