A mathematical model that successfully reproduces cellular patterns observed in sensory body tissue has been developed based on intensities of adhesion molecules as input parameters. A new numerical method for obtaining the solution of the model that can correctly handle complex cellular rearrangements preserves individual cell volumes and cell connectivity.