packie 发表于 2006-10-12 17:44

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

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

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

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
中设计指数的全部用.^,乘法用.*

packie 发表于 2006-10-12 18:23

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

toes 发表于 2006-10-12 18:30

原帖由 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个初值，就是这个：
=ode45('ode1',,,[],b3(i));

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

packie 发表于 2006-10-12 20:08

能不能帮助修改一下，谢谢
页: [1]
查看完整版本: 求助常微分方程组求解问题