function?GaussSeidel（abpdeltamax）
n=length（b）;
p1=p;
for?k=1:max
????for?i=1:n
????????x（i）=（b（i）-a（i[1:i-1]）*p（1:i-1）-a（i[i+1: