function p=blok_prob_szim(A,N,T)
# Az A matrix sorainak a szama a beerkezo folyamok szama, minden sor elso eleme a beerkezes intenzitasat mondja meg, minden sor masodik eleme a kiszolgalasi ido intenzitasat irja le
# Az N szam a kiszolgalo szereverek szamossaga
# A T szam, pedig az idokorlat amennyi idot szimulalunk
# A p erteke pedig az elutasitott csomagok aranya 
  if(columns(A)<2) 
    disp("WARNING: A-nak legyen 2 oszlopa");
    p=-1;
    break; 
  endif;
  n=rows(A); # folyamok szama
  S=[]; # szerverben levo csomagok kiszolgalasi idejenek vege
  t=0; # ido
  t_blok=0; # telitett allapotban eltoltott ido
  packet_all=0; # beerkezett csomagok szama
  packet_block=0; # blokkolt csomagok szama
  for i=1:n
    next_arr(i)=rande()/A(i,1); # a kulonbozo folyamokhoz mikor erkeznek a kovetkezo csomagok
  endfor;
  PX=[0];PY=[0]; # plothoz
  while(t<T)
    if(length(S)==0)  # ha ures a kiszolgalo
      [na,nty]=min(next_arr); # a kovetekzo elem mikor erkezik es melyik folyamhoz tartozik
      S(end+1)=na+rande()/A(nty,2); # az erkezett elem, mikor fog kiszolgalodni
      next_arr(nty)= na+rande()/A(nty,1); # az adott folyamhoz mikor erkezik a kovetkezo elem
      packet_all++; 
      diff=1; # eggyel nott a szerverben levok szama
      ne=na; # kovetkezo esemeny
    elseif(length(S)==N)
      [na,nty]=min(next_arr); # a lovetekzo elem mikor erkezik es melyik folyamhoz tartozik
      [ns,np]=min(S); # a kovetekzo elem mikor szolgalodik ki
      if(na<ns) 
        next_arr(nty)= na+rande()/A(nty,1);
        packet_block++;
        packet_all++;
        diff = 0;
        ne= na; 
      else
        S(np)=[]; # a kiszolgalt elem kivetele a szerverbol
        diff=-1;
        ne=ns;
      endif;
      t_blok+=ne-t;
    else
      [na,nty]=min(next_arr); # a lovetekzo elem mikor erkezik es melyik folyamhoz tartozik
      [ns,np]=min(S); # a kovetekzo elem mikor szolgalodik ki
      if(na<ns) 
        next_arr(nty)= na+rande()/A(nty,1);
        packet_all++;
        S(end+1)=na+rande()/A(nty,2); # az erkezett elem, mikor fog kiszolgalodni
        diff = 1;
        ne= na; 
      else
        S(np)=[];
        diff=-1;
        ne=ns;
      endif;
    endif
      t=ne;
      PX(end+1)=t;
      PY(end+1)=PY(end);
      PX(end+1)=t;
      PY(end+1)=PY(end)+diff;
  endwhile
  if(T<501) newplot(); plot(PX,PY); endif;
  p=packet_block/packet_all;
  disp("ErlangB:");
  disp(erlangB(sum(A(:,1)./A(:,2)),N));
  disp("telitett allapot idoaranya");
  disp(t_blok/t);
endfunction