In this article we introduce a numerical algorithm for finding harmonic mappings by using the shear construction introduced by Clunie and Sheil-Small in 1984. The MATLAB implementation of the algorithm is based on the numerical conformal mapping package, the Schwarz-Christoffel toolbox, by T. Driscoll. Several numerical examples are given. In addition, we discuss briefly the minimal surfaces associated with harmonic mappings and give a numerical example of minimal surfaces.