clear all
close all

filename = 'test25a';
%filename = 'WAMIT_Results\ESP_drifter';
%**************************************************************************

disp('Reading frequency domain hydrodynamic coefficients');
FD_coefs = ReadWamitFD(filename);

%Compute Impulse Response Functions from Frequency Domain Results
disp('Computing impulse response functions'); 
L = ComputeWamitIRF_coefs(filename);


dt = 0.1;       %timestep size
tf = 10;
N = tf/dt;      %Number of timesteps to compute

tau = -(N-1)*dt:dt:dt*(N-1);


i = 3;

L{7,i} = ppval(L{7,i},tau);

figure
plot(tau,L{7,3})
xlabel('s');