Finite Difference Method to find Steady State

Hello,

This relates to 1-D unsteady heat conduction using Finite Difference Method (FDM).

I have to solve the following equation, by iterating until the steady state aproaches.

Steady state is basically;

Following data is known. (calculated as an example)

When n=0, T_i^(n+1)=10, 73, 100, 85, 50 respectively at i=1, 2, 3, 4, 5.

When n=1, T_i^(n+1)=10, 62.2, 87.4, 79, 50 respectively at i=1, 2, 3, 4, 5.

How should I approach this? I'd appreciate your recommendations, reading materials, suggestions and/or pointers rather than a complete code at this point.

I know the question is vague but I'd welcome any suggestions for starting points :)

Reference: http://faculty.washington.edu/finlayso/ebook/pde/FD/FDcode.htm (I'm assuming the code should be something similar to this)

Previously posted question: https://www.mathworks.com/matlabcentral/answers/516366-finite-difference-method-to-find-steady-state