function value = omega ( x_data, x )

  n = length ( x_data );
  value = ones ( size ( x ) );

  for i = 1 : n
    value = value .* ( x - x_data(i) );
  end

  return
end

