Computational Physics, Python, Math, Physics.
for i in range(0,N_x+1): c[i] = celer(X[i]) ############## CALCULATION CONSTANTS ############### c_1 = c[0] c_2 = c[N_x] C2 = (dt/dx)**2 CFL_1 = c_1*(dt/dx) CFL_2 = c_2*(dt/dx) ############## PROCESSING LOOP ############### if loop_exec: # $\forall i \in {0,...,N_x}$ u_jm1 = np.zeros(N_x+1,float)