28-Jul-2021 00:21:39 diffusion_pde_test() MATLAB/Octave version 9.9.0.1467703 (R2020b) Solve the diffusion PDE. diffusion_pde_ftcs(): 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 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. 28-Jul-2021 00:22:33