%%
clear all
close all
clc

M = load_off('tr_reg_000.off');

plot_mesh(M);

%COLORMAP MODIFICATA

n = 100

red = [linspace(0, 1, n/2)', ones(1, n/2)'];
green = [linspace(0, 1, n/2)', linspace(1, 0, n/2)'];
blue = [ones(1, n/2)', linspace(1, 0, n/2)'];

bluewhitered = [red(:), green(:), blue(:)];

colormap (bluewhitered)

%DRAW

axis equal

shading interp


camlight head
lighting gouraud

drawnow


%CALCOLO DISTANZA EUCLIDEA

for i = 1:1:size(M.VERT, 1)
    d = mean(pdist2(M.VERT(i), M.VERT(:), 'euclidean'));

end