Computational Physics, Physics, Math, Python.
everywhere, except along the boundary y=1 Hopefully you can see how the relaxed solution and this initial condition are related. Now, run the iterative solver with a target L2-norm difference between successive iterations of 10^−8 p = laplace2d(p.copy(), 1e-8) plot_3D(x,y,p) The Jacobi method is