To solve this equation numerically using MATLAB/OCTAVE, I typically discretize the domain into small segments and approximate the derivatives using finite difference methods through codes. Plotted the temperature distribution along the rod. Analyzed the results as needed.