In the section where the assembly of the linear system is explained,
the RHS of the linear system should be the residual r(u) and not the
primary variable u itself. It should read J(u)*x=r(u).