function skeleton6 ( )

%*****************************************************************************80
%
%% skeleton6 plots a stream function psi from a velocity data file.
%
%  Discussion:
%
%    A fluid flow problem was solved on a 41x41 grid.
%
%    The file 'xy_inout.txt' lists pairs of (x,y) values.
%
%    The file 'uv_inout.txt' lists pairs of (u,v) values.
%
%    The ordering of the nodes is suggested by this array:
%
%      1641 1642 1643 1644 ... 1680 1681
%      .... .... .... .... ... .... ....
%        83   84   85   85 ...  122  123
%        42   43   44   45 ...   81   82
%         1    2    3    4 ...   40   41
%
%    In other words, the first node has the lowest values of X and Y.
%
%    This program reads the data, computes the stream function, and plots it.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    25 June 2020
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'skeleton6\n' );
  fprintf ( 1, '  Read files of XY and UV components.\n' );
  fprintf ( 1, '  Compute the stream function and plot it.\n' );
%
%  Read the XY coordinates, two vectors of 41*41 pairs of (x,y) values.
%
  xy = load ( 'xy_inout.txt' );
  x = xy(:,1);
  y = xy(:,2);
%
%  Reshape the vectors into arrays.
%
  X = ?
  Y = ?
%
%  Read the UV values, two vectors of 41*41 pairs of (u,v) values.
%
  uv = load ( 'uv_inout.txt' );
  u = uv(:,1);
  v = uv(:,2);
%
%  Reshape the vectors into arrays.
%
  U = ?
  V = ?
%
%  Compute the stream function;
%
  PSI = ?
%
%  Plot the stream function.
%
  figure ( );
  contourf ( X, Y, PSI );
  xlabel ( '<-- x -->' );
  ylabel ( '<-- y -->' );
  title ( 'INOUT Stream function \psi(x,y)' );
  
  filename = 'hw8.png';
  print ( '-dpng', filename );
  fprintf ( 1, '  Graphics saved as "%s"\n', filename );

  return
end

