function dx = fnode2(t,x)

dx = zeros(2,1);
dx(1) = x(1) - x(1)^3/3-x(2);
dx(2) = 0.08*(x(1)+0.7-0.8*x(2));