function force = Force(x,prm)

  % The derivative is different for negative and positive positions
  if x < 0
     force = -2*prm*(x+1);
  else
     force = -2*prm*(x-1);
  end
end

Published with MATLAB® R2012b