for (j=1; j<i; j++)

           if ((s[j]!=-1) && (elp[u][i-j]!=-1))

              z[i] ^= alpha_to[(elp[u][i-j] + s[j])%n];

        // переводим z[i] в индексную форму

        z[i] = index_of[z[i]];

     }

     // вычисление значения ошибок в позициях loc[i]

     //----------------------------------------------------------------------------------------------

     for (i = 0; i<n; i++)

     {

        err[i] = 0;

        // переводим recd[] в полиномиальную форму

        if (recd[i]!=-1) recd[i] = alpha_to[recd[i]]; else recd[i] = 0;

     }

     // сначала вычисляем числитель ошибки

     for (i = 0; i < l[u]; i++)

     {

        err[loc[i]] = 1;

        for (j=1; j<=l[u]; j++)

           if (z[j]!=-1)

              err[loc[i]] ^= alpha_to [(z[j]+j*root[i])%n];

        if (err[loc[i]]!=0)

        {

           err[loc[i]] = index_of[err[loc[i]]];

           q = 0 ; // формируем знаменателькоэффициента ошибки

           for (j=0; j<l[u]; j++)

              if (j!=i)

                   q+=index_of[1^alpha_to[(loc[j]+root[i])%n]];