We perform an in silico investigation of the formation of multiple intense zebra stripes by extending the domain with an appropriate extending speed. The common zebra has alternating dark and light stripes, creating a two phase pattern. However, some Equus burchelli zebras have an intermediate gray color stripe situated between the dark and light stripes. To numerically investigate the formation of multiple intense zebra stripes, we first find the equilibrium state of the governing system in the one-dimensional (1D) static domains using various frequency modes. After finding the equilibrium state for the governing system in the 1D static domains, we stack numerical data. Then, we load the stacked numerical data to use as an initial state for finding the growth rate that forms the multiple intense zebra stripe formation in the 1D extended domains. Next, convergence experiments are conducted to verify the convergence of the numerical method for the governing system. Finally, numerical simulations are performed to confirm the formation of multiple intense zebra stripes in two-dimensional extending domains and on evolving curved surfaces.