概要
|
このファイルの ベクター画像 (SVG) が利用できます。 使う目的に対し、元画像よりもSVGがより優れている場合、SVG画像を使用して下さい。
|
|
Maxima source code
/* batch file for maxima
uses :
- symmetry around horizontal ( 0X ) axis
- Psi_M function to map conjugate plane to parameter plane
- jungreis algorithm to
time :
263 sec for jMax:100; and iMax:800;
4000 sec for jMax:400
*/
start:elapsed_run_time ();
jMax:300; /* precision = proportional to details and time of computations */
iMax:800; /* number of points to draw */
/* ------------ computes b coefficient of Jungreis function ----------------*/
betaF[n,m]:=block
(
[nnn:2^(n+1)-1],
if m=0
then 1.0
else if ((n>0) and (m < nnn))
then 0.0
else (betaF[n+1,m]- sum(betaF[n,k]*betaF[n,m-k],k,nnn,m-nnn)-betaF[0,m-nnn])/2.0
)$
b[m]:=betaF[0,m+1]$
/* -------------------------------*/
/* Power of w to j */
wn[w,j]:= if j=0 then 1 else w*wn[w,j-1]$
/* ---------Jungreis function ; c = Psi_M(w) ----------------------------- */
Psi_M(w):=w + sum(b[j]/wn[w,j],j,0,jMax)$
/* exponential for of complex number with angle in turns */
GiveCirclePoint(t):=R*%e^(%i*t*2*%pi)$ /* gives point of unit circle for angle t in turns */
/* -------------------------------------- w-plane -----------------------------------------*/
/*-------------- unit circle ------------*/
R:1;
circle_angles:makelist(i/iMax,i,0,iMax/2)$
CirclePoints:map(GiveCirclePoint,circle_angles)$
/* -------------------------parameter plane = c plane -----------------------------------*/
MPoints:map(Psi_M,CirclePoints); /* Mandelbrot set points */
/* add points below horizontal axis */
for w in CirclePoints do CirclePoints:cons(conjugate(w),CirclePoints);
for c in MPoints do MPoints:cons(conjugate(c),MPoints);
/* time */
stop:elapsed_run_time ();
time:fix(stop-start);
/* ---------------- draw *--------------------------------------------------------------------------*/
load(draw); /* Mario Rodríguez Riotorto http://www.telefonica.net/web2/biomates/maxima/gpdraw/index.html */
draw(file_name = "jung300",
dimensions=[1000, 500],
terminal = 'png,
columns = 2,
gr2d(title = " unit circle {w:abs(w)=1} ",
point_type = filled_circle,
points_joined =true,
point_size = 0.34,
color = red,
points(map(realpart, CirclePoints),map(imagpart, CirclePoints))
),
gr2d(title = " Image of circle under Psi_M(w) ",
points_joined =true,
point_type = filled_circle,
point_size =0.34,
color = blue,
points(map(realpart, MPoints),map(imagpart, MPoints))
)
);
References
- ↑ Drawing Mc by Jungreis Algorithm
ライセンス
この作品の著作権者である私は、この作品を以下のライセンスで提供します。
- あなたは以下の条件に従う場合に限り、自由に
- 共有 – 本作品を複製、頒布、展示、実演できます。
- 再構成 – 二次的著作物を作成できます。
- あなたの従うべき条件は以下の通りです。
- 表示 – あなたは適切なクレジットを表示し、ライセンスへのリンクを提供し、変更があったらその旨を示さなければなりません。これらは合理的であればどのような方法で行っても構いませんが、許諾者があなたやあなたの利用行為を支持していると示唆するような方法は除きます。
- 継承 – もしあなたがこの作品をリミックスしたり、改変したり、加工した場合には、あなたはあなたの貢献部分を元の作品とこれと同一または互換性があるライセンスの下に頒布しなければなりません。
https://creativecommons.org/licenses/by-sa/3.0CC BY-SA 3.0 Creative Commons Attribution-Share Alike 3.0 truetrue
|
この文書は、フリーソフトウェア財団発行のGNUフリー文書利用許諾書 (GNU Free Documentation License) 1.2またはそれ以降のバージョンの規約に基づき、複製や再配布、改変が許可されます。不可変更部分、表紙、背表紙はありません。このライセンスの複製は、GNUフリー文書利用許諾書という章に含まれています。http://www.gnu.org/copyleft/fdl.htmlGFDLGNU Free Documentation Licensetruetrue
|
あなたは上記のライセンスから、どれか一つ以上を選択できます。