1.00001 + 1
10 ^ 309
10 ^ -309
(1 + 10^-16 - 1) * 10^16
(1 - 1 + 10^-16) * 10^16
(1:8) .* (2:9)
2 .^ (0:7)
round(cosd((0:7) * 90))
m = [1 2 3; 4 5 6; 7 8 9];
m * [0; 1; 0]
ei = [0; 1; 0]; ej = [0; 0; 1];
ei' * m * ej
[1 2; 3 4] .\ [2 4; 6 8]