function R = ged( l, c, N)

ed = @edginess;
cons = @byaverage;
n = @norme;
moy = @moyenne;
gradient = @grad;

Q = [];
R = [];
for i = 1:l-1
    for j = 1:c-1
        e(i,j) = feval( ed, N, i, j);
        d(i,j) = feval( cons, N , i, j);        
        g(i,j) = feval( gradient, N, i, j, 1);
        R(i,j) = e(i,j)*d(i,j)*g(i,j);
        if S(i,j) > 0.0025*k3
    end     
end