I thought for sure something like this would have already been done in Maple (mapleprimes).  I see it's already been done over at Mathematica.  I must admit I am a bit jealous at how much Mathematica users have done.  So anyways ... something simple

for i from 1 to 10 do
  for j from 1 to 10 do
     A[i,j]:=lcm(i,j):
  end do:
end do:

a:=convert(A,array):
b:=convert(a,Matrix):

I had to convert to array to get it into a Matrix, funny Maple won't convert directly to Matrix

plots[matrixplot](a,axes=boxed,shading=zhue,orientation=[-90,0];

Or without the grid


plots[matrixplot](aa, axes = boxed, shading = zhue, orientation = [-90, 0], style = patchnogrid)

And a 3d histogram

plots[matrixplot](aa, axes = boxed, shading = zhue, orientation = [-160, 57], heights = histogram, projection = .5)


Please Wait...