Te=10e-7
ome=6000
Q=10
H0=1
S=table(0)
S[1]=0
S[2]=0
S[n]=2*S[n-1]-(1-ome^2*Te^2)*S[n-2]-Te/Q*ome*(S[n-1]-S[n-2])+ome*Te/Q*(EA0[n-1]-EA0[n-2])