Solve() gives wrong answer (r10077)

[parisse]

The majority of my regressions tests are with exact data, this explains these kinds of quirks. I believe I can fix

Code:

fsolve(0.04e-16*x^2-2.4e-16*x+30e-16);

For the other example, it's intrinsically not possible to solve it, as explained by Han. Binary floats can only represent exactly rationals with a denominator dividing a power of 2, that's not the case for 0.04=1/25, therefore rounding errors might transform a 0 discriminant to something negative. csolve will return complex roots in this situation.