圣维南方程求解。您好，我自己编的程序为什么矩阵无解啊？怎样设置五对角矩阵的参数啊？

问题遇到的现象和发生背景

问题相关代码，请勿粘贴截图

运行结果及报错内容

我的解答思路和尝试过的方法

我想要达到的结果@赖亦无

就目前我的认知，圣维南方程在干涸河道的应用要考虑露滩的问题。
河道初始水位一般取河道断面监测到的水位即可。