function [etheta_beam1,edelta_beam1,etheta_beam2,edelta_beam2,etheta_dspe1,edelta_dspe1,etheta_dspe2,edelta_dspe2]=bearfrming1good(default)
global SNR msensk bearing1 spread bearing2 delta
 
if (exist('default') ~= 1) default = 0;  else default = 1; end
def_msens = 8;%陣元數(shù)
def_dspace = 0.5;%
def_nsamp =500;%采樣數(shù)
def_msource =2;%

%extended=0.0175;
in_type = 'lap';
msens=def_msens;
nsamp=def_nsamp;
nvar = 1;
nsource=def_msource;
ar = [1 0 0.81];
ma = 1;
bearing1=10;
spread=4;
bearing2=25;
delta=6;

bearing1=bearing1*pi/180;
bearing2=bearing2*pi/180;
for j=1:2
    smat(:,j) =10.*rpiid(def_nsamp,in_type);
end  

%Gaussian distribution
for msensk=1:msens
    [b(msensk)]=quad('myjifenmusic2D',-pi/2,pi/2);
end

amat2=[b];


for msensk=1:msens
    [b1(msensk)]=quad('myjifenmusic2D1',-pi/2,pi/2);
end

amat3=[b1];



amat=[amat2.' amat3.'];
msens = def_msens; dspace = def_dspace; nsamp = def_nsamp;
smat=smat*amat.';    
  % --- Convert ARMA noise to equivalent MA form:
  if (nvar > 0)
     cvec = cumtrue(ma,ar,2,msens);
     q = (length(cvec) - 1) /2 ;
     if (q < msens)
        cvec = [zeros(msens-q,1); cvec; zeros(msens-q,1)];
     end
     cmat = toeplitz(cvec(msens+1:2*msens));
     amatn = chol (cmat);

% --- Generate white noise, then color it:
     gmat = reshape(rpiid(nsamp*msens,'nor'),nsamp,msens);
     gmat = gmat * amatn;
     lby21 = fix((nsamp-1)/2);
     mask  = [1; ones(lby21,1)*2;1; zeros(lby21,1)];
     for i=1:msens
        gmat(:,i) = ifft( fft(gmat(:,i)).*mask );        % complex noise
	g = gmat(:,i) - mean(gmat(:,i));
	q = sqrt( real(g'*g) / (nsamp-1) );
	gmat(:,i) = gmat(:,i) / q * sqrt(nvar);
     end
     smat = smat + gmat;
  end
  ymat=smat;
  cmat=zeros(msens,msens);
  bb=ones(nsamp,1)*mean(ymat);
  
  ymat=ymat-bb;
  
  cmat=conj(ymat'*ymat)/nsamp;
 
  
  
  
  
  
  
  
  
  
  
  

 return