JPD-Coloring of the Monohedral Tiling for the Plane

Box 3

> restart: with(linalg):

Probability_Values:=proc(n::posint,i::posint,j::posint,k::posint,l::posint,m::posint)

local A,b,S,r,v,C; global PV;

A[S]:=matrix(1,1,1,1,1,1): b:=vector(1): C:=sum(v[w],w=1..5):

linsolve(A[S],b,): r: linsolve(A[S],b,,v):

PV[S[n,i,j,k,l,m]]:=subs(v1=i/n,v2=j/n,v3=k/n,v4=l/n,v5=m/n,v6=1-C,%):

end proc:

n:=6: i:=1:j:=1:k:=1:l:=1:m:=1: PV[S[n,i,j,k,l,m]]:=Probability_Values(n,i,j,k,l,m);

n:=7: i:=1:j:=2:k:=1:l:=1:m:=1: PV[S[n,i,j,k,l,m]]:=Probability_Values(n,i,j,k,l,m);

n:=8: i:=1:j:=2:k:=1:l:=2:m:=1: PV[S[n,i,j,k,l,m]]:=Probability_Values(n,i,j,k,l,m);

n:=9: i:=1:j:=2:k:=1:l:=2:m:=1: PV[S[n,i,j,k,l,m]]:=Probability_Values(n,i,j,k,l,m);

n:=18:i:=1:j:=2:k:=3:l:=4:m:=6: PV[S[n,i,j,k,l,m]]:=Probability_Values(n,i,j,k,l,m);

n:=25:i:=1:j:=2:k:=3:l:=4:m:=5: PV[S[n,i,j,k,l,m]]:=Probability_Values(n,i,j,k,l,m);

n:=25:i:=1:j:=2:k:=3:l:=4:m:=5: PV[S[n,i,j,k,l,m]]:=Probability_Values(n,i,j,k,l,m);

, , ,

, , ,