In [1]:

    

import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline  

import matplotlib.style as style
style.use('bmh')


    

In [2]:

    

p = np.array([9, 4]) # p[0] > 2p[1]                                                                                                                                                                             
f = np.array([1, 1. / p[1], 1. / p[0]])
n = 3
p3 = np.array([2] * n)
x = np.append(p.reshape(2, 1) * f.reshape(1, n), p3.reshape(1, n), axis=0)


    

In [5]:

    

w = .3
idx = np.arange(len(x[0]))
opacity = 0.4

colors = ['r', 'g', 'b']
labels = ['$p_{MOX, a}$', '$p_{UOX, a}$', '$p_{UOX, b}$']

for i in range(3):
    plt.bar(idx + i * w, x[i], w, color=colors[i], label=labels[i], alpha=opacity)
plt.ylabel('Preference (Arbitrary Units)')
plt.xticks(idx + 1.5 * w, ('$t < t_0$', '$t_0 \leq t < t_1$', '$t_1 \leq t$', '$t_2 \leq t < t_3$', '$t \geq t_3$'))
plt.legend()
plt.savefig('../figs/tariff_prefs.pdf')
plt.savefig('../figs/tariff_prefs.png')


    

In [12]:

    




    

In [ ]: