Solve the problem with Matlab

PLIAlh = A.S1. Problem 3 Consider a 2D reservoir with 7 by 7 grid blocks. Based on some simplifying assumptions, you know that you can use the steady state equation below to solve for pressure in the reservoir. a2p alp ax2 ay2 = The equation above can be shown as below when finite difference approximation is used. 8l Li — 2Pi + 4. j÷, 2Pij + 0 Ox)2 (AY)2 Solve for this equation to get pressure distribution in the reservoir. Make a 2D pressure map (contours) by using MATLAB plotting features. Do not forget about color bars and units. Pressure at the boundaries is 1000 psia (orange grid blocks) and the well block pressure (shown with blue color) is set at 100sia. P