level 1
bixinyi7
楼主
b = 10.0; n = 2.0; m = 20;
a = {}; n0 = n; n1 = b;
Do[ak = 0; n2 = n0;
While[n2 >= n1, n2 = n2/n1; ak++],
AppendTo[a, ak]; n0 = n1; n1 = n2;
If[N[Abs[n2 - 1]] < N[10^(-m)], Break[]], {k, 1, m}];
len = Length[a];
y = a[[len]]; Do[y = a[[i]] + 1/y, {i, len - 1, 1, -1}];
Print[a];
Print[{m, y, N[y, m]}];
N[Log[10, 2], m]
Do::nliter: 位置 2 处的非列表迭代器(non-list iterator) AppendTo[a,ak];n0=n1;n1=n2;If[N[Abs[n2-1]]<N[10^-m],Break[]] 没有算出实数值.
2023年06月21日 01点06分
1
a = {}; n0 = n; n1 = b;
Do[ak = 0; n2 = n0;
While[n2 >= n1, n2 = n2/n1; ak++],
AppendTo[a, ak]; n0 = n1; n1 = n2;
If[N[Abs[n2 - 1]] < N[10^(-m)], Break[]], {k, 1, m}];
len = Length[a];
y = a[[len]]; Do[y = a[[i]] + 1/y, {i, len - 1, 1, -1}];
Print[a];
Print[{m, y, N[y, m]}];
N[Log[10, 2], m]
Do::nliter: 位置 2 处的非列表迭代器(non-list iterator) AppendTo[a,ak];n0=n1;n1=n2;If[N[Abs[n2-1]]<N[10^-m],Break[]] 没有算出实数值.
