How could you find the index number of an answer output?

I have a code that looks like:

for ii=2:length(dep);

    dz(ii) = dep(ii)-dep(ii-1);

    drho(ii) = sigth(ii)-sigth(ii-1);    

end

for ii=1:length(dz);

      drhodz(ii) = drho(ii)/dz(ii);  
  end

max(drhodz)

Is there a way to go back and then see which ii this maximum answer corresponds to? (Does that question make sense? After finding the maximum drhodz I'm trying to work backwards to see at which dz (change in depth) this max occurs.)

Answer 1

[Arnaud Miege]

[max_value, max_idx] = max(drhodz)

The index ii corresponding to the maximum value max_value is given by max_idx.

HTH,

Arnaud

Answer 2

[Sean de Wolski]

Arnaud's reply is correct

To find the indices of multiple maxima:

find(drhodz==max(drhodz));

Also, your for-loops can be easily vectorized with:

dz = diff(dep);
drho = diff(sigth);
drhodz = dz./drho;

[Chelsey]

Thank you - a much more elegant solution than what I was trying to do.

[Arnaud Miege]

Good point!!

Answer 3

[Chelsey]

Thank you both!