一维稳态导热问题数值模拟