Mathematica 9 is now available
Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
2006
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*November
*December
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 2006

[Date Index] [Thread Index] [Author Index]

Search the Archive

Another older one that is a treat : a triangular von Koch type fractal

  • To: mathgroup at smc.vnet.net
  • Subject: [mg65308] Another older one that is a treat : a triangular von Koch type fractal
  • From: Roger Bagula <rlbagulatftn at yahoo.com>
  • Date: Sat, 25 Mar 2006 05:17:44 -0500 (EST)
  • References: <e002hj$pp1$1@smc.vnet.net>
  • Sender: owner-wri-mathgroup at wolfram.com

My aspect comes out two to one so that when roated 90 degrees it is 
distorted, but this picture comes out fine:

(* Cartoon/ von Koch as Peak : Besicovitch - Ursell function*)
f[x_] := 0 /; 0 <= x <= 1/3
f[x_] := -2 + 6*x /; 1/3 < x <= 1/2
f[x_] := 4 - 6*x /; 1/2 < x <= 2/3
f[x_] := 0 /; 2/3 < x <= 1
ff[x_] := f[Mod[Abs[x], 1]]
Plot[f[Mod[Abs[x], 1]], {x, 0, 2}]
s0 = Log[2]/Log[3]
(* Cartoon/ as Sigmoid : Besicovitch - Ursell function*)
g[x_] := 0 /; 0 <= x <= 1/3
g[x_] := (3*x - 1) /; 1/3 < x <= 2/3
g[x_] := 1 /; 2/3 < x <= 1
gg[x_] := g[Mod[Abs[x], 1]]
ParametricPlot[{f[t], gg[t]}, {t, 0, 1}, Axes -> False]
Plot[gg[t], {t, 0, 2}]
hh[x_] = Sum[gg[3^k*x]/3^(s0*k), {k, 0, 20}];
kk[x_] = Sum[ff[3^k*(x)]/3^(s0*k), {k, 0, 20}];
a = Table[{kk[n/30000], hh[n/30000]}, {n, 1, 30000}];
ga = Show[Graphics[{PointSize[0.003], Point /@ a}], Axes -> False]
> 
> 
> Roger L. Bagula { email: rlbagula at sbcglobal.net or rlbagulatftn at yahoo.com }
> 
> 11759 Waterhill Road,
> Lakeside, Ca. 92040    telephone: 619-561-0814
> 


  • Prev by Date: Re: Listing the partitions of a set
  • Next by Date: Re: Listing the partitions of a set
  • Previous by thread: constraints and NMinimize
  • Next by thread: Re: infinite product