program matmult
c  multiply 2 matrices together and display
c
	dimension a(3,3), b(3,3),c(3,3)
c
	write(*,*) 'Enter the first matrix:'
	write(*,*) 'Enter the first row of the first matrix:'
	read(*,*) a(1,1),a(1,2),a(1,3)
	write(*,*) 'Enter the second row of the first matrix:'
	read(*,*) a(2,1),a(2,2),a(2,3)
	write(*,*) 'Enter the third row of the first matrix:'
	read(*,*) a(3,1),a(3,2),a(3,3)
c
	write(*,*) 'Enter the second matrix:'
	write(*,*) 'Enter the first row of the second matrix:'
	read(*,*) b(1,1),b(1,2),b(1,3)
	write(*,*) 'Enter the second row of the second matrix:'
	read(*,*) b(2,1),b(2,2),b(2,3)
	write(*,*) 'Enter the third row of the second matrix:'
	read(*,*) b(3,1),b(3,2),b(3,3)
c
	write(*,*) 'Thank you!  Here is the product matrix:'
		write(*,*)
c
	do 10, i=1,3
		do 10, j=1,3
			c(i,j) = 0.
			do 10, k=1,3
				c(i,j)=c(i,j)+a(i,k)*b(k,j)
10				continue
c
	write(*,20) c(1,1),c(1,2),c(1,3)
	write(*,20) c(2,1),c(2,2),c(2,3)
	write(*,20) c(3,1),c(3,2),c(3,3)
20	format(3(2x,f8.4))
	call exit
	end