function stop=grafico3Dan(nada,state,flag)
%opciones de simulated annelling: opsa =saoptimset('PlotFcns',{@grafico3Dan});
%llamada del algoritmo:[xo,fmin,flag,opciones]=simulannealbnd(@FuncionAlpina,[0,0],[0 0],[10 10],op)

% opciones de algoritmos genéticos: opga =gaoptimset('PlotFcns',{@grafico3D},'MutationFcn', @mutationadaptfeasible,'PopInitRange',[-10;15]);
% llamada del algoritmo:
% [xo,fmin,flag,opciones]=ga(@FuncionAlpina,2,[],[],[],[],[0 0],[10 10],[],opga)


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
X=state.x;
Y=FuncionAlpina(X);
% [X1,X2]=meshgrid(X1,X2);
% meshc(X1,X2,Y,'Marker','*');
plot3(X(1),X(2),Y,'r*');
plot3(X(1),X(2),zeros(length(X),1),'b*');
% if mod(state.iteration,6)==0
%     pause
% end
stop=0;