#hw26.txt Cole Franks 01/05/2014


Help:=proc() print(`ExpD2(M,N)`): end:


ExpD2:=proc(M,N) local eq,var,f,i,var1:



var:=[seq([seq(f[i,j],i=1..M-1)], j = 1..N-1)]:
eq:={seq(seq(f[i,j]=(1/4)*f[i+1,j]+(1/4)*f[i-1,j] +(1/4)*f[i, j-1] + (1/4)*f[i,j+1]+1,i=1..M-1),j=1..N-1 )}:
eq:=subs({seq(op({f[i,0]=0,f[i,N]=0}), i=1..M-1) ,seq(op({f[0,i]=0,f[M,i]=0}), i=0..N)},eq):
var1:=solve(eq,{seq(op(i), i in var)}):

subs(var1,var):





end: