def jacobi$($A$,$ b$,$ p$0,$ tol, maxIter$=100)$:

n$=$len$($A$)$

p $=$ p$0$

for k in range$($maxIter$)$:

p$\_$old $=$ p$.$copy$()$ # In python assignment is not the same as

copyfor i in range$($n$)$:

sumi $=$ b$[$i$]$;

for j in range$($n$)$:

if i$==$j: # Diagonal elements are not included in

Jacobi

continue;

sumi $=$ sumi $-$ A$[$i$,$j$]*$ p$\_$old$[$j$]$

p$[$i$]=$ sumi$/$A$[$i$,$i$]$

rel$\_$error $=$ np$.$linalg.norm$($p$-$p$\_$old$)/$n # Actually $\text{'}$n$\text{'}$ should be

replace by norm of pif rel$\_$error