[data,labels]=loadDigits(2000,'train')
a=data;
b=tangent_d(a,a);
options.dims=2;
[Y,R,E]=Isomap(b,'k',8,options);
j=0,k=0;
for i=1:2000
if(labels(i)==1)
j=j+1;
output5{1,j}=Y.coords{1,1}(1,i);output5{2,j}=Y.coords{1,1}(2,i);
elseif(labels(i)==7)
k=k+1;
output6{1,k}=Y.coords{1,1}(1,i);output6{2,k}=Y.coords{1,1}(2,i);
end
end
output5=cell2mat(output5);
output6=cell2mat(output6);

plot(output5(1,:),output5(2,:),'go');
hold on
plot(output6(1,:),output6(2,:),'bo');