摘要
授權條款
我,本作品的著作權持有者,決定用以下授權條款發佈本作品:
- 您可以自由:
- 分享 – 複製、發佈和傳播本作品
- 重新修改 – 創作演繹作品
- 惟需遵照下列條件:
- 姓名標示 – 您必須指名出正確的製作者,和提供授權條款的連結,以及表示是否有對內容上做出變更。您可以用任何合理的方式來行動,但不得以任何方式表明授權條款是對您許可或是由您所使用。
- 相同方式分享 – 若要根據本素材進行再混合、轉換或創作,則必須以與原作相同或相容的授權來發布您的作品。
https://creativecommons.org/licenses/by-sa/3.0CC BY-SA 3.0 Creative Commons Attribution-Share Alike 3.0 truetrue
Maxima CAS src code
/*
golden ratio conjugate
= ((sqrt(5)-1)/2 = 0.618033988749895
It is approximated by finite continued fractions :
[0;1,1,1,....]
*/
kill(all);
iMax : 10;
/* continuead fraction - goldem mean */
f(i_Max):=
(
[a,i],
i:1,
a:[0,1,1],
while i<i_Max do
(a:endcons(1,a),
i:i+1),
float(cfdisrep(a))
)$
/* save the values to 2 lists */
xx:makelist (1, i, 1, 1); /* list of positive integers */
yy:makelist (f(1), i, 1, 1); /* list of cf */
for i:2 thru iMax step 1 do
(
xx:cons(i,xx),
y:float(f(i)),
yy:cons(y,yy)
);
load(draw);
draw2d(
file_name = "golden_mean",
terminal = 'png,
dimensions = [1000,1000],
title= "Finite continued fraction aproximation for golden mean",
key = "n-fraction",
xlabel = "n",
ylabel = "n-continued fractions",
point_type = filled_circle,
point_size = 1.0,
points_joined = true,
color = red,
points(xx,yy),
color = blue,
key = "golden mean",
explicit((sqrt(5)-1)/2,x,1,iMax)
);