diffusion_ftcs_pde:
  MATLAB/Octave version 9.9.0.1467703 (R2020b)

  Solve the diffusion PDE in 1D,
    du/dt - mu d2u/dx2 = 0
  over the interval:
    0.0 <= x <= 1.0
  with periodic boundary conditions:
    u(0) = u(1)
  and diffusion coefficient
    mu = constant
  and initial condition
    u(0,x) = (10x-6)^2 (8-10x)^2 for 0.6 <= x <= 0.8
           = 0 elsewhere.
  and NX equally spaced nodes in X,
  and NT equally spaced points in T,
  using the FTCS method:
    FT: Forward Time  : du/dt = (u(t+dt,x)-u(t,x))/dt
    CS: Centered Space: d2u/dx2 = (u(t,x+dx)-2u(t,x)+u(t,x-dx))/dx^2

  Number of nodes NX = 101
  Space ranges from 0 to 1
  Number of time steps NT = 201
  Time ranges from 0 to 0.01
  Constant diffusion coefficient MU = 0.5
  CFL coefficient ( must be < 0.5 ) = 0.25
  Graphics saved as "diffusion_ftcs_pde_initial.png"
  Graphics saved as "diffusion_ftcs_pde_final.png"
  Movie saved as "diffusion_ftcs_pde.avi"
  Graphics saved as "diffusion_ftcs_pde_conserved.png"

diffusion_ftcs_pde
  Normal end of execution.