aveflx1 = aveflx(1:ipeak1);

aveflx2 = aveflx(ipeak1:ipeak2);

aveflx3 = aveflx(ipeak2:utslen); 

[f1 fi1] = min(abs(aveflx1-peakfall(1)));
[f2 fi2] = min(abs(aveflx2-peakfall(1)));
fi2 = fi2 + ipeak1;
[f3 fi3] = min(abs(aveflx2-peakfall(2)));
fi3 = fi3 + ipeak1;
[f4 fi4] = min(abs(aveflx3-peakfall(2)));
fi4 = fi4 + ipeak2;

f = [fi1 fi2 fi3 fi4];