function value = wiggle ( x )
  value = cos ( 100.0 * x - 700.0 ) - 4.0 * erf ( 30.0 * x - 31.0 );
  return
end

