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

plot(output3(1,:),output3(2,:),'go');
hold on
plot(output4(1,:),output4(2,:),'bo');