Python, Simulation, Numerical Analysis, Beginner, Fluid Simulation.
continued disturbance in center for 100 steps if (t < 100): u[t, Nx // 2, Ny // 2] = numpy.sin(t / 10) u[t+1, x, y] = c**2 * dt**2 * ( ((u[t, x+1, y] - 2*u[t, x, y] + u[t, x-1, y])/(dx**2)) + ((u[t, x, y+1] - 2*u[t, x, y] + u[t, x, y-1])/(dy**2)) ) + 2*u[t, x, y] - u[t-1, x, y] While this solution