def setup(self): self.integrator = ode(self.constitutive_relation) self.integrator.set_integrator("lsoda") # Main integrator function. Output vector: # 0: friction, 1+: state parameters