perustulo.py

from pylab import *
from scipy.stats import *

# tilastokeskus 2011
num_adults = 4530e3

# tyottomyyspvrh 2014 * 5 pv / vk * 52 vk
basic_income = 32.66*5*52

# tilastokeskus 2011
labor_force = 4059e3

##
# 2009 luvut
# elake         21,4 mrd
# sairaus       3,83 mrd
# tyottomyysp   3,12 mrd
# peruspv       1,4 mrd
# ulkomaille    2,16 mrd
##
current_wealth_transfers = 31.91e9

def basic_income_cost_benefit():
    direct_costs = num_adults * basic_income
    administrative_cost_per_person = norm(250,75)
    administrative_costs = num_adults * administrative_cost_per_person.rvs()
    
    return direct_costs + administrative_costs

N = 1024*32
bi = zeros(shape=(N,), dtype=float)

for k in range(N):
    bi[k] = basic_income_cost_benefit()

subplot(211)
width = 9e9
height=50*N/1024

title("Basic Income")
hist(bi, bins=50)
axis([0,width,0,height])

show()