Te=
f0=
omega0=2*3.14*f0
Q=
K=
S=Table(0)
S=((2+(omega0*Te)/Q)*S[n-1]-S[n-2]+(K*Te*omega0)/Q*(EA0[n]-EA0[n-1]))/(1+(omega0*Te/Q)+omega0^2*Te^2)