
2D elliptic PDE on unit square
   d^2 u / dx^2 + d^2 u / dy^2 = u^3 - u + 2.0
 + homogeneous Dirichlet boundary conditions

Solution method: Modified Newton with band linear solver
Problem size: 31 x 31 =  961

Computed solution (||F|| = 1.3498e-12):

            0.03125  0.12500  0.21875  0.31250  0.40625  0.50000  0.59375  0.68750  0.78125  0.87500  0.96875  

0.03125     0.00405  0.01165  0.01617  0.01896  0.02051  0.02100  0.02051  0.01896  0.01617  0.01165  0.00405  
0.12500     0.01165  0.03772  0.05461  0.06530  0.07126  0.07318  0.07126  0.06530  0.05461  0.03772  0.01165  
0.21875     0.01617  0.05461  0.08098  0.09813  0.10780  0.11093  0.10780  0.09813  0.08098  0.05461  0.01617  
0.31250     0.01896  0.06530  0.09813  0.11989  0.13229  0.13631  0.13229  0.11989  0.09813  0.06530  0.01896  
0.40625     0.02051  0.07126  0.10780  0.13229  0.14632  0.15089  0.14632  0.13229  0.10780  0.07126  0.02051  
0.50000     0.02100  0.07318  0.11093  0.13631  0.15089  0.15564  0.15089  0.13631  0.11093  0.07318  0.02100  
0.59375     0.02051  0.07126  0.10780  0.13229  0.14632  0.15089  0.14632  0.13229  0.10780  0.07126  0.02051  
0.68750     0.01896  0.06530  0.09813  0.11989  0.13229  0.13631  0.13229  0.11989  0.09813  0.06530  0.01896  
0.78125     0.01617  0.05461  0.08098  0.09813  0.10780  0.11093  0.10780  0.09813  0.08098  0.05461  0.01617  
0.87500     0.01165  0.03772  0.05461  0.06530  0.07126  0.07318  0.07126  0.06530  0.05461  0.03772  0.01165  
0.96875     0.00405  0.01165  0.01617  0.01896  0.02051  0.02100  0.02051  0.01896  0.01617  0.01165  0.00405  

Final Statistics.. 

nni      =      7    nfe     =      8 
nbcfails =      0    nbacktr =      0 
nje      =      1    nfeB    =     63 

lenrw    =   4822    leniw   =     27 
lenrwB   = 150877    leniwB  =    961 
