Re: coupled map lattice problem

*To*: mathgroup at smc.vnet.net*Subject*: [mg81337] Re: coupled map lattice problem*From*: B^3 <bbbld at rediffmail.com>*Date*: Wed, 19 Sep 2007 05:38:58 -0400 (EDT)*References*: <fcnlg4$sal$1@smc.vnet.net><fco7k2$3ts$1@smc.vnet.net>

On Sep 18, 3:00 pm, Jean-Marc Gulliet <jeanmarc.gull... at gmail.com> wrote: > B^3 wrote: > > <snip> > > > Please advice me.. > > You could start by writing some Mathematica code that you could then > post to this newsgroup (which is dedicated to the Computer Algebra > System (CAS) called Mathematica). > > -- > Jean-Marc I am trying something as per Mathematica Guide Book for Graphics but not sure is it the right way or not................ (time taken to run at even rang=5 is tooo long) Module[{rang = 500, k = 0.9, omega = 0.068, epsilon = 0.3, incr = 0.2}, step[l_, epsilon_, f_] := (1 - epsilon)f(l) + (epsilon/2) (f[RotateRight[l]] + f[RotateLeft[l]]); f = FractionalPart[# - k/(2Pi)Sin[2Pi#] + omega] &; lembda = Abs[Sin[7^(1/7)Range[rang]]]; tabul = Table[lembda = step[lembda, epsilon, f]; {#, j} & /@ Flatten[Position[Abs[Subtract @@@ Partition[lembda, 2, 1]], _?(# > incr &)]], {j, rang}]; Show[Graphics[Map[Rectangle[# - .5, # + .5] &, tabul, {2}]],AspectRatio -> Automatic]] Please advice......... Bhaskar Bidyut B.