mathematica 引用错误提示的结果
比如我解方程sol=NDSolve[y1'[x]==f1[x,y1,y2],(1-r[y1,y2])y2'[x]==f2[x,y1,y2],y1[0]==y10,y2[0...
比如我解方程
sol = NDSolve[y1'[x] == f1[x,y1,y2], (1-r[y1,y2]) y2'[x] == f2[x,y1,y2], y1[0]==y10, y2[0]==y20, {y1,y2}, {x,xmin,xmax}];
其中y10,y20,f1,f2,r已知,且f1,f2含参数a(可以给定)。这个方程组有奇点 r==1,解会有间断,我需要知道间断位置。
如果我就这么解,就会抛出错误提示:
NDSolve::ndsz: "At x == x_jump`, step size is effectively zero; singularity or stiff system suspected.
这里x_jump就是间断位置,和方程参数有关。(也试过)
我想画出x_jump关于a的变化规律,所以想引用这个错误提示里的x_jump。
有没有大佬能给个解决方案? 展开
sol = NDSolve[y1'[x] == f1[x,y1,y2], (1-r[y1,y2]) y2'[x] == f2[x,y1,y2], y1[0]==y10, y2[0]==y20, {y1,y2}, {x,xmin,xmax}];
其中y10,y20,f1,f2,r已知,且f1,f2含参数a(可以给定)。这个方程组有奇点 r==1,解会有间断,我需要知道间断位置。
如果我就这么解,就会抛出错误提示:
NDSolve::ndsz: "At x == x_jump`, step size is effectively zero; singularity or stiff system suspected.
这里x_jump就是间断位置,和方程参数有关。(也试过)
我想画出x_jump关于a的变化规律,所以想引用这个错误提示里的x_jump。
有没有大佬能给个解决方案? 展开
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询