(* 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]