%%% Clean up workspace 
clear vars 
 
%%%variables 
tspan = [0:0.01:1];
K = 10; 
y0 = 2;
r = 2; 
 
%%% Population growth ODE solution 
[t,y] = ode45(@(t,y) y0*K/(y0+(K-y0)*exp(-r*t)), tspan, y0);
 
%%% Plot of solution 
plot(t,y,'-o')