“…can be written asC[m, n + 1] = C[m, n] − δ t 2 mδ R Jf [m, n] + J f [m, n](79)whilst initial and boundary conditions (10)-(13) can be written asC[m, 0] = C0 , 0 ≤ m ≤ M ; (80) Jf [0, n] = 0, 0 ≤ n ≤ N ; (81) Jf [M, n] = − J0 , 0 ≤ n ≤ N ; (82) Tr [M, n] = 0, 0 ≤ n ≤ N ; (83) ũ[0, n] = 0, 0 ≤ n ≤ N ;(84)respectively. By applying finite central, forward or backward (depending if it is computed at the inner or boundary of the domain) difference derivatives, at a fixed time step n, the evolution equation (79) can be solved by implementingC[m, n + 1] = C[m, n] − 1, n] , 0 < m < M − 1; (86) C[m, n + 1] = C[m, n] + δ t − 1, n] , m = M − 1; (87) C[m, n + 1] = C[m, n] + δ t δ R 2 m J0 + J0 + Jf [m − 1, n] , m = M,(88)by using(26) or(66) for the computations of the hydrostatic stress in Jf . SinceT r (R, t) = 2F(R, t),(89)by a backwards recursive method it holdsTr [m − 1, n] = Tr [m, n] − 2δ R F[m, n], 0 ≤ m ≤ M,(90)taking into account condition (83), while for the hoop stress it holds Tθ [m, n] = Tr [m, n] + H[m, n], 0 ≤ m ≤ M,…”