function phi3_test ( )

%*****************************************************************************80
%
%% PHI3_TEST tests PHI3.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    04 June 2018
%
%  Author:
%
%    John Burkardt
%
%  Reference:
%
%    William Press, Brian Flannery, Saul Teukolsky, William Vetterling,
%    Numerical Recipes in FORTRAN: The Art of Scientific Computing,
%    Third Edition,
%    Cambridge University Press, 2007,
%    ISBN13: 978-0-521-88068-8,
%    LC: QA297.N866.
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'PHI3_TEST\n' );
  fprintf ( 1, '  PHI3 evaluates the thin plate spline radial basis function.\n' );

  nd = 5;

  r = [ 1; 2; 3; 4; 5 ];
  r8vec_print ( nd, r, '  R:' );

  r0 = 0.25;
  fprintf ( 1, '  R0 = %g\n', r0 );

  v = phi3 ( r, r0 );
  r8vec_print ( nd, v, '  PHI3(R,R0) = thin plate spline:' );

  return
end

