function y = f_ode(x,a,K)

y = a*(K*x - x.^2);