Home

Illustrations

A few examples of what to expect with periodic boundary conditions, 𝑑𝑚 = 0, 𝑑𝑥 = 0, 𝑏𝑚 = 0, 𝑏𝑥 = 1, 𝑠𝑚 = 1, 𝑠𝑥 = 1, 𝑙𝑥 = 2.5 (and ℎ = 16) for the constant rates, and t_𝑖𝑛𝑖𝑡 = 0, t_𝑒𝑛𝑑 = 1000000, 𝑑𝑡 = 0.01, 𝑡 = np.arange(t_init, t_end, dt), 𝑁𝑡 = len(t) for the time parameters.

Figure 1: basic simulation of a single layer of skin (matrix L) with an absence of chromatophores (white in the figure, 0 in the code), xanthophores (yellow, 1) and melanophores (black, 2). This layer of epidermis can be cut (red, 3)

Figure 2: simulation with a 2 pixel wide band of iridophores, located in an LI matrix containing either an absence of chromatophore (white, 0), or iridophores (blue, 1). This layer of epidermis can also be cut (red, 2). The LI matrix is visible by transparency of the matrix L (α=0.7), which explains why the blue is grayed out

Figure 3: simulation with a skin defect, visible in the lower quarter, 20 pixels wide, on a single layer of skin (L matrix)

Figure 4: simulation with a band of 2 pixels wide of iridophores as well as a skin defect, visible in the lower quarter, 20 pixels wide

Figure 5: simulation with a 2 pixel wide band of iridophores as well as a diagonal cut of width 15 pixels. Healing rates set to 1 for all cell types: 𝒉𝒙=1, 𝒉𝒎=1, 𝒉𝒊=1 and 𝒉𝒔=1