tic
[x]=lpfeasdir(A,b,c,x);
mytime = toc;
check = check_lp_opt(A,b,c,x);
check(7) = mytime;
disp(sprintf('ncase=%d, check = %f %f %f %f %f %f %f',ncase,check))
fprintf(fid,'ncase=%d, check = %f %f %f %f %f %f %f \n',ncase,check)