python fipy 求解偏微分方程

1、首先简介问题。一个二维平板，顶端1摄氏度（100也可以，只是一个系数），另外三个边缘0摄氏度，初始时刻整个板子都是0摄氏度，随之时间的推进，热量在板子上传递，最后达到平衡态，我们不仅希望知道平衡态的温度分布，也希望知道温度随时间是如何变化的。如果这篇经验过于复杂者不够直接或者fipy无法安装，可以参考另一篇python偏微分方程经验 jingyan.baidu.com/article/46650658d71872f548e5f852.html

3、先看一张最后的结果过个瘾，也比较容易想象，最后整个板子的温度分布大致呈现怎样

6、创建微分方程

8、创建画图

10、建议在命令行里面运行，命令行里面可以获得动图，ipython里面只有最后一张图

12、主要参考文章为fipy的沣举泉零官方示例 www.ctcms.nist.gov/fipy/examples/diffusion/爿讥旌护generated/examples.diffusion.mesh20x20.html www.ctcms.nist.gov/fipy/examples/diffusion/generated/examples.diffusion.mesh1D.html#equation-eq:diffusion:mesh1D:constantD