%% We define the parameter and their values
n = 1;
K = 0.1;
alpha = 1;
beta = 1;

% Other possible parameter values we may want to use:
% n = 2;
% K = 1;
% alpha = .1;
% beta = 1000;

%% We define a function that returns the derivative of X
% with respect to time dX/dt
dXdt = @(t, X) beta/(1+(X/K)^n) - alpha*X;

% We call Matlab's ODE solver, to solve the dynamics of X 
% between time = 0 and time = 30, starting from X=0 at time 0.
[t,X] = ode45(dXdt, [0 30], 0);

% We plot the dynamics and label the axes
plot(t, X, '-');
xlabel('time');
ylabel('X');
