Exercise3

A typical problem in thermodynamics is to find the steady-state temperature distribution inside a thin plate if we know the temperature around the boundary. Let \(T_1, T_2, \ldots, T_6\) be the temperatures at the six nodes inside the plate as shown below.

<<SVG image is unavailable, or your browser cannot render it>>

The temperature at a node is approximately the average of the four nearest nodes: for instance,

\begin{equation*} T_1 = (10 + 15 + T_2 + T_4)/4 \text{,} \end{equation*}

which we may rewrite as

\begin{equation*} 4T_1 - T_2 - T_4 = 25 \text{.} \end{equation*}

Set up a linear system to find the temperature at these six points inside the plate. Then use Sage to solve the linear system. If all the entries of the matrix are integers, Sage will compute the reduced row echelon form using rational numbers. To view a decimal approximation of the results, you may use

A.rref().numerical_approx(digits=4)
      

In the real world, the approximation becomes better the closer the points are together, which happens as we add more and more points into the grid. This is a situtation where we would want to solve linear systems having millions of equations and millions of unknowns.

in-context