1

我想在 org-mode 文件中嵌入简单的数学表达式，例如：

sqrt(p * (1-p) * N)

其中 p=0.8, N=10,000 并且根据简单的桌面计算器的答案是 40。我没有运气让 org-babel 给我这个答案...这是一个 org 文件：

Some weird rounding errors with Python:

#+begin_src python :var N=100000 :var p=0.8
from math import sqrt
return sqrt(p * (1 - p) * N)
#+end_src

#+results:
: 126.491106407


Some other weird rounding errors, this time with Emacs lisp:

#+begin_src elisp :var N=100000 :var p=0.8
(sqrt (* p (- 1 p) N))
#+end_src

#+results:
: 126.49110640673517


Note sure which number calc is unhappy with:

#+begin_src calc :var N=100000 :var p=0.8
sqrt(p * (1 - p) * N)
#+end_src

#+results:
| 5 | Expected a number |


Yet more weird rounding errors, using bc in the shell:

#+begin_src sh :var N=100000 :var p=0.8
echo "sqrt($p * (1 - $p) * $N)" | bc
#+end_src

#+results:
: 100.0


A different set of rounding errors with bc without using variables:

#+begin_src sh
echo "sqrt(0.8 * (1 - 0.8) * 10000)" | bc
#+end_src

#+results:
: 31.6


Looks like the variables are being substituted OK:

#+begin_src sh :var N=100000 :var p=0.8
echo "sqrt($p * (1 - $p) * $N)"
#+end_src

#+results:
: sqrt(0.8 * (1 - 0.8) * 100000)


Using extra precision directly gives the correct result:

#+begin_src sh
echo "sqrt(0.80 * (1 - 0.80) * 10000)" | bc
#+end_src

#+results:
: 40.0


So, just use extra precision in the p variable, eh:

#+begin_src sh :var N=100000 :var p=0.80
echo "sqrt($p * (1 - $p) * $N)" | bc
#+end_src

#+results:
: 100.0


No! Because it get stripped out:

#+begin_src sh :var N=100000 :var p=0.80
echo "sqrt($p * (1 - $p) * $N)"
#+end_src

#+results:
: sqrt(0.8 * (1 - 0.8) * 100000)


Please God, kill me now...

我只是在寻找一种方便的方法来在 org 文件中嵌入一些数学。

（来自最近的 git checkout 的 org-mode 7.8.03，来自最近的 git chekout 的 emacs 24.0.93.1）

4

2 回答 2

4

不要因为错字而自杀：）

(sqrt (* p (- 1 p)))
0.39999999999999997

添加缺少的 N..

(let ((p 0.8) (N 10000))
  (sqrt (* p (- 1 p) N)))

40.0

对于 bc，尝试使用

scale=10
于 2012-03-04T22:05:20.840 回答
2

你说：

sqrt(p * (1-p) * N)

where p=0.8, N=10,000 and the answer according to a simple desk calculator is 40.

这是正确的，但在您使用的许多代码中

:var N=100000

相反，即减少 10 倍。

>>> from math import sqrt 
>>> p = 0.8
>>> sqrt(p*(1-p)*10000)
40.0
>>> sqrt(p*(1-p)*100000)
126.49110640673517
于 2012-03-04T22:11:56.967 回答