       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>
>
>
> 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]]