fid = fopen('ans1.txt','w');
 fprintf(fid,'%6.2f \n',ans);