FDM for the time dependent wave equation in Julia

I am writing a test code of our new time-dependent wave equation in Julia.

For now, it is equivalent to the normal wave equation and it looks fine.

I am not sure Julia is that fast and easy to write, but one of the merits is that the visualization functions can be called within Julia's framework.