Matlab Math Cleve Morler 著 陈文斌() 复旦大学2002年http:/www.stanford.edu/class/cs138/mlmath.html黄金分割 p = 1 1 1; r = roots(p); r = -0.61803398874989 1.61803398874989 C(1)*XN + . + C(N)*X + C(N+1) = 0 r = solve(1/x = x-1); r = 1/2*5(1/2)+1/2 1/2-1/2*5(1/2) 符号工具箱 Mapple Mapple phi = r(1) vpa (phi,50) 1/2*5(1/2)+1/2 1.6180339887498948482045868343656381177203091798058 phi = double(phi) 1.61803398874989图形的方法 f = inline(1/x-(x-1) ezplot(f,0,4) phi = fzero(f,1) hold on plot(phi,0,o)0 0.5 1 1.5 2 2.5 3 3.5 4 -3