subploting a simple ode45 predation model with global parametersin matlab in matlab

main.m
% Define global parameters
global r1 r2 a b

r1 = 0.5; % Growth rate of species 1
r2 = 0.3; % Growth rate of species 2
a = 0.2;  % Predation rate
b = 0.5;  % Conversion efficiency

% Define the system of ODEs
dydt = @(t,y) [r1*y(1)*(1 - y(1)/K1) - a*y(1)*y(2); 
               r2*y(2)*(1 - y(2)/K2) + b*a*y(1)*y(2)];

% Define initial conditions and time span
y0 = [50; 100]; % Initial population sizes
tspan = [0 100]; % Time span of simulation

% Solve the system of ODEs using ode45
[t,y] = ode45(dydt, tspan, y0);

% Plot the population dynamics of the two species over time
subplot(2,1,1)
plot(t, y(:,1), 'r', 'LineWidth', 2)
xlabel('Time')
ylabel('Population size')
title('Population dynamics of species 1')

subplot(2,1,2)
plot(t, y(:,2), 'b', 'LineWidth', 2)
xlabel('Time')
ylabel('Population size')
title('Population dynamics of species 2')
850 chars
32 lines