Computational Physics, Python, Math, Physics.
= V(X[i],Y[j]) U[:,:,0] = u_n.copy() #init cond - at t = 1 #without boundary cond u_np1[1:N_x,1:N_y] = 2*u_n[1:N_x,1:N_y] - (u_n[1:N_x,1:N_y] - 2*dt*V_init[1:N_x,1:N_y]) + Cx2*( 0.5*(q[1:N_x,1:N_y] + q[2:N_x+1,1:N_y ])*(u_n[2:N_x+1,1:N_y] - u_n[1:N_x,1:N_y]) - 0.5*(q[0:N_x -1,1:N_y] + q[1:N_x,1:N_y