求助常微分方程组求解问题

packie

二阶常微分方程组求解问题，我编制了个程序，运行总是有错误，恳请高手解决，谢谢

[ 本帖最后由 lxq 于 2006-10-12 20:10 编辑 ]

zfx620

dy(2)=-a*(y(1))^(1./2)-a*(y(1))^(3./2);
dy(3)=y(4);
dy(4)=b1*a*(y(1))^(1/2)*y(2)-(b1/b3)*a*(y(3))^(1/2)*y(4)+b1*a*(y(1))^(3
中设计指数的全部用.^,乘法用.*

packie

不对啊
不是这个问题的
期待大家的帮助
谢谢

toes

原帖由 zfx620 于 2006-10-12 18:00 发表
dy(2)=-a*(y(1))^(1./2)-a*(y(1))^(3./2);
dy(3)=y(4);
dy(4)=b1*a*(y(1))^(1/2)*y(2)-(b1/b3)*a*(y(3))^(1/2)*y(4)+b1*a*(y(1))^(3
中设计指数的全部用.^,乘法用.*

主要不是这个问题。
阶次，乘号前的.应该都可以不要。

主要是你的主程序中只给了2个初值，就是这个：
[T,Y1]=ode45('ode1',[0,20],[0;10],[],b3(i));

y(3) y(4) 没赋初值。

packie

能不能帮助修改一下，谢谢