% 
% Written by:
% -- 
% John L. Weatherwax                2005-08-14
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----

addpath( '../../Code/eda_data' );
addpath( '../../Code/eda_toolbox' );  

close all; drawnow; 
clc; 
clear; 

randn('seed',0); rand('seed',0); 

load scurve; % <- data is stored in the matrix "X" [3,2000]
[p,n] = size(X); 

xdist = pdist(X.','euclidean'); 
id    = idpettis(xdist,n); 
fprintf('the intrisic dimensionality of the scurve data is %10.6f\n', id ); 


load swissroll;  % <- data is stored in the matrix "X" [3,800]
[p,n] = size(X); 

xdist = pdist(X.','euclidean'); 
id    = idpettis(xdist,n); 
fprintf('the intrisic dimensionality of the swissroll data is %10.6f\n', id );