solvers.m

% solvers for DEs in ode.m

% TODO Insert relevant values into `range` and `initial_conditions`
% NOTE initial_conditions is a vector (ordered collection) of input values that will
%      unpacked in ode.m file into the X, L, S, N, and H variables respectively.
range = [1:50]
initial_conditions = [600, 3, 0, 3, 3]

# Use ode45 solver
[time_solution, var_solution] = ode45(@ode_sys, range, initial_conditions)

% Plot solutions, if necessary