1d steady state heat conduction matlab code