Exemplo n.º 1
0
 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
Exemplo n.º 2
0
 def integrate(self, ts):
     rv = rk4int(self.ode, self.y, ts)
     self.y = rv[-1, :]
     return np.sin(rv)
Exemplo n.º 3
0
 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
Exemplo n.º 4
0
 def integrate(self,ts):
     rv = rk4int(self.ode,self.y,ts)
     self.y = rv[-1,:]
     return np.sin(rv)