采用有限差分的方法，利用matlab实现，采用逐次更新矩阵的形式求解。