Нахождение экстремумов числового ряда

input int      bars=500;                 //  окно для поиска экстремумов
input double   delta_points=160;         //  размах перепада, определяющий минимальное расстояние между пиком и впадиной в пунктах
input double   first_extrem=0.9;         //  дополнительный  коэффициент для поиска первого экстремума
input double   reload_time=5;            //  значение временного интервала, через который производится перерасчет значений индикатора в секундах
 
double High[],Low[];
datetime Time[];
 
ArraySetAsSeries(Low,true);
int copied1=CopyLow(Symbol(),0,0,bars+2,Low);
ArraySetAsSeries(High,true);
int copied2=CopyHigh(Symbol(),0,0,bars+2,High);
ArraySetAsSeries(Time,true);
int copied3=CopyTime(Symbol(),0,0,bars+2,Time);
 
double delta=delta_points*Point();  //  размах перепада между экстремумами в абсолютных величинах
 
int j,k,l;
int j2,k2,l2;
double  j1,k1,l1;
int min[6];  // массив, определяющий впадины, значение соответствует номеру бара для найденного экстремума
int max[6];  // массив, определяющий пики, значение соответствует номеру бара для найденного экстремума
 
int mag1=bars;
int mag2=bars;
int mag3=bars;
int mag4=bars;
 
j1=SymbolInfoDouble(Symbol(),SYMBOL_BID)+(1-first_extrem)*delta_points*Point();
// при поиске первого экстремума дополнительный коэффициент определяет минимальную цену, ниже которой должна располагаться первая впадина
 
j2=0; // на первой итерации поиск ведется начиная с последнего бара истории
 
for(j=0;j<=15;j++) // цикл, определяющий первую впадину - min[1]
  {
   min[1]=minimum(j2,bars,j1);
   //определяется ближайшая впадина на указанном интервале
 
   j2=min[1]+1;     // на следующей итерации поиск ведется от уже найденной впадины min[1]
   j1=Low[min[1]]+delta;
   //минимальная цена для впадины, найденной на последующей итерации должна быть ниже, чем минимальная цена для впадины, найденной на текущей итерации 
 
   k1=Low[min[1]];
   //минимальная цена для впадины при поиске последующей вершины определяет максимальную цену, выше которой должна располагаться эта вершина
 
   k2=min[1];         //поиск пика, располагающегося за впадиной, ведется от найденной впадины min[1]
 
   for(k=0;k<=12;k++) // цикл, определяющий первый пик - max[1]
     {
      max[1]=maximum(k2,bars,k1);
      //--- определяется ближайший пик на указанном интервале
      k1=High[max[1]]-delta;
      //максимальная цена для пика, найденного на последующей итерации должна быть выше, чем максимальная цена для пика, найденного на текущей итерации 
 
      k2=max[1]+1;   // на следующей итерации поиск ведется от уже найденного пика max[1]
 
      l1=High[max[1]];
      //максимальная цена для вершины при поиске последующей впадины определяет минимальную цену, ниже которой должна располагаться эта впадина
      l2=max[1];     // поиск впадины, располагающейся за пиком, ведется от  найденного пика max[1]
      for(l=0;l<=10;l++) // цикл, определяющий вторую впадину - min[2], и второй пик max[2]
        {
         min[2]=minimum(l2,bars,l1);
         //---определяется ближайшая впадина на указанном интервале
         l1=Low[min[2]]+delta;
         //минимальная цена для впадины, найденной на последующей итерации должна быть ниже, чем минимальная цена для впадины, найденной на текущей итерации 
 
         l2=min[2]+1;     // на следующей итерации поиск ведется от уже найденной впадины min[2]
         max[2]=maximum(min[2],bars,Low[min[2]]);
         //определяется ближайший пик на указанном интервале
 
         if(max[1]>min[1] && min[1]>0 && min[2]>max[1] && min[2]<max[2] && max[2]<mag4)
           //отфильтровываются совпадающие экстремумы и особые случаи
           {
            mag1=min[1];   // на каждой итерации в случае выполнения условия, положения найденных экстремумов запоминаются
            mag2=max[1];
            mag3=min[2];
            mag4=max[2];
           }
        }
     }
  }
min[1]=mag1; // экстремумы определены, в ином случае всем переменным присваивается значение bars
max[1]=mag2;
min[2]=mag3;
max[2]=mag4;
int minimum(int a,int b,double price0)
//функция определяет ближайшую впадину на указанном интервале, которая располагается ниже цены price0 на расстоянии большем, чем размах перепада
  {
   double High[],Low[];
   ArraySetAsSeries(Low,true);
   int copied4=CopyLow(Symbol(),0,0,bars+2,Low);
 
   int i,e;
   e=bars;
   double pr=price0-delta_points*Point();    // значение цены ниже  которой должна располагаться впадина уже с добавлением размаха перепада
   for(i=a;i<=b;i++)                         // поиск впадины  в окне заданном параметрами a и b
     {
      if(Low[i]<pr && Low[i]<Low[i+1])       // определяется ближайшая впадина,  после которой начинается рост цены 
        {
         e=i;
         break;
        }
     }
 
   return(e);
  }
int maximum(int a,int b,double price1)
//--- функция определяет ближайший пик на указанном интервале, который располагается выше цены price1 на расстоянии большем, чем размах перепада
  {
   double High[],Low[];
   ArraySetAsSeries(High,true);
   int copied5=CopyHigh(Symbol(),0,0,bars+2,High);
 
   int i,e;
   e=bars;
   double pr1=price1+delta_points*Point();   // значение цены выше  которой должен располагаться пик уже с добавлением размаха перепада
   for(i=a;i<=b;i++)                         // поиск пика в окне заданном параметрами a и b
     {
      if(High[i]>pr1 && High[i]>High[i+1])   // определяется ближайший пик, после которого начинается падение цены     
        {
         e=i;
         break;
        }
     }
   return(e);
  }