function h=hcasetta(p,geps)
IndiciTramezze=find(dtramezze(p)<geps);
IndiciScala=find(dscala(p)<30*geps);
IndiciTetto=find(dtetto(p)<geps);

h=6*ones(size(p,1),1);
h(IndiciTetto)=htetto(p(IndiciTetto,:));
h(IndiciTramezze)=htramezze(p(IndiciTramezze,:));
h(IndiciScala)=hscala(p(IndiciScala,:));
