function dydt = myode1(t,y)
dydt = 48*y-y.*y; % Evalute ODE at time t