that multiplies matrix A by B: A(N1,N2) x B(N2,N3)
(2) Write a subroutine MATVEC(N,M,A,C)
that multiplies a matrix A(N,M) by the vector C(M)
(3) For examples that involve small matrices, check that this
works by writing a main program to call these subroutines, and
check the answers by hand.
(4)
(a) Successive random numbers can be generated using the built-in routine RAND
iseed = 123547
dum = rand(123547)
do i=1,10000
print*,rand(0)
end do
Use this approach to generate 10000 pairs of numbers X,Y randomly.
Plot these out to check they are uniformly distributed
(b)
Generate a 6x6 random matrix, invert it using the subroutine
MATINV and use matrix multiplication to confirm your
program is working.
(5) For a matrix A(N,N), write a program to determine the
transpose of the matrix, AT(N,N). Evaluate the product AT x A
and A x AT