Function MatInv_2

src/CPnumerics.cpp:71–80

}

void MatInv_2(double A[2][2], double B[2][2]) {
    double Det = NAN;
    //Using Cramer's Rule to solve

    Det = A[0][0] * A[1][1] - A[1][0] * A[0][1];
    B[0][0] = 1.0 / Det * A[1][1];
    B[1][1] = 1.0 / Det * A[0][0];
    B[1][0] = -1.0 / Det * A[1][0];
    B[0][1] = -1.0 / Det * A[0][1];
}

void solve_cubic(double a, double b, double c, double d, int& N, double& x0, double& x1, double& x2) {
    // 0 = ax^3 + b*x^2 + c*x + d