function define_zlim(matriz)

zmax = max(max(matriz));
zmin = 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);
zlim([zmin, zmax]);
