我尝试使用以下代码编写用于近似pi的Newton / Euler近似程序，但未给出所需的结果。我是递归的新手，所以这可能是问题所在。方法是pi / 2 = 1 + 1/3（1 + 2/5（1 + 3/7（1 + ......
代码：
function [pi] = approx_pi_newtoneuler(N)
if N==1
piq = 1;
else
piq = 1 + (N / (2 * N + 1)) * approx_pi_newtoneuler(N-1);
end
pi = piq * 2;
end