All these are special (particular) cases.
It is not possible to express the general solution in terms of real functions.
Just think at the much more simple case of the Laplace equation  diff(f(x,y),x,x) + diff(f(x,y),y,y) = 0.
Or, even better, at the Cauchy-Riemann equations: 
diff(f(x,y),x) - diff(g(x,y),y) = 0, diff(f(x,y),y) + diff(g(x,y),x) = 0.

@Carl Love 

The 9th hex digit is not always correct (I have checked this). I'd recommend to return only 8 hex digits in the procedure BBP (i.e. ModuleApply).

Edit: unfortunately the 8th digit is also sometimes incorrect. Please revise the code.

Edit2: I think I found the problem. ModuleApply must have

Digits:= trunc(evalhf(Digits)) + 3;

Note that +2  is not enough.
Now the results are correct at least for d <= 5000.
The 9th digit is still wrong.
I am not sure about d>5000  :-(

@Carl Love 

Nice! Vote up!

P.S. Why do you say that  16 &^ 0  mod  4;   is a bug?  It is indeed  0 &^ 0  mod 4; 
(so, a convention, even if 0^0 = 1).

The fact that the hypercube [0,1]^d is densifiable is obvious: take the (m+1)^d points having all coordinates in the set {0,1/m,2/m,...,m/m} and take any polygonal curve passing thru all of them. With this approach, the "t" can be found at once.

Of course, the "cosine" curve g is more elegant, but the problem is that the "proof" you mentioned is only a sketch. It is not very hard to obtain a real proof such that given any of the m^d small hypercubes H, find (effectively) t in [0,1] for which g(t) \in H.
This is the way to go, and the Maple code is not needed then.

So, you have a kind of "almost" space filling curve  (see wiki).
- What is the source for this curve?
- Do you know to prove mathematically the claimed property? Or you want some kind of "incomplete" proof using Maple?
Or maybe you want to use this curve for something else?

@tomleslie 

It is not difficult to explain.
When SolveEquations solves f(x)=0, it actually minimizes an associated function F(x)  (such as F(x) = f(x)^2).

A minimization with constraints g(x)<=0  is much more difficult (mainly for several variables). It is not possible in general just to minimize F and then select the solutions satisfying g(x)<=0  because in most cases the set of such solutions could be empty. It is faster for the user to select himself the desired solutions if he knows that such solutions exist and they are not too many.

@Earl 

You have mentioned the tangency. In this case of a not differentiable curve it should be interpreted as having at least a common point and disjoint interiors. This can be done but it's of course more complicated.

So, without shades of grey and without thickness? Then make a custom legend using textplot, line, pointplot.

@weidade37211 

@Earl 

If you really want exact tangency with a circle  you should solve the problem separately.
Since the curve is not differentiable for t=k*Pi (k in Z), these points need special atention, fsolve/DS only is not enough.

@tomleslie 

Combining yours with mine:

r1:=t->piecewise(t<Pi, sqrt(exp((t+(1/4)*Pi)*(1/2))),
                t<2*Pi, sqrt(exp((9*Pi/4-t)*(1/2)))
               ):
r:=t->r1(t - 2*Pi*floor(t/(2*Pi))):
plots:-animate(plot, 
[[[r(t-u)*cos(t),r(t-u)*sin(t),t=0..2*Pi] ,
[[0,r(Pi/2-u)*sin(Pi/2)],[0,2+r(Pi/2-u)*sin(Pi/2)]]   ],thickness=[1,6]],
u=0..2*Pi, scaling=constrained);

@Adam Ledger 

For this specific s()  it does not matter too much; it's only a shift of the sequence with one position. 
Of course for a recurrence  s(n) = f(n, s(n-1)) rather than  s(n) = f(s(n-1)),  the starting index could be important.

Just a remark: this way the climber will never reach the top of the mountain :-)   [it will take an ethernity!]

You should isolate the problem. You have a large execution group in 2D input which is difficult (time consumming) to split/manage.

@tsunamiBTP 

Your equation Q1=0  is polynomial of degree m in X = cos(2*Pi*x/T).
So, it cannot be solved analitically except possibly in some very special situations; you will have to compute the Galois group in order to find these situations. This is very complicated and useless for this problem.

P.S. Beying polynomial, you don't need RootFinding; you can simply find with solve all the roots of the polynomial (numerically).

Anyway the splash screen is essential because the Java interface is slow (even on decent computers) and the user has to see something until the worksheet appears.