JPD-Coloring of the Monohedral Tiling for the Plane

Box 1

> restart: with(linalg):

Probability_Values:=proc(n::posint,i::posint,j::posint)

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

AT:= matrix(1,1,1): b:= vector(1): C:=sum(Vw,w=1..2):

linsolve(AT,b,): : linsolve(AT,b,,v):

PVTn,i,j:=subs(v1=i/n,v2=j/n,v3=1-C,%):

end proc:

n:=3:i:=1:j:=1: PVTn,i,j:=Probability_Values(n,i,j);

n:=4:i:=1:j:=2: PVTn,i,j:=Probability_Values(n,i,j);

n:=5:i:=1:j:=2: PVTn,i,j:=Probability_Values(n,i,j);

n:=6:i:=1:j:=2: PVTn,i,j:=Probability_Values(n,i,j);

n:=7:i:=1:j:=2: PVTn,i,j:=Probability_Values(n,i,j);

n:=8:i:=1:j:=2: PVTn,i,j:=Probability_Values(n,i,j);

n:=8:i:=1:j:=3: PVTn,i,j:=Probability_Values(n,i,j);

, , , ,

, ,