Generate the initial conditions matrix for a finite difference

Question

I want to make a initial condition matrix for this finite difference relations:

w_{i+1,j} - 4w_{i,j} + w_{i-1,j} + w_{i,j+1} + w_{i,j-1} = h^2 g(x_i, y_j)

with the