[data,labels]=loadDigits(2000,'train')
a=data;
b=L2_distance(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;
output1{1,j}=Y.coords{1,1}(1,i);output1{2,j}=Y.coords{1,1}(2,i);
elseif(labels(i)==7)
k=k+1;
output2{1,k}=Y.coords{1,1}(1,i);output2{2,k}=Y.coords{1,1}(2,i);
end
end
output1=cell2mat(output1);
output2=cell2mat(output2);

plot(output1(1,:),output1(2,:),'go');
hold on
plot(output2(1,:),output2(2,:),'bo');