(* ab3sigma.wl -- Algebraic Derivation of AB3 *)

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