pleas solve it useing matlab or ess

Consider a two-dimensional heat conduction in a square cross section (4cm4cm,k=40W/mC, =6.694106m2/s ) as shown in the figure. The entire block is initially at 100C. The right and left surfaces of the block are maintained at constant prescribed temperatures of 300C and 100C, respectively. The bottom surface is exposed to a constant heat flux of 1000W/m2. The top surface is in contact with a convective environment (h=45W/m2K) at 20C. Using a uniform mesh size of x =y=1cm. Obtain the finite difference formulation of this problem for transient two-dimensional heat transfer, and using a reasonable time step, determine the temperature of each node with respect to time and plot them in one diagram. a) Determine the time needed to reach steady state conditions and determine the steady state temperatures. b) Obtain the finite difference formulation of this problem for steady two-dimensional heat transfer, determine the unknown nodal temperatures by solving those equations and compare them to the temperatures determined from part (b), and c) Determine the rate of heat transfer from the right surface by convection in the steady state conditions