
clear;
close all;
load clustdata.mat;

r=4;
beta=0.1;

X=X1;
for k=2:5
  y = spectral(X,k,r,beta);
  refresh;
  eval(sprintf('print -depsc hw5prob3_7_X1_%d.eps',k));
end

X=X2;
for k=2:5
  y = spectral(X,k,r,beta);
  refresh;
  eval(sprintf('print -depsc hw5prob3_7_X2_%d.eps',k));
end

X=X3;
for k=2:5
  y = spectral(X,k,r,beta);
  refresh;
  eval(sprintf('print -depsc hw5prob3_7_X3_%d.eps',k));
end
