Sorry, I don't see any bug.

@Markiyan Hirnyk 

I don't understand what workaround you are talking about.
The only issue is that in the help _C is not mentioned; probably it is considered obvious.

@Markiyan Hirnyk 

@Preben Alsholm 

But the comparison is irrelevant because RK4 has a fixed step while in RK45 it is adaptive (as implemented in Maple).

@one pound 

Formally the code works for q := Int(1/ln(t), t = 0 .. x)   too
but it is mathematically incorrect (the new integrals do not exist!).

@one pound 

You may use integration by parts but for
q := Int(1/ln(t), t = exp(1) ..  x);

and then get rid of the constant (which does not matter).

@rlopez 

The problem with this approach is that all the integrals (except q itself) are divergent (even as Cauchy Principal Value) for x>1.

@Kitonum 

No, it is in the source code of GAMMA.

@Kitonum 

This is by design because GAMMA increases very fast.
If you really want such huge numbers then set (not recommended):

`GAMMA/magic`:=100,100,infinity,50;

and use s(100);

@J F Ogilvie 

For many users the xxx__yyy typesetting convention seems to be just enough. I don't see any reason to sacrifice the table functionality for the sake of typesetting. And even without xxx__yyy: why should it be so difficult to use names with more than one letter; must they be exactly the same as in a published article? This is anyway impossible because in many cases in articles the symbols are context-sensitive.

@tomleslie 

But

simplify(subs(z = -z, sin(z)));  # ==>  - sin(z)

(automatic simplification)

It is only (4.)  which intrigues me because simplify usually distributes over operands.

map(simplify, (sqrt(1-I)*sqrt(1+I)-sqrt(2))*(t+1));  # ==> 0

Nice remark. An even shorter procedure is

MaxL:= (a,b,c,r)-> sqrt(a^2+b^2+c^2)-2*r*sqrt((a^2+b^2+c^2)/min(a,b,c)^2-1);

What about 34234-2*17117 ?

BTW, F(a) = Kitonium's f(a),  but you think f is better.

@Kitonum 

This one is simpler and works for symbolic entries too:

F := a -> piecewise(is(a,odd),{(-a+1)/2,(-a-1)/2},{-round(a/2)});

Do you use a random generator for functions of two variables?

@guras