%
% reads a epsilon vs. phi^2 file
%
fil=input('Enter input bihar file:')
fid = fopen(fil);
t=[];
t=fscanf(fid, '%g %g', [2,inf]);
fclose(fid);

fil=input('Enter input bispec file:')
fid = fopen(fil);
s=[];
s=fscanf(fid, '%g, %g', [2,inf]);
fclose(fid);

hold off;
plot(t(1,:),t(2,:),'b');
hold on;
xlabel('\epsilon');
ylabel('\Psi^2');
plot(s(1,:),s(2,:),'r');