07-Jan-2022 19:00:53

diffusion_pde_test()
  MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2
  Solve the diffusion PDE.

diffusion_pde_ftcs():
  MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2

  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

  Parameter values:
    mu    = 0.5
    t0    = 0
    tstop = 0.01
    xmin  = 0
    xmax  = 1

  Number of nodes NX = 101
  Number of time steps NT = 201
  CFL coefficient ( must be < 0.5 ) = 0.25
  Graphics saved as "diffusion_pde_ftcs_initial.png"
  Graphics saved as "diffusion_pde_ftcs_final.png"
  Movie saved as "diffusion_pde_ftcs.avi"
  Graphics saved as "diffusion_pde_ftcs_conservation.png"

diffusion_pde_ftcs():
  Normal end of execution.

diffusion_pde_test()
  Normal end of execution.
07-Jan-2022 19:01:45