发送短信 : A numerical solution for Duffing-Van der Pol oscillators using a backward difference formulation