How to numerically simulate the robust topological edge state animation in Matlab in a square lattice?

Question

Topological physics has enabled intriguing possibilities in photonics by hosting disorder-immune robust edge state propagation. I want to know how can we numerically simulate th