%%HP: T(3)A(D)F(.);
[[ 'A*X^4+B*X^3+C*X^2+D*X+E' ]
 [ 'X=Y+T' ]
 [ 'Y^4+(4*A*T+B)/A*Y^3+(6*A*T^2+3*B*T+C)/A*Y^2+(4*A*T^3+3*B*T^2+2*C*T+D)/A*Y+(A*T^4+B*T^3+C*T^2+D*T+E)/A' ]
 [ 'T=-(B/(4*A))' ]
 [ 'Y^4+P*Y^2+Q*Y+R' ]
 [ 'P=(8*A*C-3*B^2)/(8*A^2)' ]
 [ 'Q=(8*A^2*D-4*A*B*C+B^3)/(8*A^3)' ]
 [ 'R=(256*A^3*E-64*A^2*B*D+16*A*B^2*C-3*B^4)/(256*A^4)' ]
 [ 'Y^4+P*Y^2+Q*Y+R=(Y^2+J*Y+K)*(Y^2+M*Y+N)' ]
 [ 'J+M=0' ]
 [ 'J=-M' ]
 [ 'P=K+N-M^2' ]
 [ 'Q=M*(K-N)' ]
 [ 'R=K*N' ]
 [ 'K^2+N^2+2*K*N=(P+M^2)^2' ]
 [ 'K^2+N^2-2*K*N=(Q/M)^2' ]
 [ '4*K*N=(4*R=(P+M^2)^2-(Q/M)^2)' ]
 [ 'M^6+2*P*M^4+(P^2-4*R)*M^2-Q^2=0' ]
 [ 'N=(P+M^2-Q/M)/2' ]]