同学,你直接贴代码吧,要找错肯定得试进行,这么一张图,那么长的代码,我重新输会输死人的啊。
顺便,在Mathematica里使用For循环是很不明智的(它的速度很慢),建议使用Table或者Do或者While来实现。
在Mma里使用大写字母打头的变量名也是不明智的,尤其是你这里还用了N,要知道 N 已经被用作系统函数名了,这种做法可能引发不可预知的错误。
再顺便,你的这种类似穷举的寻找初始条件的方法,恐怕难以避免不合理初始条件的出现,所以可能还需要在NDSolve里面方个EventLocator,具体请参看帮助。
温馨提示:内容为网友见解,仅供参考