def alpha_for_N(comps,Influx=10,Outflux=1000):
    """calculates alpha based on required outflux for given influx"""
    c=(Outflux/Influx)**(1/float(comps))
    alpha_sol= (c-2)/(c-1)
    return alpha_sol