function define_zlim(matriz,zmin,zmax)

if isnan(zmin) || isnan(zmax)
    if (size(matriz,1) == 1) && (size(matriz,2) == 1)
        matriz
        max(matriz)
    end    
    zmax = full(max(max(matriz)));
    zmin = full(min(min(matriz)));
    fprintf(2,'variacao real = [ %24.15e %24.15e]\n', zmin, zmax);
    zmag = max(abs(zmax),abs(zmin));
    zmax = max(0.00001,zmax+0.00001*zmag);
    zmin = min(-0.00001,zmin-0.00001*zmag);
end
zlim([zmin, zmax]);
