===================================== ODE: Double Pendulum ID: 3ca6ccb2a35117347611d11b3e55f310 ===================================== xd0 = x2 xd1 = x3 xd2 = -0.01*(200.0*par1*sin(x1)*x3*x2+100.0*par1*sin(x1)*(x3)**2-981.0*par1*sin(x0)-981.0*par0*sin(x0)+100.0*par1*cos(x1)*sin(x1)*(x2)**2+981.0*par1*sin(x0)*(cos(x1))**2+981.0*par1*cos(x1)*cos(x0)*sin(x1)+100.0*par1*sin(x1)*(x2)**2)/(-1.0*par1-1.0*par0+par1*(cos(x1))**2) xd3 = 0.01*(200.0*cos(x1)*par1*sin(x1)*x3*x2+100.0*cos(x1)*par1*sin(x1)*(x3)**2+981.0*par1*sin(x0)*(cos(x1))**2+981.0*par1*cos(x1)*cos(x0)*sin(x1)+200.0*par1*sin(x1)*x3*x2+100.0*par1*sin(x1)*(x3)**2-981.0*par1*sin(x0)-981.0*par0*sin(x0)+981.0*par1*cos(x0)*sin(x1)+200.0*par1*sin(x1)*(x2)**2+200.0*par1*cos(x1)*sin(x1)*(x2)**2+100.0*par0*sin(x1)*(x2)**2+981.0*par0*cos(x0)*sin(x1))/(-1.0*par1-1.0*par0+par1*(cos(x1))**2) ------------------------------------- IVP: 126 ------------------------------------- x0 = 2.355 x1 = -1.727 x2 = 0.43 x3 = 0.67 par0 = 1.0 par1 = 1.0 ------------------------------------- IVP: 127 ------------------------------------- x0 = [1.855,2.855] x1 = -1.727 x2 = 0.43 x3 = 0.67 par0 = 1.0 par1 = 1.0 ------------------------------------- IVP: 128 ------------------------------------- x0 = 2.355 x1 = -1.727 x2 = 0.43 x3 = 0.67 par0 = [0.5,1.5] par1 = 1.0 ------------------------------------- IVP: 129 ------------------------------------- x0 = [2.255,2.455] x1 = -1.727 x2 = 0.43 x3 = 0.67 par0 = [0.9,1.1] par1 = 1.0