(* newcubic.wl -- New Multistep Method with Cubic Order *)

rho=(w-1)*(w-1/2)*(w+1/2)
Solve[rho==0,w]
sigma=rho/Log[w]
T1=Series[sigma,{w,1,2}]
P1=Normal[T1]
Expand[P1]
Expand[rho]
