function test
clear all; close all; clc;

a=2;
b=3;
y0=1;
t=linspace(0,5,100);

[t,y]=ode45(@(t,y)myODE(t,y,a,b),t,y0);
figure(1); plot(t,y); hold on;

c1=1-b;
y=a*c1*exp(-a*t)./(b+c1*exp(-a*t)).^2;
figure(1); plot(t,y,'r'); legend('Numerische Lösung','Analytische Lösung');

end

function dy=myODE(t,y,a,b)
dy=a*(1-b*y)*y;
end