Is there any other solution to : $$\frac{\mathrm{d} g(x)}{\mathrm{d}x}=h(x)$$ $$\frac{\mathrm{d} h(x)}{\mathrm{d}x}=g(x)$$ other than $h(x)=g(x)=e^x$?
By varying $\alpha,\beta$ in
$$\frac{\mathrm{d} g(x)}{\mathrm{d}x}=\alpha h(x)$$ $$\frac{\mathrm{d} h(x)}{\mathrm{d}x}=\beta g(x)$$
is it possible to obtain $(e^x,e^x) , (\sin (x),\cos(x))$ as solutions when $\alpha = 1, \beta=1$ and $\alpha = 1, \beta=-1$ (without invoking complex analysis) is there any explanation for relationships between $\alpha,\beta$ yielding relationships between $e^x,\sin(x),\cos(x)$?
The system of differential equations can be written in matrix form as $$\frac{d\vec{u}}{dx}=A\vec{u},$$ where $$A=\begin{bmatrix}0&\alpha\\\beta&0\end{bmatrix}\text{ and } \vec{u}=\begin{bmatrix}g(x)\\h(x)\end{bmatrix}.$$ The general solution can then be expressed in terms of the eigenvalues and eigenvectors of $A$.
The eigenvalues are $\lambda_{1,2}=\pm\sqrt{\alpha\beta}$ with corresponding eigenvectors $$\vec{v}_{1,2}=\begin{bmatrix}-\alpha\\\pm\sqrt{\alpha\beta}\end{bmatrix}.$$
The general form of the solution is then given by $$ \vec{u}(x)=c_1\vec{v}_1e^{\lambda_1 x}+c_2\vec{v}_2e^{\lambda_2 x}.$$
In terms of non-exponential solutions, these can be obtained with complex eigenvalues, i.e. $\alpha\beta<0$.