close all
clear all

dt = .1;

FD_coefs = ReadData('tcyl3.out');

tau = 0:dt:10;
disp('Computing impulse response functions'); 
L = ComputeIRF(FD_coefs.w,FD_coefs.F,tau);


L_coefs = ComputeIRF_coefs(FD_coefs.w,FD_coefs.F,tau);
%N = L_coefs{1,1}.pieces+1;
%tau = -(N-1)*dt:dt:dt*(N-1);

for i = 1:6
    L_coefs{6,i} = ppval(L_coefs{6,i},tau+dt/2);
end

i = 1;
figure
plot(tau,L{6,i},tau+dt/2,L_coefs{6,i})