This is a program discribing the dynamic behavior of a N pendelum system. When I tried to comparae the result of different initial conditions, I found that every time i run it, the graph will be different.
I know that double pendelum system is a typical chaotic system, is this the reason why i have different result for each time?
Here I attach the file, but plz rename it to ***.mws file before run, b/c i wrote it in maple classic worksheet.
with(SolveTools):with( plots ): with( plottools ):
with( DEtools ): with( PDEtools ):
#initialization function of each beam,in this function, we should give 3 variables, n=the number of this #beam, la=the length of this beam, ma=the mass of this beam,q=the initial angle of the beam,and qt=the initial angle velocity of the beam
l[n]:=la: #length of beam n
m[n]:=ma: #mass of beam n
#initialize the positions of each beam
x[n]:=0: #the x position of the end point of each beam
y[n]:=0: #the x position of the end point of each beam
x[n]:=0: #the y position of the center point of each beam
y[n]:=0: #the y position of the center point of each beam
#define the forces that applied on the end of beam n