#hw26.txt Andrew Lohr May 1

ExpD2:=proc(M,N) local eq,var,f,i,var1:
var:=[seq(seq(f[i][j],i=1..N-1),j=1..M-1)]:
eq:={seq(seq(f[i][j]=(f[i+1][j]+f[i-1][j]+f[i][j-1]+f[i][j+1])/4+1,i=1..N-1),j=1..M-1)}:
eq:=subs({seq(f[i][0],f[i][M],i=1..N-1),seq(f[0][j],f[N][j],j=1..M-1)},eq):
var1:=solve(eq,convert(var,set)):
subs(var1,var):
end: