poisson方程的狄利克雷边值问题五点差分怎么用matlab编程呀？

就像图中一样写出了对应的矩阵形式，然后系数矩阵中的每个C D都表示一个矩阵，请问怎么编写程序呢？