(* ab4sigma.wl -- Algebraic Derivation of AB4 *)

rho=w^4-w^3
Solve[rho==0,w]
sigma=rho/Log[w]
T1=Series[sigma,{w,1,3}]
P1=Normal[T1]
Expand[P1]
Expand[rho]
