N=1000;
Mflop=2.0*N*N*(N/3.0+1)/1e6;
tmin=1e9;
for i=[1:10]
	A=2*rand(N,N)-1;
	x=ones(N,1);
	b=A*x;
	tic;
	y=A\b;
	t=toc;
	e=norm(y-1,1);
	fprintf("Residual %g, elapsed %g seconds, MFLOPS=%g\n",e,t,Mflop/t);
	if(t<tmin)
		tmin=t;
	end
end
fprintf("Maximum MFLOPS=%g\n",Mflop/tmin);
