function fmedia = media_celular(p,f,C)
%calcula a media celula para fins do metodo IC, dadas as posicoes p dos
%pontos de amostragem desta celula, os valores f da funcao nesses pontos e
%o centro C da celula.

np = size(p,1);
nf = size(f,2);
assert(size(f,1) == np);
assert(size(p,2) == 2);
assert(size(C,2) == 2);
assert(size(C,1) == 1);

for j = 1:nf
    fmedia(j) = sum(f(:,j))/np;
end