Computational Physics, Python, Math, Physics.
u_n[i+1,:] + (CFL_3 - 1)/(CFL_3 + 1)*(u_np1[i+1,:] - u_n[i,:]) j = 0 u_np1[:,j] = u_n[:,j+1] + (CFL_1 - 1)/(CFL_1 + 1)*(u_np1[:,j+1] - u_n[:,j]) i = N_x u_np1[i,:] = u_n[i-1,:] + (CFL_4 - 1)/(CFL_4 + 1)*(u_np1[i-1,:] - u_n[i,:]) j = N_y u_np1[:,j] = u_n[:,j-1] + (CFL_2 - 1)/(CFL_2 +