How to define a surface which includes three coordinate axes? Imagine a plane with two crossing lines along $x$ and $y$ axes, and one line $y=x$. Now bend the $y=x$ line towards $z$ axis, dragging the part of that plane with it. When this line coincides with $z$ axis, you'll get a surface similar to this:

But this one was generated from $z=\text{gaussian}_{0.3a}(x)$, i.e. if continued it'd not have $z$ full axis as its subset.
I'd like the surface to be:


*

*Differentiable (except maybe at origin)

*Symmetric in the sense that for each axis $A$ there existed a rotation which placed another axis $B$ to original location of $A$, leaving the surface unchanged.

*Appeared as a single bent & stretched plane


How could I find (analytically define) such a surface (or a family of surfaces)?

As requested by OP, an answer in stead of a comment. My answer is simply an intuitive approach to finding a single function with the desired properties, and does not at all hint at what the class of desired functions might look like.
As I know close to nothing about graphs or differentiable things, I restrict my attention to polynomials, things I am more familiar with. The symmetry requirement suggests considering a sum of a polynomial $f\in k[X,Y,Z]$ over its orbit under the action of $S_3$ or $A_3$, that is
$$F_S=\sum_{\sigma\in S_3}\sigma(f)\qquad\text{ or }\quad F_A=\sum_{\sigma\in A_3}\sigma(f).$$
The constructed polynomial is (almost) symmetric in $X$, $Y$ and $Z$, so for the graphs of $F_S$ and $F_A$ to contain the coordinate axes it suffices for $f$ to be zero on the $Z$-axis. Without loss of generality we may assume $f\in(XY)$. The simplest example would then be $f=XY$, yielding
$$F_A=XY+XZ+YZ=\tfrac{1}{2}\cdot F_S,$$
which is everywhere differentiable except at the origin. Unfortunately this does not meet the third requirement. To meet this requirement, it makes sense to consider a polynomial $f$ of odd degree. 
The odd degree will make sure that the evaluation maps $k^2\ \longrightarrow\ k$ are surjective, for any of the three functions that $f$ induces on $k^2$ (corresponding to the three variables), resulting in a smoothly connected surface. Now the simplest example is $f=X^2Y$, yielding
$$F_A=X^2Y+Y^2Z+Z^2X,$$
$$F_S=X^2Y+X^2Z+Y^2X+Y^2Z+Z^2X+Z^2Y,$$
which apparently do the job.
