What part of it troubles you? It is a simple system of PDEs with an explicit scheme to solve. You can introduce a 2D square grid of points, introduce the initial condition on the grid points, and write a for loop for the last two equations. Then save the results of each step for all points and then you will get your space-time solution.