read LEON:
ope:=expand((1+N)*(1+M+M^2+M^3+M^4)-10):
print(`A basis for the linear vector space consisting of polynomials in m,n, of degree <=5`):
print(`viewed as functions on the discrete lattice with the curious property that for every 1 by 4 discrete`):
print(`rectangle in the plane, the value of the polynomial equals the average of its value at the other 9 locations is`):

print(CanBasPol(ope,M,N,m,n,5)):

print(`This took `, time(), `seconds `):

quit: