function state=grafico3D(nada,state,flag)
x=[0:.1:10];
y=x;
for i=1: length(x)
    for j=1:length(y)
        m=[x(i),y(j)];
        z(i,j)=FuncionAlpina(m);
    end
end
clf

[X1,X2]=meshgrid(x,y);
meshc(X1,X2,z)
xlabel('x1')
ylabel('x2')

hold on
X1=state.Population(:,1);
X2=state.Population(:,2);
for i=1: length(X1)
    for j=1:length(X2)
        m=[X1(i),X2(j)];
        Y(i,j)=real(FuncionAlpina(m));
    end
end

% [X1,X2]=meshgrid(X1,X2);
% meshc(X1,X2,Y,'Marker','*');
plot3(X1,X2,Y,'r*');

if mod(state.Generation,6)==0
    pause
end
