# Load Packages #
with(LinearAlgebra):with(plots):

# Set up procedure for Digital Root #
DigRot := proc(n) description "base 10 digital root"; if irem(n,9) = 0 then return 9 else return irem(n,9); end if; end proc;

# Generate Vedic Square #
VS := Matrix(9): for i from 1 to 9 do for j from 1 to 9 do VS(i,j) := DigRot(i*j); end do: end do: VS;
# Create Matrix for highlighted plot
HighVS := proc(M::Matrix,n::integer) local i,j,tempM; tempM := Matrix(9): for i from 1 to 9 do for j from 1 to 9 do if M(i,j)=n then tempM(i,j):=0 else tempM(i,j) := 1 end if; end do: end do: tempM; end proc;

# Highlighted Plots #
for k from 1 to 9 do: k; matrixplot(HighVS(VS,k),heights=histogram,orientation=[0,0],axes=none,labels=["","",""],gap=0.2,shading=zgrayscale); end do;