I have been trying to calculate the Lyapunov exponent of the Rossler oscillator with an intention of finding it at higher precision. When i calculate the same at 16 digit precision on my 64 bit workstation or on my macbook using maple it takes much time (15 hrs have already been passed) and slows down the system badly. Can somebody plz help me out how to make my code more efficient and it runs in such a way that minimum resources of my computer are exploited.
The code basically has 6 differential equations to be integrated and then last three variables have to be stored as a vector whose log of norm has to be summed over large number of iterates. and now the elements of this normalized basis will be used as the initial condition for the last three equation. First three equations have to be integrated as usual.
Following is the code that i am using.........
for c from 8.7 to 9.0 by 0.01 do