N = 5000;

rho = 3;

Rmax = 50;

total_revenues = zeros(length(N),1);

for n = 1:N

R = zeros(24,1);

x = zeros(24,1);

P_d = zeros(24,1); P_d(1) = 50;

P_s = zeros(25,1);

W= zeros(25,1);

revenues = zeros(24,1);

for t = 1:24

% Commitment

x(t) = (R(t)+10)*rho;

% Wind

W(t+1) = unifrnd(5,15);

%

P_s t=1 = P_d t = 32;

%

= W t=1 - x t0;

= x t - W t=10

> 0

R t =1 = R t = Rmax;

t = x t * P_d t;

< R t

R t=1 = R t - ;

t = x t * P_d t;

R t=1 = 0;

t = x t * P_d t - - R t * P_s t = 1;

%

01 < 0.5

= 1;

= -1;

P_d t=1 = P_d t = ;

total_revenuesn = ;

total_revenues

--------------------------------------------------------------------------------------------------------------

%%

k_1 = [0.66 0.65 0.64 0.63 0.62 0.61 0.60 0.59 0.55];

c_1 = [0.085855 0.032191 0.020795 0.017001 0.013711 0.010851 0.008388 0.006291 0.001401];

k_2 = [1.30 1.25 1.20 1.15 1.10 1.05 1.00 0.95 0.90];

c_2 = [0.137213 0.082645 0.045060 0.028338 0.016146 0.007860 0.003277 0.001134 0.000245];

= k_1k_2;

n = 100000;

i = 1k_1

j = 1k_2

:)