>>> from datetime import date
>>> x = (date(2014, 2, 22) - date(2009, 1, 3)).days
>>> x
1876

result <- lm(log10(prices) ~ months)

def iter_price(initial_price, steps, base=2):
    for i in steps:
        yield initial_price * base ** i

def iter_mbtc_left(initial_btc, rake, steps):
    base = (1 - rake)
    for i in steps:
        yield initial_btc * base ** i

def iter_mbtc_sold(initial_btc, rake, steps):
    piece = initial_btc * rake
    base = (1 - rake)
    for i in steps:
        yield piece * base ** (i - 1)

def iter_kusd_out(mbtc_sold, price):
    for mbtc, usd in zip(mbtc_sold, price):
        yield mbtc * usd / 1000.

def iter_kusd_sum(kusd_out):
    s = 0
    for kusd in kusd_out:
        s += kusd
        yield s

def iter_mbtc_val(mbtc_left, prices):
    for mbtc, price in zip(mbtc_left, prices):
        yield mbtc * price / 1000.

def iter_total_val(mbtc_val, kusd_sum):
    for item in zip(mbtc_val, kusd_sum):
        yield sum(item)

def main():
    base = 2 # doubling
    numof_steps = 10
    initial_price = 1 # kUSD
    inital_btc_stash = 10000 # mBTC
    rake = 0.1

    steps = range(1, numof_steps + 1)
    prices = list(iter_price(initial_price, steps, base=base))
    mbtc_left = list(iter_mbtc_left(inital_btc_stash, rake, steps))
    mbtc_sold = list(iter_mbtc_sold(inital_btc_stash, rake, steps))
    kusd_out = list(iter_kusd_out(mbtc_sold, prices))
    mbtc_val = list(iter_mbtc_val(mbtc_left, prices))
    kusd_sum = list(iter_kusd_sum(kusd_out))
    total_val = list(iter_total_val(mbtc_val, kusd_sum))

    print ' Price  mB left  mB sold  k$ out  mB val.  k$ sum  total val.'
    fstr = ' {0:4} {1:6.0f} {2:8.0f} {3:7.0f} {4:8.0f} {5:7.0f} {6:8.0f}'
    for args in zip(prices, mbtc_left, mbtc_sold,
                    kusd_out, mbtc_val, kusd_sum, total_val):
        print fstr.format(*args)

if __name__ == '__main__':
    main()

prices <- c(
      0.005,   0.005,   0.005,   0.005,   0.005,   0.005,
      0.005,   0.005,   0.005,   0.005,   0.005,   0.005,
      0.005,   0.005,   0.005,   0.005,   0.005,   0.008,
      0.06,    0.065,   0.062,   0.106,   0.27,    0.24,
      0.39,    0.90,    0.85,    1.50,    6.38,   18.55,
     14.10,    9.75,    5.76,    3.30,    2.60,    3.51,
      6.11,    5.06,    4.88,    4.98,    5.06,    6.03,
      8.04,   10.88,   11.46,   11.58,   11.51,   13.33,
     16.31,   25.93,   58.14,  115.09,  112.81,  108.39,
     83.80,  105.80,  124.33,  156.50
)

months <- 1:length(prices)
result <- lm(log10(prices) ~ months)

plot(months, log10(prices))
abline(result, col="red")

a <- coef(result)[1]
b <- coef(result)[2]
r2 <- summary(result)$r.squared