[Help-gsl] Should be 4 roots for non-linear equations

[Sergey Shcherbina]

Hello!

I have been trying to use gsl for solving a non-linear system equations - 
https://www.brainfuse.com/jsp/alc/resource.jsp?s=gre&c=36961&cc=108828
I took first example from this web-page and what I have as a result:

iter =   0 x =  1.000  1.000 f(x) = -4.000e+00 -1.400e+01
iter =   1 x =  1.000  1.000 f(x) = -4.000e+00 -1.400e+01
iter =   2 x =  1.212 -5.196 f(x) = -9.726e+00 -7.592e+00
iter =   3 x =  3.965 -9.619 f(x) =  9.858e-02 -4.166e-01
iter =   4 x =  4.021 -9.979 f(x) =  1.853e-01  0.000e+00
iter =   5 x =  3.990 -10.010 f(x) = -8.619e-02  0.000e+00
iter =   6 x =  4.000 -10.000 f(x) = -1.967e-04  0.000e+00
iter =   7 x =  4.000 -10.000 f(x) =  2.097e-07  0.000e+00
iter =   8 x =  4.000 -10.000 f(x) = -5.080e-13  0.000e+00

In result we see two first numbers only - (4.0, -10.0). But I don't know what 
should be changed in attached C code for calculation other two roots (-5,19). 
Or it is limited by development?  

Regards,
Sergey.

TestxGSLxMATHx002a.c
Description: Binary data