inverse of 3X3 matrix

[praneeth]

hi i want a c program to find the inverse of 3X3 matrix
urgent
i need it
fast
anyoone reply
thanks

[Banfa]

Here is some visual basic code that does something similar, see if you can convert and adapt it.

Can you describe the process of inverting a matrix?

Private Function InternalMatrixSolve(ByRef r1() As Double, ByRef r2() As Double, ByVal size%)
    Dim Matrix(2, 2) As Double
    Dim Vector(2) As Double
 
    Rem Set up the normalised matrix and vector
 
    For y% = 0 To 2
        For x% = 0 To 2
            Matrix(y%, x%) = 0
            For r% = 1 To size
                Matrix(y%, x%) = Matrix(y%, x%) + ((r1(r%) ^ y%) * (r1(r%) ^ x%))
            Next
        Next
 
        Vector(y%) = 0
        For r% = 1 To size
            Vector(y%) = Vector(y%) + ((r1(r%) ^ y%) * r2(r%))
        Next
    Next
 
    For y% = 0 To 2
 
        Rem Get a 1 on the diagonal by dividing through the current line by the diagnoal cell
        Divisor# = Matrix(y%, y%)
 
        If Divisor# = 0 Then
            InternalMatrixSolve = "MATHS ERROR"
            Exit Function
        End If
 
        For x% = y% To 2
            Matrix(y%, x%) = Matrix(y%, x%) / Divisor#
        Next
 
        Vector(y%) = Vector(y%) / Divisor#
 
        Rem get 0s below the current diagonal by taking away multiples of the current line
        For y2% = y% + 1 To 2
            Multiplier# = Matrix(y2%, y%)
            For x% = y% To 2
                Matrix(y2%, x%) = Matrix(y2%, x%) - (Matrix(y%, x%) * Multiplier#)
            Next
 
            Vector(y2%) = Vector(y2%) - (Vector(y%) * Multiplier#)
        Next
    Next
 
    Rem we now have 1s on the diagnoal and )s below the diagonal
    Rem need to get 0s above the diagonal
 
    For y% = 2 To 0 Step -1
 
        Rem get 0s above the current diagonal by taking away multiples of the current line
        For y2% = y% - 1 To 0 Step -1
            Multiplier# = Matrix(y2%, y%)
            For x% = y% To 2
                Matrix(y2%, x%) = Matrix(y2%, x%) - (Matrix(y%, x%) * Multiplier#)
            Next
 
            Vector(y2%) = Vector(y2%) - (Vector(y%) * Multiplier#)
        Next
    Next
 
    InternalMatrixSolve = Vector
End Function