Advection equation of a square wave, upwind method
Numerical solution of the advection equation u_t + a u_x = 0 with the upwind method for spatial discretization and the RK4 method for time integration.
Periodic boundary conditions.
Parameters:
Number of volumes= 200
a = 2
dx = 0.005
CFL number = 0.9