fid=fopen('qmin.dat','a');

[xm,fm]=quadmin(2,4,1,5);
fprintf(fid, '%5.2f    %5.2f \n',xm,fm);

[xm,fm]=quadmin(2,4,-3,4);
fprintf(fid, '%5.2f    %5.2f \n',xm,fm);

[xm,fm]=quadmin(2,4,-10,-6);
fprintf(fid, '%5.2f    %5.2f \n',xm,fm);

fclose(fid);