Simulation of unit operations on Matlab using iterative differential equation solution methods, Monte Carlo methods and statistics