![]() plot_surface ( Y, X, U_exact, cmap = cm. H is an iteration matrix that depends on A and B. Where x k + 1 and x k are approximations for the exact root of Ax B at (k + 1)th and kth iterations. The general iterative formulas can be given as: x k + 1 Hx k k 1, 2, 3. # enable 3D plotting import matplotlib.pyplot as plt from mpl_toolkits.mplot3d import Axes3D from matplotlib import cm from import spsolve u_exact = spsolve ( A, f ) print ( 'norm of the residual = grid points'. Iterative methods Jacobi and Gauss-Seidel in numerical analysis are based on the idea of successive approximations.
0 Comments
Leave a Reply. |
AuthorWrite something about yourself. No need to be fancy, just an overview. ArchivesCategories |