% salgstall fra BROD, pr måned 
data2017 = [  726320  745551   996583  890146 ...
	     1037795 1160839  1113303  977774 ...
	      924832 1103728  1145418  838077 ] ;
data2018 = [  674622  706994   860976  835847 ...
	     1038115  998085  1002583  930458 ...
	      836671  936679   830854       0 ] ;  

% andelen av Lervigs produksjon som går til eksport i 2018
eksport  = 50 ; 

fprintf('Ølsalg småskala i 2017 og 2018\n') ; 
fprintf('%5.0f ', data2017/1000)
fprintf('\n'); 
fprintf('%5.0f ', data2018/1000) ;   
fprintf('\n'); 


% først ekstrapolerer vi des 2018 med tall fra 2017
data2018(12) = data2017(12) ;
% så finner vi månedenes prosentvise andel av salget i 2017
pct17 = data2017/sum(data2017)*100 ;
fprintf('\nDen prosentvise fordeling av småskala salg i 2017\n');
fprintf('%5.1f ', pct17) ; 
fprintf('\n'); 

% Lervigs 2018-salg med antatt månedlig fordeling à la 2017
lervig = 1750000 * (1-eksport/100) * pct17/100 ;
fprintf('\nLervigs antatte salg i 2018, med fordeling som småskala salg i 2017\n');
fprintf('%5.0f ', lervig/1000) ; 
fprintf('\n'); 

% første approx, legg til salgstall etter 2017-fordeling av salg
data2018b = data2018 ; 
data2018b(6:11) = data2018b(6:11) + lervig(6:11) ;
fprintf('\nTall for 2018 med Lervig, gitt fordeling som i 2017\n') ; 
fprintf('%5.0f ', data2018b/1000) ; 
fprintf('\n'); 

% andre approx, bruk dette som salgsfordeling for 2018, istf 2017-fordelingen
pct18 = data2018b/sum(data2018b)*100 ;
fprintf('\nDen prosentvise fordeling av antatt totalt salg i 2018 med estimater for Lervig\n');
fprintf('%5.1f ', pct18) ; 
fprintf('\n'); 

lervig = 1750000 * (1-eksport/100) * pct18/100 ; 
fprintf('\nLervigs antatte salg i 2018, med fordeling som antatt småskala salg i 2018\n');
fprintf('%5.0f ', lervig/1000) ; 
fprintf('\n'); 

data2018c = data2018 ; 
data2018c(6:11) = data2018c(6:11) + lervig(6:11) ;  
fprintf('\nTall for 2018 med Lervig, gitt fordeling som for antatte tall for 2018\n') ; 
fprintf('%5.0f ', data2018c/1000) ; 
fprintf('\n'); 

diff2018b = (data2018b-data2017)./data2017*100 ;
diff2018c = (data2018c-data2017)./data2017*100 ;
fprintf('\nProsentvis endring i salg fra 2017 til de to estimatene for 2018\n') ; 
fprintf('%5.1f ', diff2018b) ; 
fprintf('\n'); 
fprintf('%5.1f ', diff2018c) ; 
fprintf('\n');