[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)==4)
j=j+1;
output11{1,j}=Y.coords{1,1}(1,i);output11{2,j}=Y.coords{1,1}(2,i);
elseif(labels(i)==9)
k=k+1;
output12{1,k}=Y.coords{1,1}(1,i);output12{2,k}=Y.coords{1,1}(2,i);
end
end
output11=cell2mat(output11);
output12=cell2mat(output12);

plot(output11(1,:),output11(2,:),'go');
hold on
plot(output12(1,:),output12(2,:),'bo');