mathematical models of several applications in Engineering, Economy, Biology and different
branches of Physics. These kinds of problems present several analytical and numerical
difficulties related, for example, to time evolution and moving boundary. In this work we
implement a numerical method based on the finite difference scheme for time evolution and
nonlinear complementarity algorithm (FDA-NCP) for solving the problem at each time step …