import?numpy?as?np
import?matplotlib.pyplot?as?plt
a=np.zeros（（101101））
for?i?in?range（0101）:
????a[i][0]=0
????a[i][100]=0
for?j?in?range（1100）:
????a[100][j]=100
????a[0][j]=0
b=np.zeros（（101101））
k=0
e=1.0
t=0
while（e>0.001）:
????k=k+1
????print（k）
????e=0
????for?i?in?range（1100）:
????????for?j?in?range（1100）:
????????????b[i][j]=