subploting a simple ode45 predation model with global parametersb in matlab

main.m
% Define global parameters
global r_a r_b a b

% Initialize parameter values
r_a = 1.5;
r_b = 0.5;
a = 0.6;
b = 0.3;

% Define the ODE system
f1 = @(t,x) r_a*x(1)*(1 - x(1)/a) - x(1)*x(2)/(1 + x(1)^2);
f2 = @(t,x) r_b*x(2)*(1 - x(2)/b) + x(1)*x(2)/(1 + x(1)^2);

% Define the time interval
tspan = [0 20];

% Define the initial conditions
x0 = [0.5; 0.5];

% Solve the ODE system using ODE45 and plot the results
[t,x] = ode45(f1, tspan, x0);
subplot(2,1,1)
plot(t, x(:,1), 'b', 'LineWidth', 1.5)
xlabel('Time')
ylabel('Prey Population')
hold on

[t,x] = ode45(f2, tspan, x0);
subplot(2,1,2)
plot(t, x(:,2), 'r', 'LineWidth', 1.5)
xlabel('Time')
ylabel('Predator Population')
676 chars
33 lines