Write a finite difference scheme for $u=f(x,y).$ Write a MATLAB code for the numerical solution of the following $2D$ boundary value problem

$-u=12xyin[0,1][0,1]$

$u(0,y)=u(x,0)=u(1,y)=u(x,1)=0.$

Demonstrate MATLAB code by plotting the solution at the respective node points.