The specification for a0 and a1 can be input as solution commands in two forms
1. RAYLeigh,, a0, a1 where you compute the values
2. RAYLeigh FREQuency zeta omega1 omega2 where omega1 to omega2 is the range you want to approximate and zeta the damping ratio.
a1 = 2*zeta/omega1 + omega2
a0 = a1*omega1*omega2