function dxdt=fnode3_ora(t,x, hip)

dxdt=zeros(2,1);
if t<hip
    dxdt(1)=x(1)-x(1)^3/3-x(2)-0.6;
else
    dxdt(1)=x(1)-x(1)^3/3-x(2);
end
dxdt(2)=0.08*(x(1)+0.7-0.8*x(2));