clear rows cols

maxSteps = steps;
minSteps = 1;
for stepsIndex=minSteps:maxSteps
 
    
nRange=stepsIndex;
c1Range=stepsIndex;

C=max(qdot(1:nRange,1:c1Range));
D=max(C);
[rows(stepsIndex), cols(stepsIndex)] = find(qdot(1:nRange, 1:c1Range) == D);

end

figure(1)
plot(rows, cols)
title('Peak Qdot')
xlabel('n')
ylabel('c1')