I'm trying to solve the following FEM equation for 1D unsteady heat conduction using MATLAB. I'm new to FEM and to the mathematical approach in MATLAB. So, I'd appreciate any help/suggestion.
What I need help with is to get the converged solution.
At i=0, T=10 and at i=5, T=50. When n=0, T_i^(n+1)=10, 73, 100, 85, 50 respectively at i=1, 2, 3, 4, 5.
dt=0.1 and (alpha/dx)=2
Again, I'd appreciate any and every suggestion here.
Thank you in advance!