用python科学计算解偏微分方程代码

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]=