a bettter Siegel disk program in Mathematica

From: Roger Bagula (tftn_at_earthlink.net)
Date: 11/23/04


Date: Tue, 23 Nov 2004 22:35:18 GMT


I tried to make it work better.
You really have to watch getting the starting point too large!

Clear[x,y,a,b,s,g,a0,t]
(*Siegel Disk Julia map*)
gm=N[(1+Sqrt[5])/2];
a=Cos[2*Pi*gm];
b=Sin[2*Pi*gm];
digits=400;
x[n_,t_]:=x[n,t]=(x[n-1,t]*a-y[n-1,t]*b+x[n-1,t]^2-y[n-1,t]^2)
y[n_,t_]:=y[n,t]=(x[n-1,t]*b+y[n-1,t]*a+2*x[n-1,t]*y[n-1,t])
x[0,t_]:=0.50/(1+t/2);y[0,t_]=0.01/(1+t/2);
a=Flatten[Table[Table[{x[n,t],y[n,t]},{n,0, digits}],{t,1,10}],1];
ListPlot[a, PlotRange->All]
b=Table[a[[n]][[2]],{n,1,digits+1}];
ListPlot[b,PlotJoined->True, PlotRange->All]

Respectfully, Roger L. Bagula
tftn@earthlink.net, 11759Waterhill Road, Lakeside,Ca 92040-2905,tel: 619-5610814 :
alternative email: rlbtftn@netscape.net
URL : http://home.earthlink.net/~tftn




J_SIEGEL_DOUBLE.gif

Relevant Pages

  • a bettter Siegel disk program in Mathematica
    ... You really have to watch getting the starting point too large! ... (*Siegel Disk Julia map*) ... Respectfully, Roger L. Bagula ...
    (sci.fractals)
  • Roger DID NOT carry that 15 jacket onto court !!!
    ... "Roger Federer's cheeky little number hard to defend" ... It's truly amazing that columnist Giles Smith would write a full- ... Surely others here on rst also recorded it? ... Watch as an official comes up to him and Roger turns back and to his ...
    (rec.sport.tennis)
  • The Freaks on rogerfederer.com worse than those on rst??
    ... 1)"I feel really sad and just can't bring myself to watch the ... "Ah Roger, on to Basel! ... learned the true nature of his illness at the beginning of the year. ...
    (rec.sport.tennis)
  • Re: The Freaks on rogerfederer.com worse than those on rst??
    ... 1)"I feel really sad and just can't bring myself to watch the ... "Ah Roger, on to Basel! ... Any chance you can coax a few over to rst so I can administer some 'perspective'? ...
    (rec.sport.tennis)
  • Re: Roger Federer and his Amazing Balls
    ... Let us watch this rivalry. ... Roger Federer was hitting AmazeBalls. ... balls and he just finds these angles and blasted away. ... Did Max Mirnyi make the semis and his old doubles partner Roger Federer ...
    (rec.sport.tennis)