clear;
clc;

% case 1
A = zeros(10,2);
G = A'*A;
[V,D] = eig(G);
fprintf ('\nConstant intensity: %f %f',D(2,2),D(1,1));

% case 2
A = randn(10,2)*0.1;
G = A'*A;
[V,D] = eig(G);
fprintf ('\nConstant intensity plus noise: %f %f',D(2,2),D(1,1));

% case 3
A = repmat([4 5],10,1);
G = A'*A;
[V,D] = eig(G);
fprintf ('\nConstant gradient: %f %f',D(2,2),D(1,1));

% case 4
A = repmat([4 5],10,1) + randn(10,2)*0.1; 
G = A'*A;
[V,D] = eig(G);
fprintf ('\nConstant gradient plus noise: %f %f',D(2,2),D(1,1));

% case 5
A = rand(10,2)*4;
G = A'*A;
[V,D] = eig(G);
fprintf ('\nGradients with several (more than 2) directions: %f %f',D(2,2),D(1,1));


