Como calcular este problema com Runge Kuttade Segunda Ordem?

A equação de Van der Pol, da eletrônica, é:

y″ + (1 − y2)y′ + y = 0,

com as condições iniciais: y(0) = 0.5 e y′(0) = 0. Obtenha o valor de y, y′, y″ em t = 0.4 usando um método de Runge-Kutta de ordem 2.