JPD-Coloring of the Monohedral Tiling for the Plane

Box 2

> restart: with(linalg):

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

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

AF:=matrix(1,1,1,1): b:=vector(1): C:=sum(V[w],w=1..3):

linsolve(A[F],b,): : linsolve(A[F],b,,v):

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

end proc:

n:=4: i:=1:j:=1:k:=1: PV[F[n,i,j,k]]:=Probability_Values(n,i,j,k);

n:=5: i:=1:j:=2:k:=1: PV[F[n,i,j,k]]:=Probability_Values(n,i,j,k);

n:=6: i:=1:j:=2:k:=1: PV[F[n,i,j,k]]:=Probability_Values(n,i,j,k);

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

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

n:=8: i:=1:j:=2:k:=3: PV[F[n,i,j,k]]:=Probability_Values(n,i,j,k);

n:=8: i:=1:j:=3:k:=1: PV[F[n,i,j,k]]:=Probability_Values(n,i,j,k);

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

n:=9: i:=1:j:=4:k:=3: PV[F[n,i,j,k]]:=Probability_Values(n,i,j,k);

n:=13:i:=1:j:=2:k:=1: PV[F[n,i,j,k]]:=Probability_Values(n,i,j,k);

n:=13:i:=1:j:=3:k:=4: PV[F[n,i,j,k]]:=Probability_Values(n,i,j,k);

n:=13:i:=1:j:=1:k:=6: PV[F[n,i,j,k]]:=Probability_Values(n,i,j,k);

n:=13:i:=1:j:=1:k:=7: PV[F[n,i,j,k]]:=Probability_Values(n,i,j,k);

n:=13:i:=1:j:=1:k:=8: PV[F[n,i,j,k]]:=Probability_Values(n,i,j,k);

n:=14:i:=1:j:=1:k:=2: PV[F[n,i,j,k]]:=Probability_Values(n,i,j,k);

n:=14:i:=2:j:=3:k:=4: PV[F[n,i,j,k]]:=Probability_Values(n,i,j,k);

n:=14:i:=1:j:=5:k:=6: PV[F[n,i,j,k]]:=Probability_Values(n,i,j,k);

n:=14:i:=1:j:=2:k:=7: PV[F[n,i,j,k]]:=Probability_Values(n,i,j,k);

n:=14:i:=1:j:=2:k:=8: PV[F[n,i,j,k]]:=Probability_Values(n,i,j,k);

n:=14:i:=2:j:=2:k:=9: PV[F[n,i,j,k]]:=Probability_Values(n,i,j,k);

, , , ,

, , , ,

, , , ,

, , , ,

, , ,