接下来，给出一个采用单松弛因子模型(LBGK)的MATLAB通道流程序。该程序包括一个主程序(Poiseuille.m)和两个子程序(Stream.m和FEQf.m)，子程序分别负责迁移计算和平衡态分布函数的计算。
1.主程序Poiseuille.m
2.子程序Stream.m
3.子程序FEQ f.m
