function plot_contour_funcao(parametros,C,p,a,CD,LD,NP) 
%CD centro da regiao e TD tamanho da regiao; supoe que p eh uma grade 
NX = NP(1); assert(floor(nx) == nx);
NY = NP(2); assert(floor(ny) == ny);
xg = transf_vetor_em_matriz(NX,NY,p(:,1));
yg = transf_vetor_em_matriz(NX,NY,p(:,2));
fg = transf_vetor_em_matriz(NX,NY,a);
fprintf(2,'conferir transf_vetor_em_matriz\n');
%hold on
contour(xg,yg,fg,30);
axis equal
xlim([CD(1)-LD(1)/2, CD(1) + LD(1)/2]);
ylim([CD(2)-LD(2)/2, CD(2) + LD(2)/2]);
% for j = 1:nb
%     plot(C(1,j),C(2,j),'.');
% end
