Possibility of solving a certain 4D system with 4 equations with certain methods Can the system 
\begin{align}
x+y& =a\\
 x^2+b& =z^2\\
y^2+c&=w^2\\
dzy&=wx
\end{align} 
(with $x,y,z,w,a,b,c,d >0$) be solved for $x, y, w$ and $z$ by hand without cubics (which I could solve... but yikes) or quartics, and using only elementary functions? I was challenged to find a simple solution but after alot of manipulating I'm starting to think it's not possible. I know that once you solve for one variable you can use one of the first three equations to isolate another easily with only roots and etc., and then do that again with the variable you just found until you have them all. The trouble is in getting just one of these variables. 
I've tried squaring both sides of the fourth equation and writing in terms of x and y, then using the first equation to get just x, by that left me with a bad quartic. I also though because of the elegance of the second and third equation I'd try some hyperbolic trig subs with two auxilarily variables but that didn't get anything new either. Neither did adding some of these equations, as I assume this is because you're using inputs and information from the same two equations. 
I don't mind if Galois theory or anything is used in the answer, but I probably won't understand it (at least, not yet ;)). 
 A: Let us reduce the issue to the obtention of a root $x_0$ of a certain equation $f(x)=0$ (see figure below).
Knowing that your issue is oriented towards applications, a final step is to use a numerical method such as dichotomy (or Newton's method) to find an approximation to $x_0 \in (0,a)$ ; knowing $x_0$, one will get the unique values of $y$, $z$ and $w$. We will leave this numerical solution up to you.
Let us explain how one arrives at this function $f$. 
The different variables $y,z,w$ in system :
$$\begin{cases}
x+y& =&a\\
x^2+b& =&z^2\\
y^2+c&=&w^2\\
dzy&=&wx
\end{cases}$$
can be progressively transformed into expressions involving $x$ only (and parameters $a,b,c,d$). 
$y,z$ and $w$ can be expressed as functions of $x$ in the following way :
$$y=a-x, z=\sqrt{x^2+b}, w=\sqrt{y^2+c}=\sqrt{(a-x)^2+c}$$
(no $\pm$ sign in front of the square roots because $z$ and $w$ are positive) ; inserting these expressions into the fourth equation gives the constraint :
$$d\sqrt{x^2+b}(a-x)=x\sqrt{(a-x)^2+c}$$
which is equivalent, under the condition $a>x$, to :
$$d^2(x^2+b)(a-x)^2=x^2((a-x)^2+c)$$
or
$$d^2\dfrac{x^2+b}{x^2}=\dfrac{(a-x)^2+c}{(a-x)^2} \ \iff$$
$$d^2(1+\dfrac{b}{x^2})=1+\dfrac{c}{(a-x)^2} \ \iff $$
$$\underbrace{(d^2-1)+\dfrac{bd^2}{x^2}-\dfrac{c}{(a-x)^2}}_{f(x)}=0.\tag{1}$$
As
$$f'(x)=-\dfrac{bd^2}{x^3}-\dfrac{c}{(a-x)^3}<0 \ \ \text{for} \ \ x \in (0,a), $$
the curve of $f$ is decreasing on $(0,a)$ from $+\infty$ to $-\infty$, thus with a unique root in this interval. See the graphical representation in the example below, representative of the general case. We don't consider the other branch because $x>a$ isn't allowed.

Fig. 1 : Graphical representation of $f$ (Case $a=5;b=2;c=3;d=2$) with solution $x_0 \approx 4.07186$. Only the left branch is to be considered, because of condition $y=x-a>0$. 
Remark : $f(x)=0$ is equivalent to a 4th degree equation :
$$(d^2-1)x^2(a-x)^2+bd^2(a-x)^2-cx^2 \ = \ 0$$
A: As you didn't show your results, I'll redo the computation.
Let us substitue the RHS of the two middle equations in the fourth, squared:
$$d^2(x^2+b)y^2=(y^2+c)x^2,$$ or
$$(d^2b-x^2)y^2=(c-d^2)x^2.$$
Then, from the first equation,
$$y^2=(a-x)^2$$
and finally, eliminating $y$,
$$(c-d^2)x^2=(d^2b-x^2)(a-x)^2.$$
This confirms that the problem is of the quartic degree and has an analytical solution, but the formulas are terrible and hardly manageable by hand.
And as the polynomial has three independent coefficients, there is no reason that any simplification occurs and you won't find any simpler analytical resolution (otherwise, you would revolutionize the Galoi's theory).

The other option is to resort to numerical methods, but whether this is feasible by hand is unsure, and it is only doable for particular values of the parameters $a,b,c,d$.
