function dxdt=fnode_ora(t,x)

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