1548.sagews

def logistic(N0, K, r, t):
    return K/(1.+((K-N0)/N0*exp(-r*t)))
    
plot(logistic(1000., 500., 0.1, x),(0.,100.)) + plot(logistic(100., 500., 0.1, x),(0.,100.))