Computational Physics, Python, Math, Physics.
#Infinitesimal distance in the x direction N_y = int(L_y/dy) #Points number of the spatial mesh in the y direction Y = np.linspace(0,L_y,N_y+1) #Spatial array in the y direction #Temporal mesh with CFL < 1 - n indices L_t = 4 #Duration of simulation [s] dt = dt = 0.1*min(dx, dy) #Infinitesimal time