Please i need the code in python thank you In lib490 complet

Please i need the code in python. thank you.

In lib490, complete the Lyapunov class in lyapunov. py. When complete, this class can be used with the ODE class to calculate the largest Lyapunov exponent of a continuous dynamical system. An important consideration is how this class will be used. The key idea is that a single ODE object will be used to evaluate two trajectories of the equations of motion at the same time. This is best illustrated with an example. The trajectories of two non-dimensionalized simple harmonic oscillators, A and B, can be combined into a single system. The list of dynamical variables is then y = [y0 y_1 y_2 y_3] = [x^A y^A x^B u^B] = [x^A x^A x^B x^B] The generalized velocities are then dy/dt = [dot x^A -x^A dot x^B -x^B] = [y_1 -y_0 y_3 -y_2]. With this approach the largest Lyapunov exponent will be calculated using: I import lib490.ode as ode ? import lib490.lyapunov as lyapunov PHYS 490 Assignment 5 1/3 Fall 2016 3 4 eom = lambda t, y: [y(1], -y(0], y(3), -y[2]] 5 y^0 = [1., 1., 1.+ le-8, 1.] 6 7 lyap = lyapunov.Lyapunov () 8 9 odesolver = ode. ODE (eom) 10 odesolver. Set_initial_value(y^0) 11 odesolver.set_integrator (ode.runge_kutta_2) 12 13 # warmup 14 nsteps = 1000 15 dt = .01 16 odesolver.integrate(nsteps, dt) 17 18 # production 19 nsteps = 1000 20 odesolver.set_solout (lyap.update_distance) 21 odesolver.integrate (nsteps, dt) 22 23 print (lyap. Largest_exponent ()) Complete the function update_distance(self, time, y) This will be used as a solout function in ODE and called once per time step. It should update the log of the total distance between the two integrations of the same system and rescale the separation between the two positions. The first half of y will be the main trajectory while the second half will be the perturbed trajectory. Complete the function largest_exponent(self) This function will return the largest Lyapunov. Complete the unit test test_three_stepsto show your class computes, within numerical tolerance, the total time, log of the total distance, and the largest Lyapunov exponent for the synthetic data provided.

Solution

Can help you with this any queries please comment

from math import log
def d(series,i,j):
return abs(series[i]-series[j])
f=open(\'timeseries.txt\', \'r\')
series=[float(i) for i in f.read().split()]
f.close()
N=len(series)
eps=input(\'Initial diameter bound: \')
dlist=[[] for i in range(N)]
n=0 #number of nearby pairs found
for i in range(N):
for j in range(i+1,N):
if d(series,i,j) < eps:
n+=1
print n
for k in range(min(N-i,N-j)):
dlist[k].append(log(d(series,i+k,j+k)))
f=open(\'lyapunov.txt\',\'w\')
for i in range(len(dlist)):
if len(dlist[i]):
print>>f, i, sum(dlist[i])/len(dlist[i])
f.close()