

p = 1;


for i=2:size(PicVex,2)
	i
	A(:,i-1) = PicVex(:,i)/norm(PicVex(:,i));
end

p=2;
k=3;
l=4;
ll=5;

I = PicVex(:,p)' * A(:,k) * A(:,k) + PicVex(:,p)'*A(:,l)*A(:,l) + PicVex(:,p)'*A(:,ll)*A(:,ll);
