function cout=convpipe(t,n,dt,cin,d,l,q)
%
% vypocet odezvy cout(n) na vzruchovou funkci cin(n) pro ekvidistantni casy
% t(n).
% System je trubka o prumeru d, delky l, pri prutoku q.
% Predpoklada se konvektivni model.
%
tmean=pi*d^2*l/(4*q);
for i=1:n
    if t(i)<tmean/2
        e(i)=0;
    else
        e(i)=tmean^2/(2*t(i)^3);
    end
end
fcin=fft(cin);
fe=fft(e);
fcout=fcin.*fe;
cout=ifft(fcout)*dt;
end
