function NPrime = odefun(t,N) r = .8378367715; K = 49293; NPrime = r*N*(1 - N/K);