using namespace std;
    char** argv;

    Matrix mat;

    octave_value_list in, out;

    in(0) = "../../../matrices/10.mat";

    if( octave_main(1, argv, 1) ) {
        out = feval("PetscBinaryRead", in, 1);
        mat = out(0).matrix_value();

        clean_up_and_exit (0);
    } else {
        cout << "Octave was not started correctly." << endl;
    }