
clc;  clear;  close
load ("data/Parameters.mat",'m');

%% The equilibrium
% Case 1: competition + prohibit
m.DataPolicy = 'prohibit';  % prohibit free
m.MarketStruct = 'competition';  % competition oligopoly monopoly 
m.Privacy = 'nonePrivacy';  % nonePrivacy privacyExternality

tic
[wel1,Case1] = func_eqm_iter(m);
toc

% Case 2: competition + free
m.DataPolicy = 'free';  % prohibit free
m.MarketStruct = 'competition';  % competition oligopoly monopoly 
m.Privacy = 'nonePrivacy';  % nonePrivacy privacyExternality

tic
[wel2,Case2] = func_eqm_iter(m);
toc

% Case 3: monopoly + prohibit
m.DataPolicy = 'prohibit';  % prohibit free
m.MarketStruct = 'monopoly';  % competition oligopoly monopoly 
m.Privacy = 'nonePrivacy';  % nonePrivacy privacyExternality

tic
[wel3,Case3] = func_eqm_iter(m);
toc

% Case 4: monopoly + free
m.DataPolicy = 'free';  % prohibit free
m.MarketStruct = 'monopoly';  % competition oligopoly monopoly 
m.Privacy = 'nonePrivacy';  % nonePrivacy privacyExternality

tic
[wel4,Case4] = func_eqm_iter(m);
toc

save("data/Baseline.mat",'Case1','Case2','Case3','Case4');
