clf
% Número de pontos
N = 100;
% População Maxima
Xmax = 1.0;
Ymax = 1.0;

% Cria o domínio da solução
xx = linspace(0,Xmax,N);
yy = linspace(0,Ymax,N);
[XX,YY] = meshgrid(xx,yy);

% Calcula as derivadas nos pontos do domínio
zz = [];
DXX = [];
DYY = [];
for i = 1:N,
  zz = foo(0,[XX(i,:);YY(i,:)]);
  DXX(i,:) = zz(1,:);
  DYY(i,:) = zz(2,:);
end

% Plota o mapa de gradientes
% vetor scale
s = 2.0;
quiver(XX,YY,DXX,DYY,s)
grid
xlabel('x')
ylabel('y')
axis([0,Xmax,0,Ymax])