function [p,q,v,w]=brown(W,nr_iter,sa_nr)
[M,N]=size(W);
p=zeros(1,M);q=zeros(1,N);
sa=zeros(1,N);sb=zeros(M,1);
if nargin==3, sb(sa_nr)=0.0001; end
for i = 1:nr_iter
    %ruch A
    [v_max,ia]=max(sb);
    sa=sa+W(ia,:);
    p(ia)=p(ia)+1;
    %ruch B
    [v_min,ib]=min(sa);
    sb=sb+W(:,ib);
    q(ib)=q(ib)+1;
end
p=p/nr_iter;q=q/nr_iter;
v(1)=min(sa)/nr_iter;v(2)=max(sa)/nr_iter;
w(1)=min(sb)/nr_iter;w(2)=max(sb)/nr_iter;