使用數值方法求開普勒方程的近似解 R

m = 3; r = .5; result = c(m); i = 1; while(i < 30){ result[i + 1] = m + r * sin(result[i]); i = i + 1; } print(result); plot(1:length(result), result); 調整m與r的會發現,

在一定的範圍內,

m增大會減緩收斂速度

m減小會加快收斂速度;

r增大會加快收斂速度

r減小會減緩收斂速度