# Creating the simulation object
sim = ms.TimeDomainSimulation(D, Ze, pm=Pm, h0=S0, 
    fs=44100, piecewise_constant=False)
# changing integrator and solving 
sim.set_integrator('vode', nsteps=20000)
sim.integrate(t=tsim)