1、存在问题：微分方程函数和边界条件函数的定义，function后面没有空格，导致两个函数被误作为变量，根本没起到作用。之所以没有报错，是因为twoode和twobc作为系统提供的例子，确实有这两个函数。这也是从系统示例修改程序时应该尽量避免的做法。 2、不知道你的初始估计solinit是怎么生成的，可以考虑使用下面的做法：solinit = bvpinit(linspace(0,4,5),[0 0 0 0]);
