fid = fopen('setprob.data');
beta = fscanf(fid,'%g',1);   fscanf(fid,'%s',1);
freq = fscanf(fid,'%g',1);   fscanf(fid,'%s',1);
fclose(fid);