Explicit Derivation of Weierstrass Normal Form for Cubic Curve In page 22-23 of Rational Points on Elliptic Curves by Silverman and Tate, authors explain why is it possible to put every cubic curve into Weierstrass Normal Form. Here are relevant pages: (My question is at the end; I have put a red line across the point I am interested in)



Could someone explain to me how to "work out" the algebra part? 

This seems pretty important derivation actually. I am afraid I don't even know how to start the algebra. I would appreciate any suggestions/hints. Thank you :)
Related. There is a close question I found here on MSE. However, that particular question asks more about geometrical insight. My question is how does one explicitly go from the general cubic to the equation $xy^2+(ax+b)y=cx^2+dx+e$. By the way, the equation for the general cubic is (just so the coefficients are consistent with the above formula): $ax^3+bx^2y+cxy^2+dy^3+ex^2+fxy+gy^2+hx+iy+j=0$
 A: The calculation is straightforward using the Riemann-Roch theorem, but not so much without it. Let's look at the poles and zeroes of $x$ and $y$. Let $P$ be the third point of intersection of $Z$ and $C$, and $Q$ the third point of intersection of $X$ and $C$, and $A$ and $B$ be the two other points of intersection of $Y$ with $C$. Remark that:


*

*$Z$ has a double zero at $\mathcal O$ and a simple zero at $P$;

*$X$ has a double zero at $P$ and a simple zero at $Q$

*$Y$ has a simple zero at $\mathcal O$ and zeroes at $A$ and $B$.


From (1) and (2), $x=X/Z$ has a double pole at $\mathcal O$ (and simple zeros at $P$ and $Q$);
From (1) and (3), $y=Y/Z$ has simple poles at $\mathcal O$ and at $P$ (and simple zeros at $A$ and $B$).
It follows that the function $xy$ has a triple pole at $\mathcal O$ and no other poles. Thus, the functions $x$ and $y':=xy$ only have poles at $\mathcal O$ (of orders $2$ and $3$ respectively). Therefore, the 7 functions
$$1, x, x^2, x^3, y', (y')^2, xy'$$
all have poles of order $\leq 6$ at $\mathcal O$ and no other poles. But the vector space of such functions has dimension $6$ by the Riemann-Roch theorem. Therefore, we know that there is a linear combination of them which is $0$, i.e. a Weierstrass equation for $C$. 
In order to find the transformation explicitly, you can express the $7$ functions above as formal Laurent series around $\mathcal O$, with coefficients in $k(a,b,c,d,e,f,g,h,i,j)$, and find a combination which cancels the principal parts. Riemann-Roch is only required to prove that this will work, but if you want to avoid Riemann-Roch and you are extremely patient, it is possible to prove it directly. At some point in the calculation, you will have to divide by the discriminant $D(a,b,c,d,e,f,g,h,i,j)$... good luck with that! :)
A: By taking the derivative of the equation for the general cubic, and after some elementary manipulation, and preserving all of the original coefficients, one gets the following:
$x y^2 + [[(2b + 2c)/(3d + c)] x + (2g + f)/(3d +c)] xy$ = cubic in $x$
Now all one has to do is use a new coordinate $y$ that transforms the left hand side to $y^2$.  Silverman's change of variables, replacing $xy$ with a new $y$, does not seem to work for this form of the equation.
A: In the projective coordinates $[X, Y, Z]$, the intersection of the lines
$Y = 0$ and $Z = 0$ is $\mathcal{O} = [1, 0, 0]$, and the intersection of the
lines $X = 0$ and $Z = 0$ is $[0, 1, 0]$. Furthermore, $C$ and the line
$Z = 0$ have an intersection of multiplicity $2$ at $[1, 0, 0]$, and $C$ and
the line $X = 0$ have an intersection of multiplicity $2$ or $3$ at $[0, 1, 0]$.
In the affine space $X = 1$, $\mu_{[1, 0, 0]}(C, Z) = 2$ implies that the
equation for $C$ has, up to a constant factor, the form
$\alpha \cdot Y^3 + Y^2 + Z \cdot (\ldots) = 0$. Thus, in a homogeneous
equation for $C$, the coefficients of the monomials $X^3$ and $X^2 \cdot Y$
are zero, and the coefficient of $X \cdot Y^2$ is not zero.
In the affine space $Y = 1$, $2 \le \mu_{[0, 1, 0]}(C, X) \le 3$ implies that the
equation for $C$ has, up to a constant factor, the form
$\beta \cdot Z^3 + \gamma \cdot Z^2 + X \cdot (\ldots) = 0$, where $\beta$ and $\gamma$ are not both $0$.  Thus, in a homogeneous
equation for $C$, the coefficients of the monomials $Y^3$ and $Y^2 \cdot Z$
are zero, and the coefficients of $Y \cdot Z^2$ and of $Z^3$ are not both not zero.
Hence, in a homogeneous equation for $C$, only the following monomials
may have a non-zero coefficient: $X^2 \cdot Z$, $X \cdot Y^2$,
$X \cdot Y \cdot Z$, $X \cdot Z^2$, $Y \cdot Z^2$, and $Z^3$. Also, the
coefficient of $X \cdot Y^2$ is not zero, and the coefficients of $Y \cdot Z^2$ and $Z^3$ are not both zero. Therefore, in
the affine plane $Z = 1$, the equation for $C$ has, up to a constant factor,
the form:
$$ X \cdot Y^2 + (a \cdot X + b) \cdot Y = c \cdot X^2 + d \cdot X + e $$
where $b \ne 0$ or $e \ne 0$.
