This question is to be coded in Matlab

You will develop a 2D finite difference model to simulate the propagation of an acoustic wave. The model is to contain two layers that are separated by a reflecting / transmitting interface. The model specifications are given below.

1.) Test Cases Build suite of test cases that will be used to benchmark your code. These are to be based on analytical 'textbook' solutions. You should develop at least two test cases and one should test the behaviour of the wave-field at the internal interface. Show that the test cases are appropriate for a 2D model of acoustic wave propagation. 2.) Pseudo Code The pseudo code should contain all major elements of your code. Loops or other iterations over matrix elements should be detailed. 3.) Finite Difference Model Minimum Model Size 100 100 Time> 1000 steps Wavelet Source at one grid point (i.e. no entire boundary) Acoustic wave equation to be solved with time-marching five-point centred grid scheme. (More sophisticated schemes like nine-point centred grid or staggered-grid are permitted) One internal interface allowing reflection and transmission Absorbing boundaries (this will be considered an advanced feature) * .) Visualisation Synthetic seismogram with shot gather receiver/source configuration Snapshots of the wave-field at pertinent time steps * * 1.) Test Cases Build suite of test cases that will be used to benchmark your code. These are to be based on analytical 'textbook' solutions. You should develop at least two test cases and one should test the behaviour of the wave-field at the internal interface. Show that the test cases are appropriate for a 2D model of acoustic wave propagation. 2.) Pseudo Code The pseudo code should contain all major elements of your code. Loops or other iterations over matrix elements should be detailed. 3.) Finite Difference Model Minimum Model Size 100 100 Time> 1000 steps Wavelet Source at one grid point (i.e. no entire boundary) Acoustic wave equation to be solved with time-marching five-point centred grid scheme. (More sophisticated schemes like nine-point centred grid or staggered-grid are permitted) One internal interface allowing reflection and transmission Absorbing boundaries (this will be considered an advanced feature) * .) Visualisation Synthetic seismogram with shot gather receiver/source configuration Snapshots of the wave-field at pertinent time steps * *