%  Batch file for ggm_queue and ggm_queue_contrvar
 a1=0;a2=10; b1=10; b2=20; s=4;T=1000;cL=-0.04; cW=-0.03; N=1000;
 % Tables for original Monte Carlo scenaria
 L1=zeros(N,1); W1=zeros(N,1);
 % Tables for Control Variates Monte Carlo scenaria
 L2=zeros(N,1); W2=zeros(N,1);
 L3=zeros(N,1); W3=zeros(N,1);
 NA=zeros(N,1);
 for i = 1:N
 	[L1(i), W1(i)]=ggm_queue(a1, a2, b1, b2, s, T);
	[L2(i), W2(i),L3(i), W3(i), NA(i)]=ggm_queue_contvar(a1, a2, b1, b2, s, T,cL, cW);	
 end;
 L1m=mean(L1)
 L1s=std(L1)
 L2m=mean(L2)
 L2s=std(L2)
 L3m=mean(L3)
 L3s=std(L3)
  W1m=mean(W1)
 W1s=std(W1)
 W2m=mean(W2)
 W2s=std(W2)
 W3m=mean(W3)
 W3s=std(W3)
 NAm=mean(NA)
 NAst=std(NA)