Okay, a little update, thank you for all your support. I've decided to neglect convection as btrueblood suggested, because it doesn't affect temperatures at all. Also the reason behind using constant step is that I can adjust numerical result by using least squares method for discrete data and...