clc;
global p1 d1
x0 = [0 0.5]';
tol = 0.00000005;
x1 = mincg('mincg_f', 'mincg_pd', 'mincg_ftau', x0, tol);
disp('Press a key for plots')
pause
mincg_plot