[data,labels]=loadDigits(2000,'train')
a=data;
b=tangent_d(a,a);
options.dims=2;
[Y,R,E]=Isomap(b,'k',8,options);
