function ploskev(B,k);
 
  [n,m]=size(B);
  u=linspace(0,1,k);
  v=linspace(0,1,k);
  V=zeros(k,n);
  U=zeros(m,k);
  for i=1:n
    V(:,i)=bern(n-1,i-1,v)';
  end
  for j=1:m
    U(j,:)=bern(m-1,j-1,u);
  end
  Z=V*B*U;
  clf
  mesh(u,v,Z);
  
   
