%popul 
function dxdt=populatio(t,x, r, k)
dxdt=r*x*(1-(x/k));
end