def integrate(self, ts): #rv = eulerint(self.ode,self.y,ts,h=10**-2) rv = rk4int(self.ode, self.y, ts, h=1e-3) #rv = odeint(self.ode,self.y,ts) self.y = rv[-1, :] return rv
def integrate(self, ts): rv = rk4int(self.ode, self.y, ts) self.y = rv[-1, :] return np.sin(rv)
def integrate(self,ts): #rv = eulerint(self.ode,self.y,ts,h=10**-2) rv = rk4int(self.ode,self.y,ts,h=1e-3) #rv = odeint(self.ode,self.y,ts) self.y = rv[-1,:] return rv
def integrate(self,ts): rv = rk4int(self.ode,self.y,ts) self.y = rv[-1,:] return np.sin(rv)