Coupled ODE Solver Producing Unexpected Graph

1 visualización (últimos 30 días)
Brooke Mitchell
Brooke Mitchell el 2 de Feb. de 2020
Editada: darova el 30 de En. de 2021
Incorrect graph generating with given code. Trying to solve a couped ODE problem. Thank you in advance.
clear;
clc;
t_initial=0;
t_final=800;
x0=500;
y0=200;
tSpan=[t_initial t_final];
v0=[x0 y0];
[t,v]=ode45(@ODEfun,tSpan,v0);
figure
plot(t,v(:,1),t,v(:,2));
legend('Fox','Rabbit');
xlabel('Time (days)');
ylabel('Population');
title('Rabbit and Fox Population Overtime');
function f=ODEfun(~,v)
x=v(1);
y=v(2);
k1=0.02;
k2=0.0004;
k3=0.00004;
k4=0.04;
dx_dt=(k1*x)-(k2*y);
dy_dt=(k3*x*y)-(k4*y);
f=[dx_dt;dy_dt];
end
  2 comentarios
darova
darova el 2 de Feb. de 2020
Here is what i got for t_final = 80
img1.png
VBBV
VBBV el 29 de Nov. de 2020
I think you seem to get correct graph only for the given input conditions
Do you mean a graph like this below ?

Iniciar sesión para comentar.

Respuestas (1)

Pravin Jagtap
Pravin Jagtap el 6 de Feb. de 2020
Hello Brooke,
Please refer to following documentation for solving coupled system of ODEs

Categorías

Más información sobre Programming en Help Center y File Exchange.

Etiquetas

Productos


Versión

R2018a

Community Treasure Hunt

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

Start Hunting!

Translated by