COMPUTE FULL_RESPONSE_TRAJ

function [Aout, Zout,Znorm,varargout] = compute_full_response_traj(W_0, W_1, epsilon, t, p, om, outdof)

Full system response: $\mathbf{z}(t) = \mathbf{W}_{\epsilon}(\mathbf{p}(\Omega t), \Omega t)$

z = reduced_to_full_traj(t,p,W_0,W_1,epsilon,om);
%
% Timeseries at outdof
Zout = z(outdof,:);
%
% $$\|\mathbf{z}\|_{L_2}=\sqrt{\frac{1}{T}\int_0^T\mathbf\|\mathbf{z}(t)\|^2dt}$$
Znorm = sqrt(trapz(t, sum(z.^2))/(t(end)-t(1)));
% Aout = norm(Zout,'inf');
Aout = [];
for k=1:numel(outdof)
    amp = norm(Zout(k,:),'inf');
    Aout = [Aout amp];
end
varargout{1} = z(:,1);

end