octave利用黎曼和计算pi的近似值

黎曼和公式如下所示：

pi = lim(n->inf) sum(i=1 to n) 4/(1+(i-0.5/n)^2) / n

其中，n是区间的等分数量，i是区间的索引。在Octave中，我们可以使用循环结构来实现这个公式的计算。

以下是一个实现该公式的Octave代码示例：

function pi_approx = compute_pi(n)
  sum = 0;
  for i = 1:n
    x = (i - 0.5) / n;
    sum = sum + 4 / (1 + x^2);
  endfor
  pi_approx = sum / n;
endfunction

你可以调用这个函数，将等分数量n作为参数传递进去，即可获得pi的近似值。

例如，调用compute_pi(1000000)可以得到pi的近似值约为3.141593。