function dxdt = fnode_parc( t,x)

dxdt=zeros(length(x),1);
d = 1;

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


for i = 1:(length(x)/2-2)
    dxdt(2*i+1) = (x(2*i-1)-2*x(2*i+1)+x(2*i+3))/d^2 +x(2*i+1)-(x(2*i+1)^3)/3-x(2*i+2);
    dxdt(2*i+2) = 0.08*(x(2*i+1)+0.7-0.8*x(2*i+2));
end

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

