I am reading about quaternions and their basic properties from this article Algebra and Geometry of Hamilton’s Quaternions
quaternions are quadruples of real numbers $q = (α, β, γ, δ)$ forming a four dimensional real vector space $H$ (named after Hamilton). In terms of the canonical basis: $e = (1, 0, 0, 0), i= (0, 1, 0, 0), j= (0, 0, 1, 0)$ and $k = (0, 0, 0, 1) (3)$ ....
Since any quaternion is a linear combination of the basis quaternions $e, i, j$ and $k$, it suffices to specify their products. First, $e$ is taken to be the multiplicative identity (sometimes denoted $1$), so $eq = qe = q$ for any quaternion $q$. In addition, $i^2 = j^2 = k^2 = −e$, while $ij = −ji = k$, $jk = −kj = i$ and $ki = −ik = j. (4)$
However I do not understand how to obtain the eqns $(4)$ using the forms of $e,i,j,k$ given in $(3)$? Any clue regarding the calculations would be helpful.