%-----------------------------------------------------------
function z=sbiniz(dummy)
  global im jm km lm
  x1=0.4;
  y1=0.4;
  x2=0.6;
  y2=0.6;
  sigma=0.05;
  
  for i=1:im
    x=(i-1)/im;
    for j=1:jm
      y=(j-1)/jm;
      g(j,i)=exp(-((x-x1)^2+(y-y1)^2)/(2*sigma^2))+exp(-((x-x2)^2+(y-y2)^2)/(2*sigma^2));
    end
  end
  
  z=p2g2sa(lm,km,g);
  z(1+lm,1+km)=0;
  
