%% Fonction echantillonneur de Gibbs
%
%  Echantillonne un champ de Gibb avec modele d'Ising en partant de
%  l'initialisation Srs.
%


function Srs=echan(Srs,beta_reg)

for i=1:size(Srs,1),
    for j=1:size(Srs,2),

        
        % voisins V4
        i1=mod(i-2,size(Srs,1))+1;
        i2=mod(i,size(Srs,1))+1;
        j1=mod(j-2,size(Srs,2))+1;
        j2=mod(j,size(Srs,2))+1;
        
        % P(site k = v1)
	% energie conditionnelle locale pour la valeur 0
	% a completer
        Ureg0=double();
	
	% probabilite conditionnelle locale non divisee par Zloc
        % a completer
	p0=exp();
        

        % P(site k = v2)
	% energie conditionnelle locale pour la valeur 1
        % a completer
        Ureg1=double();

	% probabilite conditionnelle locale non divisee par Zloc
	% a completer
        p1=exp();
        

	% Tirage
        if(p0+p1==0)
            break;
        end
        if(rand()<=p0/(p0+p1))
            Srs(i,j)=;
        else
            Srs(i,j)=;
        end
    end
end
