Whilst acknowledging the very nice solution by Yves Daoust, here's an alternative solution derived earlier using only basic calculus and symmetry of an ellipse about its centre. This is posted for general information only.
1. Centre of Ellipse
Equation of ellipse:
$$3x^2-x+6xy-3y+5y^2=0$$
Differentiating w.r.t $x$:
$$6x-1+6x\frac{dy}{dx}+6y-3y\frac{dy}{dx}+10y\frac{dy}{dx}=0$$
When $\dfrac{dy}{dx}=0$, $6x+10y=3$.
Substituting into the original equation gives $24x^2+16x-9=0$.
Solving gives
$$(x,y)=\left(-\frac13\pm\frac{\sqrt{70}}2,\frac12\mp\frac{70}2\right)$$
which are the coordinates of the highest and lowest points of the ellipse. Since both these points are symmetrical about the point
$$\left(-\frac13,\frac12\right)\qquad \blacksquare$$
this must be the centre of the ellipse, given that the ellipse is symmetrical about its centrepoint.
Translate the ellipse such that the centre is at the origin, by substituting $(x,y)$ in the original equation with $(x+h,y+k)$, where $(h,k)=(-\frac13,\frac12)$, which gives the translated ellipse equation as:
$$3x^2+6xy+5y^2-\frac{7}{12}=0$$
2. Incline of Ellipse
Consider the intersection of the translated ellipse with the line $y=mx$. The incline of (the major/minor axis of) the ellipse is the value of $m$ when the distance of the intersection points from the origin is at a maximum/minimum.
Substituting $y=mx$ into the original equation gives:
$$(3+6m+5m^2)x^2=\frac7{12}$$
Let $l$ be the distance of the intersection point from the origin.
$$l^2=x^2+y^2=x^2(1+m^2)=\frac7{12}\underbrace{\left(\frac{1+m^2}{3+6m+5m^2}\right)}_{f(m)}$$
Differentiating w.r.t. m shows that $\dfrac{df}{dm}=0$ when
$$3m^2-2m-3=0\\
m=\frac{1\pm\sqrt{10}}3\quad \Rightarrow m^2=\frac{11\pm2\sqrt{10}}9$$
Substituting into the above formula shows that $l^2$ (and hence $l$) is maximum when
$$m=\frac{1-\sqrt{10}}3\qquad \blacksquare$$
hence this is the incline (or rather, $\tan\theta$, where $\theta$ is the angle of incline) of the major axis of the ellipse.
3. Semi-major and Semi-minor Axes
The semi-major and semi-minor axes (usually denoted by $a,b$ in the standard ellipse notation $\frac{x^2}{a^2}+\frac{y^2}{b^2}=1$) are given by $l_{\text{max}},l_{\text{min}}$ respectively.
Hence, semi-major axis, $a$, is
$$a=l_{\text{max}}=\sqrt{\frac7{12}\left(\frac{1+m^2}{3+6m+5m^2}\right)}\biggr|_{m
=\frac{1-\sqrt{10}}3}
=\sqrt{\frac7{72}(4+\sqrt{10})}$$
and the semi-minor axis, $b$, is
$$b=l_{\text{min}}=\sqrt{\frac7{12}\left(\frac{1+m^2}{3+6m+5m^2}\right)}\biggr|_{m
=\frac{1+\sqrt{10}}3}
=\sqrt{\frac7{72}(4-\sqrt{10})}\qquad \blacksquare$$