[data,labels]=loadDigits(3000,'train');
a=data;
b=tangent_d(a,a);
%b=L2_distance(a,a);
options.dims=1:10;
[Y,R,E]=Isomap(data,labels,b,'k',8,options);
