Hello Philippe,

I do not think this is a bug. You just bumped into some numerical problems.
For instance, try to increase the number of decimal places:

res←solve(0.05-x-x^2*1/K,x) res=mat(0.00126962,0.00040433,2,1)
x←el(res,1) 0.05-x-x^2*1/K=-0.04391025
x←el(res,2) 0.05-x-x^2*1/K=0.04020022

Not so good, isn't it. That is a problem when you have the function with a rational expression and small value in the denominator (division with K) and the roots close to zero. This is a polynomial of the second degree and when you plot this function you could see that there are two roots, one positive and one negative. Let's try the "roots" this time (guess values are important here):

res1←roots(0.05-X-X^2*1/K,X,-1)
res1=-0.00094148 0.05-res1-res1^2*1/K=2.07062052*10^{-17}
res2←roots(0.05-X-X^2*1/K,X,1)
res2=0.00092408 0.05-res2-res2^2*1/K=-1.47531712*10^{-16}

For this problematic equation, the transforming of the equation will give a slight better results (not much better, but "solve" will did its best):

res←solve(K*(0.05-x)-x^2,x) res=mat(0.00126962,0.00040433,2,1)
K*(0.05-el(res,1))-el(res,1)^2=-0.00000076
K*(0.05-el(res,2))-el(res,2)^2=0.0000007

Moreover, you could try to use the searching interval:
res←solve(K*(0.05-x)-x^2,x,0,1)
res=0.0009240784783
K*(0.05-res)-res^2=3.5504197442949*10^{-16}
(negative root could not find?)

You could try "roots" here as well. The results will be the same (-0.00094148,0.00092408). These are very close to the analytical results.

"solve" and "roots" use the different numerical algorithms and solving the nonlinear equations could be troublesome sometimes. Here you have both roots close to zero and that could make troubles.

Regards,
Radovan

Edited by user 10 November 2010 09:51:35(UTC)  | Reason: Not specified