%
% read a matrix split into vectors
%
fil=input('Enter input file for sinc biharmonic:')
fid = fopen(fil);
t=[];
t=fscanf(fid, '%g %g %g %g %g', [5,inf]);
fclose(fid);
fil=input('Enter input file for sinc bispectral:')
fid = fopen(fil);
s=[];
s=fscanf(fid, '%g %g %g %g %g', [5,inf]);
fclose(fid);

hold off;

figure(1);
semilogy(t(1,:),t(2,:));
hold on;
semilogy(s(1,:),s(2,:),'r');
xlabel('N');
ylabel('Max Error');
%axis([0 1000 1e-15 1e6]);
hold off;

figure(2);
semilogy(t(1,:),t(3,:));
hold on;
semilogy(s(1,:),s(3,:),'r');
xlabel('N');
ylabel('L2 Error');
%axis([0 1000 1e-15 1e6]);
hold off;

figure(3);
semilogy(t(1,:),t(4,:));
hold on;
semilogy(s(1,:),s(4,:),'r');
xlabel('N');
ylabel('Rel Error');
%axis([0 1000 1e-15 1e20]);
hold off;

figure(4);
plot(t(1,:),t(5,:));
hold on;
plot(s(1,:),s(5,:),'r');
xlabel('N');
ylabel('Time');
%axis([0 1000 0 55]);
hold off;