I want to solve a 1-D heat conduction PDE using Matlab which looks like $$\rho c_p \dfrac{\partial T}{\partial t} = \dfrac{\partial}{\partial z}\left( \lambda \dfrac{\partial T}{\partial z} \right),$$where the initial condition is T1 (in Kelvin) and the final condition is T2 (also in Kelvin), the difference being T1 takes different values from t=0 hours to t=24 hours (t being time) and T2 being constant (say 200 K)I checked the pdepe toolbox in Matlab and I felt the conditions that I wanted to implement might not be possible through the toolbox.Can anyone suggest a way to create a loop through which T1 can read values from an array so that a diurnal plot can be created ?
Solving a PDE using Matlab (with varying initial conditions)
matlab;pde