@xzczd 非常感谢啊 麻烦你挺晚了还看 按照基本概念你说的对啊 我试了一下 应该是没问题的 而且你这个做法对本征值为实数保持得更好 能看下你怎么写的吗 我自己写的得让eps取个非零小量来避免左边界点除0

不过说实在的 我怎么觉得我用的那个边条件 按照ODE在0点的渐近分析貌似没错啊 （不过肯定是哪里不对……） 还是说只要条件overdertermined就会坏掉…

另外这个问题不知道为什么 不论是naive FDM还是你的那种 解画出来经常遇到这种两支不同曲线混到一起的情况 你以前有遇到过类似的吗

@xiaohuamiao 我就是在你那个基础上改了一点点，主要就是#[[2 ;; -2]] &改#[[1 ;; -2]] &，再把下面构造阵的部分的bcL去掉。eps我直接取10^-3了，更正统的做法我没仔细考虑……

关于混淆的问题，我一时想不起来实例，不过会不会是变量排序不对？

还有我也觉得你的渐近分析没问题，到底哪里不对一时间搞不清楚……

嗯 我也是这么做的 不过应该不是变量排序 我反复验证过 俩code排序不同 而且pdetoae写的那个是最简单的排序 没有混合alpha beta那些分量 但是俩code都出现了…

@xzczd 能解决还是很感谢你的👍

至于混合的问题 我也反复查过 特征值还是对的 因为同一特征值解不唯一 怀疑是把不同的解混合了 比如有的明显像是奇数点一个解偶数点另一个解 虽然是不对的 但可以直观想象是两个解交错插到一起mesh扩大一倍 或许在FDM意义下点够密就对特征值并没影响…