20-Mar-2022 20:48:51

advection_pde_test()
  MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2
  Solve the advection PDE.

  Parameter values():
    c     = 5
    t0    = 0
    tstop = 0.08
    xmin  = 0
    xmax  = 1

advection_pde_ftcs():
  Solve the constant-velocity advection PDE in 1D,
    du/dt + c du/dx = 0
  over the interval:
    0.0 <= x <= 1.0
  with periodic boundary conditions:
    u(0) = u(1)
  and advection velocity
    c = 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: du/dx = (u(t,x+dx)-u(t,x-dx))/2/dx

  Number of nodes NX = 101
  Number of time steps NT = 8001
  Graphics saved as "advection_pde_ftcs_initial.png"
  Graphics saved as "advection_pde_ftcs_final.png"
  Movie saved as "advection_pde_ftcs.avi"
  Graphics saved as "advection_pde_ftcs_conservation.png"

advection_pde_test()
  Normal end of execution.
20-Mar-2022 20:51:27