PRO arr_matrix_multip

	a = indgen(3,2) + 1
	b = indgen(2,3) - 1
	print, 'a = '
	print, a
	print, 'b = '
	print, b
	
	c = a # b
	print, 'c = a # b'
	print, c
	
	print, 'transpose(c)'
	print, transpose(c)
	
	d = indgen(2,4) - 2
	e = indgen(4,2) + 2
	
	print, 'd = '
	print, d
	print, 'e = '
	print, e
	
	f = d ## e
	print, 'f = d ## e'
	print, f

END