含边值问题的分段微分方程如何求解

肉肉坨

问题描述





当-1<=x<0时
y1’=y2;
y2’=-0.7*y1+-1.9*y1^3+a3—2.61*y3^5+z1;
z1’=z2;
z2’=y1;

当0<x<1时
y1’=y2;
y2’=0.0244*y1+-0.0063*y1^3+0.0065*y3^5+z1;
z1’=z2;
z2’=y1;

边界条件：
y2(-1)=0,y2(1)=0;
z1(-1)=0,z1(1)=0;

并且在x=0上满足：
y1(0-)=y1(0+);y2((0-)=y2(0+);
z1(0-)=z1(0+); z2(0-)=y2(0+);

求y1,y2,z1,z2

看着像一个边值问题的微分方程，且是分段函数，但是在分段的界限上满足连续的条件。

各位大侠有什么好点子吗？

[ 本帖最后由 eight 于 2008-1-22 20:45 编辑 ]

xjzuo

试试bvp4c, 分段的话，自己设置判断语句即可。