// Calculate the determinant of a 3 x 3 matrix
//
// David Edwards, Sun Dec 31 04:41:59 EST 2000

double det3(double **a)
{

    double d,d1,d2,d3;
    
    d1=a[1][1]*(a[2][2]*a[3][3] - a[3][2]*a[2][3]);
    d2=a[1][2]*(a[2][1]*a[3][3] - a[3][1]*a[2][3]);
    d3=a[1][3]*(a[2][1]*a[3][2] - a[2][2]*a[3][1]);
    d=d1 - d2 + d3;

    return d;
}