function u = diffusion_initial ( x )

%*****************************************************************************80
%
%% diffusion_initial() sets the initial condition for the diffusion PDE.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    26 January 2021
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    real X(1,NX), the coordinates of the nodes.
%
%  Output:
%
%    real U(1,NX), the initial condition for the solution.
%
  u = zeros ( size ( x ) );
  k = find ( 0.6 <= x & x <= 0.8 );
  u(k) = ( 10.0 * x(k) - 6.0 ).^2 .* ( 8.0 - 10.0 * x(k) ).^2;

  return
end
