{ diff(x1, t) = k3*(x1 + (-x1^3/3) + x2 + k4),
  diff(x2, t) = (-1/k3)*(x1 + (-k1) + k2*x2) }