Computational Physics, Python, Math, Physics.
u_np1[-1,:] = 0 u_np1[:,0] = 0 u_np1[:,-1] = 0 elif bound_cond == 2: #Nuemann bound cond i,j = 0,0 u_np1[i,j] = 2*u_n[i,j] - (u_n[i,j] - 2*dt*V_init[i,j]) + Cx2*(q[i,j] + q[i+1,j])*(u_n[i+1,j] - u_n[i,j]) + Cy2*(q[i,j] + q[i,j+1])*(u_n[i,j+1] - u_n[i,j]) i,j = 0,N_y u_np1[i,j] = 2*u_n[i,j] -