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

Answers (1)

Torsten
Torsten on 15 Nov 2017
Edited: Torsten on 15 Nov 2017
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 Comments

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

Sign in to comment.

Categories

Find more on Mathematics in Help Center and File Exchange

Asked:

on 15 Nov 2017

Edited:

on 15 Nov 2017

Community Treasure Hunt

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

Start Hunting!