[zw024005]

S =
(-epsilon^4*sin(309/1250)+3*epsilon^2*sin(309/1250)+epsilon^4+epsilon^2)/(-sin(309/1250)+1-epsilon^2*sin(309/1250)+epsilon^2)
范围是-1到1
plot(epsilon,S)
谢谢高手帮忙！