resolution=0.001;
t=0:resolution:10;
nu=0;
input=cos(2*pi*nu*t);
output=identme(input,resolution);
H=max(output)

for nu=0:40 
input=cos(2*pi*nu*t);
output=identme(input,resolution);
H=[H,max(output
end