for a = 1:500
	for b = 1:500
		c(a,b) = sin(a + b^2);
	endfor
endfor

g = eig(c);
max(real(g))


%Correct result is ans =  16.915
%With MKL in Ubuntu 20.04, I get random numbers of order 10^5 - 10^6, which changes every time