Lindstedt–Poincaré method