热耦合方腔流动，求解带有温度场的平板驱动方腔流动

代码说明：

本程序是用来求解平板驱动方腔流动的，利用涡量——流函数方法，耦合温度场，泊松方程采用迭代求解，空间导数采用中心差分。