Computational Physics, Python, Math, Physics.
print("Please choose a correct right boundary condition") #Spatial mesh - i indices L_x = 1.5 #Range of the domain according to x [m] dx = 0.01 #Infinitesimal distance N_x = int(L_x/dx) #Points number of the spatial mesh X = np.linspace(0,L_x,N_x+1) #Spatial array #Temporal mesh with CFL < 1 - j