#hw26.txt; Frank Wagner; Due 5/4/2014


###################
#####Problem 1#####
###################

#in progress....




###################
#####Problem 2#####
###################

ExpD2:=proc(M,N) local i,j,L,a,var,eq,var1:

for i from 0 to M do
 L[i,0]:=0:
od:

for j from 0 to N do
 L[0,j]:=0:
 L[M,j]:=0:
od:

for i from 1 to M-1 do
 L[i,1]:=a[i]:
od:

#L[i,j]=(1/4)*(L[i-1,j]+L[i+1,j]+L[i,j-1]+L[i,j+1])+1
#4*(L[i,j]-1)-L[i-1,j]-L[i+1,j]-L[i,j-1]=L[i,j+1]

for j from 2 to N do
 for i from 1 to M-1 do
  L[i,j]:=4*(L[i,j-1]-1)-L[i-1,j-1]-L[i+1,j-1]-L[i,j-2]:
 od:
od:

var:={seq(a[i],i=1..M-1)}:

eq:={seq(L[i,N],i=1..M-1)}:

var1:=solve(eq,var):

if var1=NULL then
 RETURN(FAIL):
fi:

subs(var1,[seq([seq(L[i,j],j=1..N-1)],i=1..M-1)]):


end: