reaction advection diffusion equation matlab code

{\Delta x^2} \] Advection term: Upwind or high-order schemes can be employed. Upwind schemes are often favored for stability: \[ \frac{\partial C}{\partial x} \approx \frac{C_i^n - C_{i-1}^n}{\Delta x} \] when flow is in the positive direction.