function isaac_test ( )

%*****************************************************************************80
%
%% ISAAC_TEST does some RBF calculations.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    24 February 2015
%
%  Author:
%
%    John Burkardt
%
  timestamp ( )
  fprintf ( 1, '\n' );
  fprintf ( 1, 'RBF_ISAAC:\n' );
  fprintf ( 1, '  MATLAB version.\n' );
  fprintf ( 1, '  Do some sample tests and plots for Isaac Lyngaas.\n' );
%
%  Plots of one RBF, with several parameter values.
%
  one_rbf ( @ phi1, 'quadric' );
  one_rbf ( @ phi2, 'inverse_quadric' );
  one_rbf ( @ phi3, 'thin_plate_spline' );
  one_rbf ( @ phi4, 'gaussian' );
  pause
%
%  Plot two RBF's, with parameter 1.
%
  two_rbf ( @ phi4, 'gaussian' );
  pause
%
%  Plot three RBF's, with parameter 1, and symmetric spacing.
%
  three_rbf ( @ phi4, 'gaussian' );
  pause
%
%  Plot two RBF weight functions, as a sort of bar plot.
%
  two_weights ( @ phi4, 'gaussian' );
  pause
%
%  Plot three RBF weight functions, as a sort of bar plot.
%  I really need to make a sequence of plots as R0 varies!
%
  three_weights ( @ phi4, 'gaussian' );
  pause
%
%  Look at condition numbers for the 1D case, with 8, 16, 32, 64 nodes,
%  equally spaced, chebyshev spaced, randomly spaced, in [0,1].
%
  condition_test ( @ phi4, ');
%
%  Terminate.
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'RBF_ISAAC:\n' );
  fprintf ( 1, '  Normal end of execution.\n' );
  fprintf ( 1, '\n' );
  timestamp ( )

  return
end

