function temperature_slice_y ( y )

%% temperature_slice_y plots t(x,y) for a fixed Y value.

  xvec = linspace ( 0.0, 6.0, 101 );
  yvec = y * ones ( 1, 101 );

  zvec = temperature ( xvec, yvec );

  plot ( xvec, zvec, 'Linewidth', 3 );
  xlabel ( '<-- X -->' );
  s = sprintf ( '<-- T(X,%f) -->', y );
  ylabel ( s );
  s = sprintf ( 'Temperature field T(X,%f)', y );
  title ( s );
  grid on
  print ( '-djpeg', 'temperature_slice_y.jpg' );

  return
end