Here, $({X_n})_{n\geq1}$ is a sequence of random variables .
Result : $$ {X_n} \rightarrow X $$ if $E(X_n) \rightarrow E(X)$ and $ Var(X_n)\rightarrow Var(X)$ as $ n\rightarrow \infty $ .
Take any $ \varepsilon > 0.$
$$ \implies P[ | X_n-X| > \varepsilon ] \\ \\ \implies P [ | X_n -X|^2> \varepsilon^2 ] \leq \frac{E(X_n-X)^2}{\varepsilon^2} \\ \\ = \frac{[E(X_n)-E(X)]^2+ E[(X_n-E(X_n)]^2-E[(X-E(X)]^2}{\varepsilon^2 } . $$
I am not able to understand the last step .
I tried many ways but couldn't reach the above equation .
First, note that for any r.v.s $X,Y\in L_2$, $$ \mathsf{E}[X-Y]^2=(\mathsf{E}[X-Y])^2+\operatorname{Var}(X-Y). $$ Then the equality you're interested in is incorrect in general because $$ \operatorname{Var}(X-Y)\ne \operatorname{Var}(X)-\operatorname{Var}(Y) $$ unless $\operatorname{Var}(Y)=\operatorname{Cov}(X,Y)$. Take, for example, $X\sim N(0,1)$ and $Y=2X$. Then $$ 1=\operatorname{Var}(X-Y)\ne \operatorname{Var}(X)-\operatorname{Var}(Y)=-3. $$
Since $$ \mathsf{E}[X-X_n]^2=(\mathsf{E}[X_n-X])^2+\operatorname{Var}(X_n)+\operatorname{Var}(X)-2\operatorname{Cov}(X_n,X), $$ convergence in prob. follows if, in addition, $\operatorname{Cov}(X_n,X)\to \operatorname{Var}(X)$.