Computational Physics, Python, Math, Physics.
u_jp1[N_x] = u_j[N_x-1] + (CFL_2 -1)/(CFL_2 + 1)*(u_jp1[N_x-1] - u_j[N_x]) u_jm1[:] = u_j.copy() #go to the next step u_j[:] = u_jp1.copy() #go to the next step U[:,j] = u_j.copy() ############## PLOT ############### anim1 = viz_tools.anim_1D(X,U, dt, 10, save = False , myxlim = (0, 1.5) , myylim =