Could I have help writing question 1 into matlab code?

Consider the 2D Steady Heat Equation: g(x, y) on a 1 x 1 square. The dependent variable T represents the equilibrium temperature and q(x) is a volumetric heat source. Assume the boundaries of the square are at a fixed temperature of 300K (room temperature) Assume there is uniform electrical heat generation in the domain which results in q(x,y) 500 W/m. The thermal conductivity of the material is k 1 W/(mK). . Using the provided matlab code, calculate the heat flux through the bottom of the domain using a mesh of 64 64 segments. . Repeat the process using COMSOL and output the heat flux for a mesh of the same resolution . Calculate the percentage difference in the two results. . Repeat the process using a mesh of 128 x 128 segments