Let $x\in F\setminus\mathbb{R}1$. Since $F$ is finite-dimensional, say of dimension $N\geq1$, the family of $N+1$ vectors $(1,x,x^2,\ldots,x^N)$ is a dependent family. Hence there exists some real numbers $\alpha_0,\ldots,\alpha_N$ such that $(\alpha_0,\ldots,\alpha_N)\neq(0,\ldots,0)$ and such that
$$\alpha_N x^N+\cdots+\alpha_1x+\alpha_01=0.$$
Consider now the associated polynomial $P\in\mathbb{R}[X]$:
$$P=\alpha_N X^N+\cdots+\alpha_1X+\alpha_0.$$
Now, this polynomial can be factored as a product of irreducible polynomials in $\mathbb{R}[X]$; and by the Fundamental Theorem of Algebra we know that irreducible polynomials of $\mathbb{R}[X]$ are of degree $1$ or $2$. Without loss of generality we may hence assume that $P$ has the form:
$$P=X^2+aX+b.$$
(We can discard the case were the degree of $P$ is $1$, as this would mean that $x\in\mathbb{R}1$.)
Since $P$ is irreducible, we have $a^2/4-b>0$.
Now, define the element $j$ of $F$ as
$$j=\frac1{\sqrt{a^2/4-b}}x+\frac a{2\sqrt{a^2/4-b}}1.$$
We have:
$$j^2=\frac1{a^2/4-b}\left(x^2+ax+\frac{a^2}41\right)=\frac1{a^2/4-b}\left(-b+\frac{a^2}4\right)1=-1$$
since $P(x)=x^2+ax+b1=0$.
Observe that there are exactly two elements of $F$ that have a square equal to $-1$, namely $j$ and $-j$. Indeed, if $u\in F$, then the following factorization holds true:
$$u^2+1=(u+j)(u-j)$$
since $F$ is commutative, and since $F$ is an integral domain (since $F$ is a field),
$$u^2+1=0\iff u=-j\ \text{or}\ u=j.$$
We now conclude by showing that $F=\operatorname{Span}\{1,j\}=\mathbb{R}1\oplus\mathbb{R}j$. First observe that the family $(1,j)$ is independent since $j\not\in\mathbb{R}1$ since $j^2=-1$. Now, the previous discussion shows that given an element $y\in F\setminus\mathbb{R}1$, there exists an irreducible polynomial $Q\in\mathbb{R}[X]$ of degree $2$, say $Q=X^2+\alpha X+\beta$, such that $Q(y)=0$, and the computation we performed earlier shows that the element
$$u=\frac1{\sqrt{\alpha^2/4-b}}y+\frac\alpha{2\sqrt{\alpha^2/4-b}}1\in F$$
has a square equal to $-1$, hence this element is either $j$ or $-j$. Hence
$$y=-\frac\alpha2\,1\pm\sqrt{\alpha^2/4-b}\,j$$
and in any case belongs to $\operatorname{Span}\{1,j\}$. Hence $F\subset\operatorname{Span}\{1,j\}$ hence $F=\operatorname{Span}\{1,j\}$. It is now obvious that $F$ is isomorphic (field isomorphism) to $\mathbb{C}$.