a = 0.0;
b = 3.5;
x = linspace ( a, b, 101 );

clf ( );
hold ( 'on' );
grid ( 'on' );
for n = [ 1, 3, 5, 9, 17, 33 ]
  if ( n == 1 )
    xdata = ( a + b ) / 2.0;
  else
    xdata = linspace ( a, b, n );
  end
  w = omega ( xdata, x );
  plot ( x, w );
  label = sprintf ( 'omega(x) for n = %d', n );
  title ( label );
  pause ( 5 );
end
hold ( 'off' );

print ( '-dpng', 'omega_plot.png' );

