51c51
<     x[n]=x[0]; y[n]=y[0];
---
>     x[n]=x[0]; y[n]=x[0];
66c66
<         nx=(1.0-c2) * x[k] + c2 * x[k+1];
---
>         nx=(1.0-c1) * x[k] + c2 * x[k+1];
83c83
<         nx=(1.0-c2) * x[k] + c2 * x[(k+nm)%n];
---
>         nx=(1.0-c1) * x[k] + c2 * x[(k+nm)%n];
142c142
<       k1 %= n; k2 %= n; k3 %= n;
---
>       k1 %= n; k2 %= n;