% Simulate a Lotka-Volterra system

tspan = 0:0.05:12;  % time points
x0 = [1; 0.5];           % initial condition
k   = [3; 1.5];           % rate constants

% Define right-hand side of the ODE, i.e. f(x,k)
odefun = @(t,x) [k(1)*x(1) - k(2)*x(1)*x(2);
                            k(2)*x(1)*x(2) - x(2)];
                          
 % Solve the system
??

% Plot the solution
???