how to solve this system differential equations?

2017 11 15
1 个回答
更新时间：2017 11 15
5 次查看（30 天）

0 个投票

0.2x(t)^9(1/x(t)^8+1/y(t)^8-1)^(9/8)=diff(x(t),t)
0.3y(t)^9(1/x(t)^8+1/y(t)^8-1)^(9/8)=diff(y(t),t)
with x(0)=1 and y(0)=1

请先登录，再回答此问题。

回答（1 个）

Torsten
Torsten 2017-11-15
编辑：Torsten 2017-11-15

0 个投票

fun=@(t,z)[0.2*z(1)^9*(1/z(1)^8+1/z(2)^8-1)^1.125;0.3*z(2)^9*(1/z(1)^8+1/z(2)^8-1)^1.125];
tspan = linspace(0,5,20);
z0 = [1;1];
[T,Z] = ode15s(fun,tspan,z0);
plot(T,Z(:,1),'-o',t,Z(:,2),'-o')
Best wishes
Torsten.

2 个评论
显示 无 隐藏 无

thanos zisopoulos
thanos zisopoulos 2017-11-15
thanks..in the last line you have t replace with T??
Torsten
Torsten 2017-11-15
编辑：Torsten 2017-11-15
Yes, becasue they have different meanings.
t is a formal parameter in a function call, T is the vector of output times at which the solution Z is supplied.
Best wishes
Torsten.

请先登录，再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Mathematics 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by