int dOF = 459; /* Degrees of freedom also equation to number of equations */
MKL_INT n, m;
n = dOF;
m = n;
int c, k, i, j, rows, cols;
vector> matJacbn(dOF, vector(dOF, 0));
double* x;
/* memory allocation */
x = (double*) malloc (sizeof (double)*n);
fjac = (double*) malloc (sizeof (double)*m*n);
/*
ofstream out1("fjac-array.txt");
for (cols=0; cols