how to solve this system differential equations?

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

Antworten (1)

Torsten
Torsten am 15 Nov. 2017
Bearbeitet: Torsten am 15 Nov. 2017

0 Stimmen

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 Kommentare

thanks..in the last line you have t replace with T??
Torsten
Torsten am 15 Nov. 2017
Bearbeitet: Torsten am 15 Nov. 2017
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.

Melden Sie sich an, um zu kommentieren.

Kategorien

Mehr zu Mathematics finden Sie in Hilfe-Center und File Exchange

Gefragt:

am 15 Nov. 2017

Bearbeitet:

am 15 Nov. 2017

Community Treasure Hunt

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

Start Hunting!

Translated by