open Impact

let c = [| [| 2.; 1.; -1.; 8.|];
           [| -3.; -1.; 2.; -11.|]; 
           [| -2.; 1.; 2.; -3.|] |] 

let d = [| [| 4.; 8.; 4.|]; 
           [| -3.; -1.; 7.|] |]

let printMatrix(matrix : float[][]) =
    for rowNum in 0 .. matrix.Length-1 do
        let row = matrix.[rowNum]
        for colNum in 0 .. row.Length-1 do
            printf "  %10f" row.[colNum]
        printfn ""
    printfn ""
 
let scalarMult (vector : float [], multiplier : float) =
    [| for i in 0 .. vector.Length-1 -> vector.[i] * multiplier |]
 
let vectorAdd(v1 : float [], v2 : float []) =
    [| for i in 0 .. v1.Length-1 -> v1.[i] + v2.[i]|]

c.[0] <- scalarMult(c.[0],1.0/2.0); printMatrix c
c.[1] <- vectorAdd(scalarMult(c.[0], 3.),c.[1]); printMatrix c
c.[2] <- vectorAdd(scalarMult(c.[0], 2.),c.[2]); printMatrix c