ma2rten

Permit Number,Facility Name,Combined Owner,Regulated Operation,Permit Type,Regulated Activity,Allowable Count,Number Of Lagoons,Issued Date,Effective Date,Expiration Date,Admin Region,County Name,Location Lat Num,Location Long Num,Address 1,Address 2,City,State,Zip
AWC010002,Piedmont L/S Co Farm,Joe Jones,Cattle,Cattle State COC,Cattle - Beef Feeder,300,0,10-01-2014,10-01-2014,09-30-2019,Winston-Salem,Alamance,36.1742,-79.4994,,,,,
AWC010002,Piedmont L/S Co Farm,Joe Jones,Swine,Cattle State COC,Swine - Farrow to Wean,300,0,10-01-2014,10-01-2014,09-30-2019,Winston-Salem,Alamance,36.1742,-79.4994,,,,,
AWC010006,Covington Dairy Farm Inc,William Covington,Cattle,Cattle State COC,Cattle - Milk Cow,300,1,10-01-2014,10-01-2014,09-30-2019,Winston-Salem,Alamance,36.0442,-79.3261,3008 S Nc119,,Mebane,NC,27302
AWC010010,Triple W Farms,Harold Woody,Cattle,Cattle State COC,Cattle - Dairy Heifer,200,1,03-11-2016,03-11-2016,09-30-2019,Winston-Salem,Alamance,35.9042,-79.3,3545 E Greensboro-Chapel Hill Hwy,,Snow Camp,NC,27349

ma2rten

from __future__ import division

import sqlite3

db = sqlite3.connect('db')

def p_ben(n, history_length):
    p_remove = 1/3 # .3
    p_keep = 1 - p_remove

ma2rten

from __future__ import division
from math import log

def p_ben(n, history_length):
    p_remove = .3
    p_keep = 1 - p_remove

    if n == 1:
        return 0.
    if n == history_length:

ma2rten

def p_ben(n, history_length):
    p_remove = .3
    p_keep = 1 - p_remove

    if n == 1:
        return 0.
    if n == history_length:
        return p_keep ** (n-2)

    return p_keep ** (n-2) * p_remove