Computational Physics, Python, Math, Physics.
= 0 u_jp1[0] = -u_jm1[0] + 2*u_j[0] + C2*( 0.5*(q[0] + q[0+1])*(u_j[0+1] - u_j[0]) - 0.5*(q[0] + q[0+1])*(u_j[0] - u_j[0+1])) elif left_bound_cond == 3: #Mur bound cond #i = 0 u_jp1[0] = u_j[1] + (CFL_1 -1)/(CFL_1 + 1)*( u_jp1[1] - u_j[0]) #right bound conditions if right_bound_cond == 1: