Hola: Yo había tenido un problema similar. En mi caso me había equivocado al copiar la definición de X^k+1. Fijate que durante la iteración se saltea el valor actual. Es decir va por hasta j-1 y después arranca en j+1. El valor j viene afuera.

Saludos,
Jerónimo