Please help with only the MATLAB question!!!!!

Consider an lx w cm fin, with d mm thickness. Assume that power Wis supplied along the [a,b] interval on the left edge. According to Fourier's first law heat energy per unit time, q, is where K is the thermal conductivity constant of the material; A is the cross sectional area of the material, Vu is the gradient of the temperature. Newton's law of cooling defines q in terms of convection: where H is the convective heat transfer constant of the material:; up is the ambient temperature Based on that formulate the 2D discrete mathematical model, based on the ecliptic PDE with Robin (convective) boundary conditions. 6 marks] Write and test MatLAB code which solve the problem and plots the results as a 2D contour plot. [4 marks] Consider an lx w cm fin, with d mm thickness. Assume that power Wis supplied along the [a,b] interval on the left edge. According to Fourier's first law heat energy per unit time, q, is where K is the thermal conductivity constant of the material; A is the cross sectional area of the material, Vu is the gradient of the temperature. Newton's law of cooling defines q in terms of convection: where H is the convective heat transfer constant of the material:; up is the ambient temperature Based on that formulate the 2D discrete mathematical model, based on the ecliptic PDE with Robin (convective) boundary conditions. 6 marks] Write and test MatLAB code which solve the problem and plots the results as a 2D contour plot. [4 marks]