求助 这个超越方程的数值解的代码哪里出错了

程序思想：
1、首先根据洛夫波的频散函数，即令x5=0，它是关于波数和波速的非线性方程；
2、然后固定波数的数值，使波速从2596到3128取值，给定某一搜索步长，搜索有根区间，只要找到一个有根区间，就调用二分法的求根子函数，求得零点；
3、最后改变频率的数值重复2步，得到不同频率时候的零点值。

但是程序是有问题的，请大神帮我看看~~~~~