mathematica分形源码 - C

C 少壮不努力老二徒伤悲

http://jonahs.cgpad.com

空间首页

文章分类

最新评论

SPAN @ 2008-07-02

SPAN @ 2008-06-20

SPAN @ 2008-06-17

jonahs @ 2008-06-16

jonahs @ 2008-06-16

文章存档

2008年08月 (1)

2008年06月 (3)

mathematica分形源码

分类：发表：2008-06-23 19:44:06 点击：2689 回复：5

Fx[x_, y_, cx_, cy_, n_] := Block[{z, ct = 0}, z = x + y*I;
While[(Abs[z] < 2.0) && (ct < 50), ++ct; z = z^n + cx + cy*I];
Return[ct];]
Ht[cx_, cy_, n_, pu_List, po_List, pl_List] :=
Block[{kok},
kok = DensityPlot[
    Fx[xx, yy, cx, cy, n], {xx, pu[[2]], pu[[3]]}, {yy, po[[2]],
     po[[3]]}, pl, Mesh -> False]; Return[kok]]
Julia = Ht[0.54496, 0.45559,
3, {x, -1.5, 1.5}, {y, -1.5, 1.5}, {PlotPoints -> 200,
   PlotLabel -> "Julia.3"}]
Show[Julia,
Graphics[Line[{{0.4, -0.3}, {0.4, 0}, {0.8,
     0}, {0.8, -0.3}, {0.4, -0.3}}]]]

Fx[x_, y_, cx_, cy_, n_] := Block[{z, ct = 0}, z = x + y*I;
While[(Abs[z] < 2.0) && (ct < 50), ++ct; z = z^n + cx + cy*I];
Return[ct];]
Ht[cx_, cy_, n_, pu_List, po_List, pl_List] :=
Block[{kok},
kok = DensityPlot[
    Fx[xx, yy, cx, cy, n], {xx, pu[[2]], pu[[3]]}, {yy, po[[2]],
     po[[3]]}, pl, Mesh -> False]; Return[kok]]
Julia = Ht[0.69455, 0.28586,
4, {x, -1.5, 1.5}, {y, -1.5, 1.5}, {PlotPoints -> 200,
   PlotLabel -> "Julia.4"}]
Show[Julia,
Graphics[Line[{{0.4, -0.3}, {0.4, 0}, {0.8,
     0}, {0.8, -0.3}, {0.4, -0.3}}]]]

Fx[x_, y_, cx_, cy_, n_] := Block[{z, ct = 0}, z = x + y*I;
While[(Abs[z] < 2.0) && (ct < 50), ++ct; z = z^n + cx + cy*I];
Return[ct];]
Ht[cx_, cy_, n_, pu_List, po_List, pl_List] :=
Block[{kok},
kok = DensityPlot[
    Fx[xx, yy, cx, cy, n], {xx, pu[[2]], pu[[3]]}, {yy, po[[2]],
     po[[3]]}, pl, Mesh -> False]; Return[kok]]
Julia = Ht[0.340652, 0.7033651,
5, {x, -1.5, 1.5}, {y, -1.5, 1.5}, {PlotPoints -> 200,
   PlotLabel -> "Julia.5"}]
Show[Julia,
Graphics[Line[{{0.4, -0.3}, {0.4, 0}, {0.8,
     0}, {0.8, -0.3}, {0.4, -0.3}}]]]

Fx[x_, y_, cx_, cy_, n_] := Block[{z, ct = 0}, z = x + y*I;
While[(Abs[z] < 2.0) && (ct < 50), ++ct; z = z^n + cx + cy*I];
Return[ct];]
Ht[cx_, cy_, n_, pu_List, po_List, pl_List] :=
Block[{kok},
kok = DensityPlot[
    Fx[xx, yy, cx, cy, n], {xx, pu[[2]], pu[[3]]}, {yy, po[[2]],
     po[[3]]}, pl, Mesh -> False]; Return[kok]]
Julia = Ht[0.73251, 0.414193,
5, {x, -1.5, 1.5}, {y, -1.5, 1.5}, {PlotPoints -> 200,
   PlotLabel -> "Julia.6}]
Show[Julia,
Graphics[Line[{{0.4, -0.3}, {0.4, 0}, {0.8,
     0}, {0.8, -0.3}, {0.4, -0.3}}]]]

标签：

附件：

声明：CGPAD文章版权属于作者，受法律保护。没有作者书面许可不得转载。

评论或跟帖

#1 SPAN (张友邦) 2008-07-02 13:10:43

非常不错，谢谢分享。mathematic由于其数学基础的关系，非常适合做研究性指导的工作，但在产品级的渲染方面应该还是差一筹。

#2 views63 (Ze) 2008-12-01 12:36:00

网页上的代码最后一个漏了“，图片上的是正确的。

#3 xiu (秀) 2008-12-01 14:55:57

枫叶的眼睛真尖~~~！

#4 zhaowt (zhaowt) 2009-07-28 11:38:24

和matlab差不多了

#5 lapulandebaby (水果篮子) 2010-04-05 22:17:40

感谢分享啊，有学mathematica的加我QQ14458307，本人是初学者

请登陆后发表评论

搜索 CGPAD.COM