; File: matrix.pro
; Author: Erik Brisson

A = dblarr(2, 4)
for i = 0, 1 do begin 
    for j = 0, 3 do begin 
        a(i, j) = 10 * i + j  
    endfor
endfor
print, 'A'
print, A

v1 = [1, 2, 0]
v2 = [1, 0, 0]
v3 = [4, 5, 6]
A = [[v1], [v2], [v3]]
print, 'A'
print, A

Atrans = transpose(A)
print, 'Atrans = transpose(A)'
print, Atrans

v = [1, 2, 3]
print, 'v'
print, v

w = A ## v
print, 'w = A ## v'
print, w

w = v ## A
print, 'w = v ## A'
print, w

d = determ(float(A))
print, 'd = determ(float(A))'
print, d

Ainv = invert(A)
print, 'Ainv = invert(A)'
print, Ainv

B = A ## Ainv
print, 'B = A ## Ainv'
print, B

b = float([1, 2, 16])
A = float(A)
x = cramer(A, b)
print, b
print, 'x = cramer(A, b)'
print, x

v4 = crossp(v1, v2)
print, "cross product crossp v4 = crossp,  v1, v2 = ", v4

end