Re: newbie: programmatic sequence of plots?

*To*: mathgroup at smc.vnet.net*Subject*: [mg96817] Re: newbie: programmatic sequence of plots?*From*: Albert Retey <awnl at gmx-topmail.de>*Date*: Wed, 25 Feb 2009 04:07:03 -0500 (EST)*References*: <go0j4b$n0o$1@smc.vnet.net>

Tom Roche wrote: > How to get a programmatic sequence of plots? > > As advertised by the Subject:, I'm new to Mathematica (specifically > 7.0.0 for Students for 32-bit Windows on XP Pro) so please excuse any > lapses in terminology. I'd also appreciate pointers to specific docs > rather than just an RTFM. > > Since I'm currently unable to model the logistic map y=a y (1-y) (see > previous post) I'm trying to do a series of plots of it, to see how > behavior changes over time. I can get a plot of a single instance with > > Clear[a,y,y0,nRecurr]; > a=3.5; > y0=0.4; > nRecurr=50; > ListPlot[ > RecurrenceTable[{y[n+1]==a y[n] (1-y[n]),y[0]==y0}, > y,{n,1,nRecurr}],Joined->True > ] > > I can get a sequence of tables of values with > > Clear[a, y, y0, nRecurr, yinit, yinc, yfinal]; > a = 3.5; > nRecurr = 50; > yinit = 0.1; > yinc = 0.01; > yfinal = 0.2; > For[y0 = yinit, y0 <= yfinal, y0 += yinc, Print[y0, ":\n", > RecurrenceTable[{y[n + 1] == a y[n] (1 - y[n]), y[0] == y0}, y, {n, > 1, nRecurr}]]] > > But when I try to plot those values with > > Clear[a, y, y0, nRecurr, yinit, yinc, yfinal]; > a = 3.5; > nRecurr = 50; > yinit = 0.1; > yinc = 0.01; > yfinal = 0.2; > For[y0 = yinit, y0 <= yfinal, y0 += yinc, > ListPlot[ > RecurrenceTable[ > {y[n + 1] == a y[n] (1 - y[n]), y[0] == y0}, y, {n, 1, nRecurr} > ] > ,Joined->True] > ] > > I get absolutely nothing: no plot, no error message, no nada. > Similarly, when I try to plot the difference between two runs of the > map with > > seqDiff[a_,y_,inc_]:=Module[ > {r=y,s=y+inc}, > ListPlot[ > Table[r=a r (1-r),{50}]-Table[s=a s (1-s),{50}], > Joined->True,PlotRange->{-1,1} > ] > ] > seqDiff[4.0,0.3,0.1] > > I get a plot. But when I try to put that inside a for loop > > Clear[a, i,init,inc,final]; > a = 4.0; > init=0.3; > inc=0.1; > final=0.5; > For[i=init,i<=final,i+=inc,seqDiff[a,i,inc]] > > I get nothing: no plot, no error message. > > What am I doing wrong? How to get a series of plots? Just do what you did with the tables: Print the plots. You should understand that plots are no different than any other expression in Mathematica and only when returned as output, Mathematica knows how to format them. Since your For loop never returns anything, nothing is shown, just as you would expect for other expressions, too. Note that this has changed from version 5.2 to version 6. Another hint: Using Do instead of For leads to code which is easier to understand and less error prone. It also has the advantage that the syntax is exactly like that of Table which would return a list of these plots additionaly or instead of just printing them, so that you can use ListAnimation or TabView or whatever you like for displaying them. hth, albert