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

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

