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 ExchangeMathematics についてさらに検索

質問済み:

2017 年 11 月 15 日

編集済み:

2017 年 11 月 15 日

Community Treasure Hunt

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

Start Hunting!

Translated by