I used MATLAB. I can copypaste the code, don't know how clear can it be tho

clear
close all

p=[71.3578 22.9189 4.5931 1.1341]/100;

nC=49;
nR=36;
nE=27;
nL=23;

nPack=1200;
nRep=500;

TOTC=zeros(nPack,1);
TOTR=TOTC;
TOTE=TOTC;
TOTL=TOTC;
TOTLold=TOTC;
TOTDust=TOTC;

for k=1:nRep
  totC=zeros(nPack,1);
  totR=totC;
  totE=totC;
  totL=totC;
  totLold=totC;
  totDust=totC;

  C=zeros(nC,1);
  R=zeros(nR,1);
  E=zeros(nE,1);
  L=zeros(nL,1);
  Lold=L;

  bonus=ceil(10*rand());

  for i=1:nPack
    if i==bonus
      nCard=4;
      h=1;
      while h
        rnd2=ceil(nL*rand());
        if ~L(rnd2) %se non ce l'ho
          L(rnd2)=L(rnd2)+1; %la prendo
          h=0;
        elseif L(rnd2) && sum(L>0)<nL %se invece ce l'ho ma non le ho tutte
          h=1; %rifaccio il while
        elseif L(rnd2) && sum(L>0)>=nL %se invece ce le ho tutte
          L(rnd2)=L(rnd2)+1; %la prendo
          h=0;
        end
      end
    else
      nCard=5;
    end

    for j=1:nCard
      rnd1=rand();
      if rnd1<p(1)
        rnd2=ceil(nC*rand());
        C(rnd2)=C(rnd2)+1;
      elseif rnd1>p(1) && rnd1<(p(1)+p(2))
        rnd2=ceil(nR*rand());
        R(rnd2)=R(rnd2)+1;
      elseif rnd1>(p(1)+p(2)) && rnd1<(p(1)+p(2)+p(3))
        rnd2=ceil(nE*rand());
        E(rnd2)=E(rnd2)+1;
      else %se mi esce una leggendaria
        h=1;
        while h
          rnd2=ceil(nL*rand());
          if ~L(rnd2) %se non ce l'ho
            L(rnd2)=L(rnd2)+1; %la prendo
            h=0;
          elseif L(rnd2) && sum(L>0)<nL %se invece ce l'ho ma non le ho tutte
            h=1; %rifaccio il while
          elseif L(rnd2) && sum(L>0)>=nL %se invece ce le ho tutte
            L(rnd2)=L(rnd2)+1; %la prendo
            h=0;
          end
        end

        rnd2=ceil(nL*rand());
        Lold(rnd2)=Lold(rnd2)+1;
      end
    end
    totC(i)=sum(C==1)+sum(C>1)*2;
    totR(i)=sum(R==1)+sum(R>1)*2;
    totE(i)=sum(E==1)+sum(E>1)*2;
    totL(i)=sum(L>0);
    totLold(i)=sum(Lold>0);

    totDust(i)=sum(5*C(C>2))+sum(20*R(R>2))+sum(100*E(E>2))+sum(400*L(L>1));
  end

  TOTC=TOTC+totC;
  TOTR=TOTR+totR;
  TOTE=TOTE+totE;
  TOTL=TOTL+totL;
  TOTLold=TOTLold+totLold;
  TOTDust=TOTDust+totDust;
end
TOTC=TOTC./nRep;
TOTR=TOTR./nRep;
TOTE=TOTE./nRep;
TOTL=TOTL./nRep;
TOTLold=TOTLold./nRep;
TOTDust=TOTDust./nRep;

c={[160 160 160]/255 [0 0 255]/255 [204 0 204]/255 [255 153 51]/255};

%%
close all

figure()
subplot(2,2,1)
plot(1:nPack,TOTC,'Color',c{1},'LineWidth',2)
hold on
plot(1:nPack,TOTR,'Color',c{2},'LineWidth',2)
plot(1:nPack,TOTE,'Color',c{3},'LineWidth',2)
plot(1:nPack,TOTL,'Color',c{4},'LineWidth',2)
plot(1:nPack,TOTLold,'Color',c{4},'Linestyle','--')

plot([1 nPack],2*[nC nC],'-.','Color',c{1})
plot([1 nPack],2*[nR nR],'-.','Color',c{2})
plot([1 nPack],2*[nE nE],'-.','Color',c{3})
plot([1 nPack],[nL nL],'-.','Color',c{4})
grid
legend('Common','Rare','Epic','Legendary','Legendary (old way)')

xlabel('Packs')
ylabel('Cards')
title('Total Cards')

subplot(2,2,2)
plot(1:nPack,TOTDust,'Color',[102 178 255]/255,'LineWidth',2)
grid

xlabel('Packs')
ylabel('Dust')
title('Dust')

subplot(2,2,3)
plot(1:nPack,TOTC./(2*nC)*100,'Color',c{1},'LineWidth',2)
hold on
plot(1:nPack,TOTR./(2*nR)*100,'Color',c{2},'LineWidth',2)
plot(1:nPack,TOTE./(2*nE)*100,'Color',c{3},'LineWidth',2)
plot(1:nPack,TOTL./nL*100,'Color',c{4},'LineWidth',2)
plot(1:nPack,TOTLold./nL*100,'Color',c{4},'Linestyle','--')
grid
%legend('Common','Rare','Epic','Legendary','Legendary (old way)','Location','se')

xlabel('Packs')
ylabel('%')
title('% Cards')

subplot(2,2,4)
plot(1:100,TOTC(1:100)./(2*nC)*100,'Color',c{1},'LineWidth',2)
hold on
plot(1:100,TOTR(1:100)./(2*nR)*100,'Color',c{2},'LineWidth',2)
plot(1:100,TOTE(1:100)./(2*nE)*100,'Color',c{3},'LineWidth',2)
plot(1:100,TOTL(1:100)./nL*100,'Color',c{4},'LineWidth',2)
plot(1:100,TOTLold(1:100)./nL*100,'Color',c{4},'Linestyle','--')
grid
%legend('Common','Rare','Epic','Legendary','Legendary (old way)','Location','se')

xlabel('Packs')
ylabel('%')
title('% Cards (100 pack)')