136
136
// Now calculate the Hessian value and store it
137
137
real hessval = (fiplus + fiminus - 2*fi)/(2*_H*_H);
138
138
MatStore(hess, i, i, hessval, TRUE);
139
printf("Stored value %i,%i\n", i, i);
139
//printf("Stored value %i,%i\n", i, i);
140
140
printf("First derivative is: %.8f", (fiplus-fiminus)/(2*_H));