Please write a matlab code for this question:

The concept of thermal resistance can be used to find the temperature distribution in a flat square plate shown in the following figure. The plates edges are insulated so that no heat can escape, except at two points where the edge temperature is heated to Ta and Tb respectively. Modeling the temperature distribution in the flat square plate by using the finite difference approach, results in the following set of linear equations.

Solve the set of equations for the case where Ta and Tb are 150 and 20⁰C, respectively.