function bvp ( )

%*****************************************************************************80
%
%% bvp() is a sample two-point boundary value problem.
%
%  Reference:
%
%    Ward Cheney & David Kincaid
%    Numerical Mathematics and Computing, Fifth Edition
%    Brooks/Cole Publishing Company, 2003.
%
  solinit = bvpinit ( linspace(1,2,5), [0 1] );
  soln = bvp4c ( @bvpfcn, @bvpbc, solinit )
  x = linspace(1,2,12)
  y = deval(soln,x)
  yy = exp(x) - 3*cos(x)

  plot ( x, y(1,:), 'ro', x, yy(1,:), 'b-' )
  title('BVP problem solution')
  xlabel('x')
  ylabel('Solution y')
  grid ( 'on' );
  legend ( 'Computed', 'Exact' );
  filename = 'bvp_solution.png';
  print ( '-dpng', filename );
  fprintf ( 1, '  Graphics saved as "%s"\n', filename );

  return
end


