|
肃杀的冬日里,还未见一瓣雪花,圣诞就这样悄然降临。说不清是爱是恨,又哪管宗教国度,只知道活在当下,活在圣诞的幻梦中。喜悦感动或是寂寥孤独,都不要忘记《图形公式》就在这里。
昨夜寄给你的苹果可收到了么?
吃了?
没事儿,我在Mathmatica里键入ParametricPlot3D[{(1 + Cos[n]) Cos[m] + 0.07 Cos[5 m] + (0.96n/Pi)^95, (1 + Cos[n])*Sin[m], 5.7 Sin[n] + 2 Cos[n] - 0.73 Log[1 - n/Pi]}, {m,-Pi, Pi}, {n, -Pi, Pi}, Mesh -> False, BoxRatios -> 1,PlotStyle ->{Red, Specularity[Green, 10]}, ImageSize -> 400, Boxed -> False]让鲜亮不腐的苹果象征我们永恒的……
今夜,或许你会沉浸在《胡桃夹子》的旋律中睡去,愿你像Mary一样游一回雪国,去享受身在糖果王国的甜蜜,只是陪伴在你身边的不会是王子。但你在梦里也不会忘记,有一棵满缀幸福与欢乐的圣诞树,那是我用Mathmatica画出的s[prob_, A_,init_, max_] := FoldList[#2.{#[[1]], #[[2]], 1} &, init, RandomChoice[prob-> A, max]];
Christmas = {{{0.03, 0}, {0, 0.1}},{{0.848, 0}, {0, 0.848}}, {{0.8, 0}, {0, 0.8}}, {{0.2, -0.08}, {0.15, 0.22}},{{-0.2, 0.08}, {0.15, 0.22}}, {{0.25, -0.1}, {0.12, 0.25}}, {{-0.2,0.1}, {0.12,0.2}}};
tree = Map[List, {{0, 0}, {0, 1.5}, {0,1.5}, {0, 0.848}, {0, 0.848}, {0, 0.3}, {0, 0.4}}, {2}];
{A, prob, init, max} = {N@Join[Christmas,tree, 3], {2, 60, 10, 7, 7, 7, 7}/100., {0., 2.}, 10^5};
pts = ifs[prob, A, init, max]; //AbsoluteTiming
Graphics[{{Darker@Green, PointSize@Tiny,Point@pts}, {Hue@Random[], PointSize@Large, Point@#} & /@ RandomChoice[pts,200]}, AspectRatio -> 1.48]
每一颗光点都是对你的祝福。
圣诞了,Mathmatica不会言语,但他仍想画出苹果和圣诞树只为祝你平安喜乐。
关注“图形公式不烦恼”公众号,让我们用Mathmatica的神奇给你带来诗意与感动。
注:1)余王昕同学配文;2)程序修改自网络
拿起手机扫描二维码码,关注我的课程公众号吧!
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-24 12:44
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社