function dx = fnode_parc2(t,x)

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

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

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

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