Computational Physics, Python, Math, Physics.
np.zeros((N_x+1,N_y+1),float) q = np.zeros((N_x+1, N_y+1), float) #init cond - at t = 0 for i in range(0, N_x+1): for j in range(0, N_y+1): q[i,j] = c[i,j]**2 for i in range(0, N_x+1): for j in range(0, N_y+1): u_n[i,j] = I(X[i],Y[j]) for i in range(0, N_x+1): for j in range(0, N_y+1): V_init[i,j]