I've come across a problem I can't solve. The following differentiation is inconsistently failing to round.
The following answer should have a power of 0.2 NOT 0.1999999999:
This fails to work for x powers of 0.8 and 0.9, but works for 0.1 to 0.7 (see example below):
Has anyone experienced this and can anyone suggest a solution?
a0:ev(decimalplaces(rand(0.8),1)+0.1,simp); a:decimalplaces(a0,1); f(x):=x^a; g(x):=diff(f(x),x);
Many thanks in advance.