Computational Physics, Python, Math, Physics.
indices L_t = 4 #Duration of simulation [s] dt = 0.01*dx #Infinitesimal time with CFL (Courant–Friedrichs–Lewy condition) N_t = int(L_t/dt) #Points number of the temporal mesh T = np.linspace(0,L_t,N_t+1) #Temporal array #Velocity array for calculation (finite elements) c = np.zeros(N_x+1, float)