function dx = fnode3(t,x,s)

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