您还没有登陆到CGPAD.COM   现在登录   注册新用户
搜索 CGPAD.COM
 
 
 
 浏览论坛    搜索论坛    规章制度    帮助  
社区成员: 23570   主题总数: 1374   回复总数: 3307   帖子总数: 4681   论坛跳转:  
 
 
论坛 数学软件 Mathematica (一个广泛使用的数学相关软件) 本版斑竹:招募中...
 
   帖子列表  原创作品  
 
发表新主题
帖子搜索:  
 
mathematica分形源码
查看:2874  |  回复:5  |  创建:2008-06-23 19:44:06
 
jonahs (C)
注册: 2008-06-14
积分: 191 分
等级:
少壮不努力 老二徒伤悲

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

alt

 

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

alt

 

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

alt

 

 

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

alt





删除
 
编辑
 
标签:
 
附件:请登陆后查看附件内容!
 
声明:CGPAD文章版权属于作者,受法律保护。没有作者书面许可不得转载。
 
 <<上一页 1 下一页>>   

 注册: 2008-04-24
 积分: 12335 分
 等级:
 尘世如潮人如水 只叹江湖几人回


  2008-07-02 13:10:43 #1
非常不错,谢谢分享。mathematic由于其数学基础的关系,非常适合做研究性指导的工作,但在产品级的渲染方面应该还是差一筹。




 

 注册: 2008-12-01
 积分: 0 分
 等级: 等级修炼中
 您的签名信息


  2008-12-01 12:36:00 #2

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






 

 注册: 2008-05-06
 积分: 1399 分
 等级:
 


  2008-12-01 14:55:57 #3
枫叶的眼睛真尖~~~!




 

 注册: 2009-07-28
 积分: 0 分
 等级: 等级修炼中
 您的签名信息


  2009-07-28 11:38:24 #4
和matlab差不多了




 

 注册: 2010-04-05
 积分: 0 分
 等级: 等级修炼中
 物理爱好者


  2010-04-05 22:17:40 #5
感谢分享啊,有学mathematica的加我QQ14458307,本人是初学者




 
 <<上一页 1 下一页>>   
 
 
版权所有 © 2005-2008 CGPAD.COM,湘ICP备07500998号,兼容浏览器:IE6IE7FireFoxOperaSafariChrome
Total Requests: 5275983, Total Visits: 3085184, Processing Time: 37ms,