function swrap(data, labels)

% size of data
[r, c] = size(data);

nearest_nbrs = 20;

% open output fid
fid = fopen('specbrain.txt', 'w');
fprintf(fid,'kNN\tFP\tFN\tFPr\tFNr\tTOT\tBASE\n');

output = [];
% iterate over different k-NN
for k=1:nearest_nbrs
    y = spectral(data,2,k);
    yinv(find(y==1)) = 2;
    yinv(find(y==2)) = 1;
    err = specerr(y, labels);
    errinv = specerr(yinv, labels);
    if (sum(errinv) < sum(err)),
        err = errinv;
    end;
    fprintf(fid,'%d\t%d\t%d\t%2.3f\t%2.3f\t%d\t4\n', k, err(1), err(2), err(1)/r, err(2)/r, err(1) + err(2));
    err = [k, err];
    output = [output; err];
    out = [y, labels];
    output
end;
fclose(fid);