function mydx = diffloudApfiltA(t,x)

% Differential equation for device defined in loudApfiltA_com  

global P0 PP


k = fix(t*P0.fs)+1;
% Output of input filter
eP=P0.e;
uf = real(PP.s * eP(k,:) * (P0.r .* exp(P0.p *t)));


x01 = x(1);
x02 = x(2);
x03 = x(3);
x04 = x(4);
x05 = x(5);
x06 = x(6);
x07 = x(7);

mydx = zeros(3,1);

% Device
loudApfiltA_com
