Computational Physics, Python, Math, Physics.
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,N_y+1), float) for i in range(0,N_x+1): for j in range(0,N_y+1): c[i,j] =