14-May-2025 23:08:29

fd1d_advection_diffusion_steady_test():
  MATLAB/Octave version 6.4.0
  Test fd1d_advection_diffusion_steady().

FD1D_ADVECTION_DIFFUSION_STEADY:
  MATLAB version

  Solve the 1D steady advection diffusion equation:
    v du/dx - k d2u/dx2 = 0
  with constant, positive velocity V and diffusivity K,
  over the interval:
    0.0 <= x <= 1.0
  with boundary conditions:
    u(0) = 0, u(1) = 1.

  Use finite differences:
    d u/dx  = (u(x+dx)-u(x-dx))/2/dx
    d2u/dx2 = (u(x+dx)-2u(x)+u(x-dx))/dx^2

  Diffusivity K = 0.05
  Constant velocity V = 1
  Number of nodes NX = 41
  DX = 0.025
  Maximum safe DX is 0.1

fd1d_advection_diffusion_steady():
  Normal end of execution.
  Graphics saved as "fd1d_advection_diffusion_steady_test.png".

fd1d_advection_diffusion_steady_test():
  Normal end of execution.

14-May-2025 23:08:29